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SUMMARY 

The  U.S.  Army  Mobility  Equipment  Research  and  Development  Command 
(MERADCOM)  has  been  involved  in  a program  for  development  of  an  off-road 
mine  detection  system.  The  Georgia  Institute  of  Technology  Engineering 
Experiment  Station  (GIT/EES)  was  awarded  a nine-month  contract  to  gather  and 
analyze  wideband  radar  spectral  data  to  determine  detection  and  classifica- 
tion potential  for  buried  objects  of  interest.  In  particular,  the  portion 
of  the  spectrum  which  is  most  useful  in  the  detection,  discrimination,  and 
classification  of  buried  targets  was  to  be  identified.  Later,  the  contract 
was  modified  to  include  additional  measurement  and  analysis  tasks  and  to 
extend  the  effort  to  18  months. 

Two  sets  of  data  were  taken.  Wideband  CW  radar  data  were  collected 
by  personnel  of  the  National  Bureau  of  Standards  (NBS)  in  Boulder,  Colorado, 
using  a computer-driven  automatic  network  analyzer.  For  the  initial  measure- 
ments used  in  this  program,  a pair  of  broadband,  non-dispersive  antennas 
was  mounted  on  a cart  and  pointed  at  the  ground.  Wideband  pulse  radar 
data  were  collected  by  the  Georgia  Institute  of  Technology  Engineering  Experi- 
ment Station  using  a Government  Furnished  Equipment  (GFE)  short-pulse  radar 
owned  by  MERADCOM.  The  wideband  target  signatures  obtained  using  the  MERADCOM 
short-pulse  radar  were  used  to  verify  the  utility  of  the  processing  algorithms 
which  were  developed  from  wideband  CW  radar  data.  The  MERADCOM  radar 
is  mounted  on  a General  Electric  Electro-Trak  tractor  with  a broadband,  non- 
dispersive  antenna  which  is  adjustable  in  height.  The  tractor-mounted  radar 
can  be  easily  moved  from  one  location  to  another  to  accommodate  measurements. 
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The  measurement  and  analysis  tasks  were  completed,  but  mixed  success 
was  achieved  in  finding  exploitable  target  discriminants.  The  performance 
of  each  processing  step — detection,  discrimination,  and  classification — 
was  acceptable  for  the  specified  scenario  when  the  soil  was  relatively  dry 
(i.e.,  zero  to  ten  percent  moisture  content).  However,  when  the  soil  was 
more  moist  than  this  (e.g.,  ten  to  thirty  percent),  the  performance  of  the 
processor  deteriorated  noticeably. 

The  processing  time  required  with  a modest  amount  of  processing  equip- 
ment was  determined  to  be  about  40  milliseconds  for  each  mine  which  was 
correctly  detected,  discriminated,  and  classified. 

It  is  suggested  that  only  the  spectrum  from  0.5  to  1.5  GHz  needs  to 
be  used  to  adequately  detect,  discriminate,  and  classify  anti-vehicular, 
non-metallic  mines  in  relatively  dry  soil.  Thus,  an  antenna  smaller  and 
lighter  than  the  ones  specified  with  the  MERADCOM  tractor-mounted  radar  can 
be  utilized  providing  a more  portable  mine  detection  system. 
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PREFACE 

The  research  on  this  program  was  conducted  by  the  Radar  Analysis  Divi- 
sion, Radar  and  Instrumentation  Laboratory,  Engineering  Experiment  Station, 
Georgia  Institute  of  Technology,  Atlanta,  Georgia,  with  Mr.  James  A.  Scheer 
serving  as  Project  Director,  and  Dr.  J.  D.  Echard  serving  as  Associate  Project 
Director. 

This  report  was  prepared  by  the  Engineering  Experiment  Station  at  the 
Georgia  Institute  of  Technology  for  the  U.S.  Army  Mobility  Equipment  Research 
and  Development  Command  under  Contract  DAAG53-76-C-0112.  The  Technical  Monitor 
was  Mr.  L.  Mittelman,  Jr.  For  the  purposes  of  internal  control  at  Georgia 
Tech,  the  effort  was  designated  project  A-1828.  The  measurements  reported 
herein  were  conducted  by  Mr.  G.  Counas  under  the  direction  of  Mr.  F.  Clague 
of  the  U.S.  Department  of  Commerce,  National  Bureau  of  Standards. 

The  final  report  summarizes'  the  work  performed  under  this  contract  and 
gives  recommendations  for  a follow-on  program.  It  is  divided  into  two  volumes. 
Most  of  the  work  is  presented  in  Volume  I;  Volume  II  contains  a brief  descrip- 
tion of  the  inert  mines  used  in  the  data  collection  phase  and  identifies  the 
type  A,  B,  and  C mines  referred  to  in  Volume  I.  This  report  covers  work 
which  was  performed  from  March  1976  through  September  1977.  In  addition  to 
this  final  report,  an  interim  report  was  Issued  in  December  1976  entitled 
"Spectral  Characteristics  of  Buried  Objects." 

The  authors  are  grateful  to  Messrs.  Clague  and  Counas  of  the  National 
Bureau  of  Standards  for  their  care  and  perserverance  in  conducting  the  mea- 
surements, often  under  difficult  conditions,  and  to  Mr.  Clague  for  supplying 
the  photographs  reproduced  herein. 
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The  guidance  of  Mr.  J.  L.  Eaves  and  Dr.  E.  K.  Reedy,  Associate 
Director  and  Director,  respectively,  of  the  Radar  and  Instrumentation 
Laboratory  Is  appreciated.  We  also  appreciate  the  assistance  of  Mr.  M.  C. 
McGee  In  proofreading  these  two  volumes  and  suggesting  Improvements  to  the 
presentation  of  this  material.  Our  final  appreciation  goes  to  the  Radar 
Analysis  Division  secretary,  Mrs.  S.  A.  Busby,  who  typed  this  report. 
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SECTION  I 
INTRODUCTION 

The  U.S.  Array  Mobility  Equipraent  Research  and  Developraent  Command 
(MERADCOM)  has  been  involved  in  a program  for  development  of  an  off-road 
mine  detection  system.  In  March  of  1976,  the  Georgia  Institute  of  Technology 
Engineering  Experiment  Station  (GIT/EES)  began  a nine-month  contract  to  gather 
and  analyze  wideband  radar  spectral  data  to  determine  detection  and  classifi- 
cation potential  for  buried  objects  of  interest.  In  particular,  the  portion 
of  the  spectrum  which  is  most  useful  in  the  detection,  discrimination,  and 
classification  of  buried  targets  was  to  be  Identified.  In  January  of  1977, 
the  contract  was  modified  to  include  additional  measurement  and  analysis 
tasks  and  to  extend  the  effort  to  18  months.  The  measurement  and  analysis 
tasks  were  completed,  but  mixed  success  was  achieved  in  finding  exploitable 
target  discriminants.  This  final  report  describes  the  work  accomplished  under 
contract  DAAG53-76-C-0112 . 

Two  sets  of  data  were  taken.  The  wideband  radar  data  were  collected 
personnel  of  the  National  Bureau  of  Standards  (NBS)  in  Boulder,  Colorado, 
using  a computer-driven  automatic  network  analyzer.  The  equipment  is  housed 
in  a mobile  van  and  can  therefore  be  driven  to  any  site  accessible  by  vehicle; 
a trailer-mounted  generator  set  augments  this  mobility.  This  equipment  in- 
cludes a rubber-tired  cart  which  can  be  commanded  by  the  computer  to  advance 
to  any  of  a set  of  predetermined  positions  and  to  then  remain  stationary  while 
the  data  are  collected.  For  the  initial  measurements  used  in  this  program,  a 
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pair  of  broadband,  non-dispersive  antennas  was  mounted  on  the  cart  and  pointed 
at  the  ground. 

Wideband  pulse  radar  data  were  collected  by  the  GIT/EES  using  a Govern- 
ment Furnished  Equipment  (GFE)  short-pulse  radar  owned  by  MERADCOM.  The  wide- 
band target  signatures  obtained  using  the  MERADCOM  short-pulse  radar  were  used 
to  verify  the  utility  of  the  processing  algorithms  which  were  developed  from 
wideband  CW  radar  data.  The  MERADCOM  radar  is  mounted  on  a General 
Electric  Electro-Trak  tractor  with  a broadband,  non-dispersive  antenna  which 
is  adjustable  in  height.  The  tractor-mounted  radar  can  be  easily  moved  from 
one  location  to  another  to  accommodate  measurements. 

GIT/EES  personnel  analyzed  the  data  to  identify  distinguishing  character- 
istics in  the  temporal  and  spectral  signatures  of  the  targets.  These  distin- 
guishing characteristics  were  exploited  to  develop  detection,  discrimination, 
and  classification  algorithms.  The  Fast  Fourier  Transform  was  used  for  trans- 
forming the  data  into  the  time  domain  or  the  frequency  domain,  as  required. 

The  analysis  also  included  direct  comparison  of  measured  spectral  and  temporal 
data  with  data  derived  using  a mathematical  model . 

Techniques  for  detecting  the  presence  of  potential  targets  buried  in 
soil  of  various  types  and  moistures  were  investigated.  Discrimination  algo- 
rithms to  differentiate  between  mines  and  non-mine  objects  were  employed  and 
evaluated.  In  addition,  classif loation  techniques  were  used  to  distinguish 
the  various  types  of  mines.  Each  of  these  processes  were  evaluated,  and  per- 
formance curves  were  generated  and  are  included  in  this  report  with  accompany- 
ing discussions. 
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The  final  report  summarizes  the  work  performed  under  this  contract  and 

i 

gives  recommendations  for  a follow-on  program.  It  is  divided  into  two  volumes. 
Most  of  the  work  is  presented  in  Volume  I;  Volume  II  contains  a brief  descrip- 
tion of  the  mines  used  in  the  data  collection  phase  and  identifies  the  type 
A,  B,  and  C mines  referred  to  in  Volume  I. 

Section  II  of  Volume  I details  the  data  collection  effort  with  both  the 
NBS  measuring  equipment  and  the  MERADCOM  short-pulse  radar.  The  objects  used 
in  the  measurement  program  were  two  types  of  inert  non-metallic  mines,  a 
root,  two  rocks,  and  a flat  metal  plate,  all  buried  in  various  types  of  soil 
with  various  degrees  of  moisture  content. 

I Since  one  of  the  objectives  for  the  mine  detection  program  was  to  develop 

algorithms  to  automatically  process  real-time  daca,  it  was  necessary  to  visu- 
ally observe  and  compare  the  temporal  and  spectral  responses  for  the  various 
objects  sensed  with  the  measurement  system.  The  data  were  processed  off-line 
with  Georgia  Tech's  CYBER  74  computer  to  provide  plots  of  time  and  frequency 
responses.  The  processing  equipment  required  for  this  purpose  is  discussed 
:1  in  Section  III. 

When  the  measured  data  were  processed  and  examined,  it  was  noted  that 
radar  returns  from  the  same  target  at  the  same  depth  varied  with  soil  moisture. 
Since  consistency  in  target  return  signatures  seemed  lacking  in  the  measured 
i data,  a generalized  mathematical  model  was  developed  for  the  metal  plate  and 

[ the  type  A,  B,  and  C targets  buried  in  soil.  This  model  provided  an  under- 

■■  standing  of  the  target  and  the  environment.  The  mathematical  model  was  pro- 

grammed on  a general  purpose  computer  to  yield  plots  of  the  temporal  and 
spectral  characteristics  of  radar  returns  from  various  types  of  targets. 
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Parameters  such  as  target  depth  and  soil  moisture  were  varied,  and  changes  in 

target  return  characteristics  were  observed  and  compared.  A discussion  of  * 

the  mathematical  model  and  its  describing  equations  are  presented  in  Section  IV. 

Section  V discusses  the  selection  and  development  of  the  discrimination/ 
classification  algorithms.  Section  VI  presents  the  algorithms'  performance 
and  compares  the  various  techniques  that  were  considered.  The  real-time  hard- 
ware requirements  for  the  recommended  discrimination/ classification  techniques 
are  discussed  in  Section  VII.  Finally,  the  conclusions  and  recommendations 
are  given  in  Section  VIII. 

Three  appendices  present  plots  of  the  collected  data  and  plots  of  the 
performance  for  each  processing  algorithm  considered.  A brief  effort  was  , 

expended  on  the  homomorphic  deconvolution  approach  and  this  is  reported  in 
Appendix  D. 
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SECTION  II 

I' 

DATA  COLLECTION 

In  the  following  paragraphs,  the  two  measurement  systems  used  in  ob- 
taining the  radar  data  are  briefly  described,  along  with  some  of  the  measure- 
ments. The  objects  used  in  the  measurement  program  were  two  types  of  inert 
non-metallic  mines,  a root,  two  rocks,  and  a flat  metal  plate,  all  buried  in 
various  types  of  soil  with  various  degrees  of  moisture  content. 

The  first  task  of  the  subsurface  target  data  collection  effort  was 
to  decide  what  and  how  many  measurements  were  to  be  made.  The  decision  was 
influenced  in  part  by  the  equipment  available,  the  nature  of  the  obstacles 
of  interest,  and  the  range  of  soil  characteristics  desired;  the  amount  of  time 
' and  funds  that  were  available  for  the  task  also  influenced  the  decision. 

Two  types  of  GFE  measuring  equipment  were  made  available.  Initially, 
a Hewlett-Packard  computer-controlled  automatic  network  analyzer,  installed 
in  a 19-foot  step  van  and  assembled  in  the  S-parameter  configuration,  was 
used  to  obtain  reflectivity  measurements  from  various  types  of  objects  buried 


in  soil  of  various  textures.  This  equipment  was  owned  by  the  U.S.  Army  and 


operated  by  the  National  Bureau  of  Standards  in  Boulder,  Colorado.  Later  the 
Army  furnished  to  Georgia  Tech  a short-pulse  radar  (1.2  nanosecond  pulse 
length)  originally  developed  by  the  CALSPAN  Corporation.  It  was  used  as  a 
tool  to  obtain  wideband  target  signatures  for  use  in  verifying  the  utility 
of  the  processing  algorithms  which  have  been  developed  with  the  network 
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taken  from  the  two  measurement  systems  will  be  presented. 

NBS  MEASUREMENTS 

The  NBS-operated  equipment  used  in  collecting  the  data  will  first  be 
discussed  and  the  measurement  sites  described.  The  selection  of  the  test 
matrix  and  some  of  the  data  parameters  will  then  be  explained.  This  will  be 
followed  by  a discussion  of  some  of  the  salient  features  found  in  the  NBS  data. 
Measurement  System 

The  measurement  equipment  consists  of  a government-furnished  Hewlett- 
Packard  computer-controlled  automatic  network  analyzer  installed  in  a 19-foot 
step  van  for  mobility.  Only  its  salient  features  will  be  discussed.  The 
equipment  and  operation  thereof  will  be  described  in  more  derail  in  a manual 
currently  being  prepared  by  NBS  personnel.  [1]  The  data  were  collected  by 
means  of  an  automatic  network  analyzer  whose  basic  function  is  to  measure 
the  scattering  (S)  parameter  matrix  of  a two-port  network  (or  device)  at  a 
sequence  of  discrete  frequencies.  The  S-matrix  is  a set  of  four  complex 
numbers  that  describe  the  network  behavior*  representing  the  reflection  and 
coupling  coefficients  relating  incident  and  reflected  vaves.  [2] 

Specifically,  we  may  represent  incident  voltage  or  electric  field 
waves  injected  into  the  two  ports  by  a^  and  as  shown  in  Figure  1,  and 
waves  emerging  from  these  ports  by  and  a^  and  a^  are  considered 

independent  quantities  applied  externally  (and  simultaneously,  if  necessary) 
while  b^  and  represent  the  resulting  response.  The  response  involves 


the  S-parameters  via  the  network  equations 


bi  - Sj_^  a^  + S^2  ^2 


(1) 


S21  + S22 


(2) 
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where  and  5^2  are  the  reflection  coefficients  associated  with  ports  1 
and  2,  while  Sj^2  ^ad  822^  are  port-to-port  transmission  coefficients.  In 
general,  it  requires  four  measurements  to  determine  the  four  coefficients, 
but  for  linear,  passive  devices  82^^  = S2^2»  only  three  measurements  are 
needed.  In  the  measurements  program,  identical  antennas  were  used;  hence, 
it  was  assumed  that  82,  was  essentially  the  same  as  and  only  two  com- 

plex quantities  measured  were,  and  52^^.  All  four  parameters  can  be 
determined  by  measuring  b^  and  b2  with  either  or  a2  set  to  zero. 


a? 


-<■ 


Figure  1.  Two-port  network. 


The  S-parameters  are  functions  of  frequency  and  the  automatic  network 
analyzer  can  be  indexed  through  discrete  frequencies  from  about  100  MHz 
through  4 GHz.  A minicomputer  with  16,000  words  of  memory  is  used  to  con- 
trol the  frequency  stepping,  as  well  as  to  position  the  antenna  cart  along 
a fixed  path  and  to  record  the  measured  Information  on  tape.  In  its  factory- 
delivered  configuration,  the  automatic  network  analyzer  was  restricted  to 
a rapacity  of  53  frequencies  due  to  software  stored  in  core,  but  NB8  personnel 
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were  later  able  to  expand  this  to  a total  of  128  frequencies.  As  will  be 

« 

discussed  later,  the  number  of  frequencies  must  be  an  Integral  power  of  2 In 
order  to  exploit  the  features  of  the  Fast  Fourier  Transform  (FFT) . 

The  device  connected  to  the  automatic  network  analyzer  consisted  of 
a pair  of  non-dlsperslve,  broadband  antennas  aimed  at  the  ground.  The  Input 
connector  of  one  antenna  became  port  #1  and  that  of  the  other  antenna  became 
port  #2.  Since  the  electromagnetic  fields  are  not  confined  within  a closed 
system,  the  environment  surrounding  the  antennas.  Including  antenna  coupling, 
cable  coupling,  the  ground,  and  any  nearby  burled  obstacles,  becomes  an  Inte- 
gral part  of  the  "two-port  device."  The  network  analyzer  does  not  know  this, 
of  course,  and  although  the  physical  arrangement  of  this  "network"  Is  com- 
plicated, the  network  analyzer  obligingly  measures  its  characteristics. 

The  antennas  were  connected  to  the  equipment  Inside  the  van  by  a pair  ♦ 

of  coaxial  cables,  each  about  50  feet  long  and  supported  by  a boom  mounted 
on  the  van.  The  system  had  to  be  periodically  calibrated  (every  half  hour 
or  so,  perhaps  longer  under  good  conditions)  by  the  substitution  of  standard 
devices  for  the  device  to  be  measured.  These  must  be  connected  at  the  far 
end  of  the  cables  and  take  the  form  of  precision  short  circuits,  matched 
loads,  feed-throughs  and  the  like.  The  cables  themselves  change  character- 
istics with  changing  temperature  and  NBS  personnel  have  found  that  wrapping 
them  with  metal  foil  reduces  their  sensitivity  to  wind  and  sunlight. 

Measurement  Sites 

The  measurements  were  conducted  on  35-foot  test  lanes  prepared  at  two 
different  sites  in  and  near  Boulder,  Colorado.  Three  soil  types  were 
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represented  and,  with  appropriate  test  lane  flooding,  at  least  three  levels 
of  moisture  content  were  achieved.  The  first  site  was  on  U.S.  Department 
of  Commerce  property  in  the  city  of  Boulder  and  the  second  was  on  Gun  Barrel 
Hill,  located  some  15  miles  east  and  north  of  Boulder. 

The  NBS  site  contained  two  types  of  soil,  one  being  a sandy  loam 
trucked  in  from  off  the  site  several  years  ago  to  build  up  a thick  layer — 
perhaps  4 feet — of  homogeneous  soil.  The  sandy  loam  constituted  the  first 
20  feet  of  the  test  lane.  The  remaining  15  feet  contained  soil  of  essentially 
the  same  texture  but  Interspersed  with  rocks  and  pebbles  whose  exact  sizes 
and  distributions  are  not  known  with  certainty.  Near  the  surface,  these 
rocks  were  a few  centimeters  in  diameter. 

The  soil  at  the  Gun  Barrel  Hill  site  was  a class  of  bentonite  inter- 
spersed with  small  pebbles  and  rocks.  The  site  had  been  leased  by  the  Depart- 
ment of  Commerce  from  a local  farmer  who  had  been  raising  crops  on  the  land. 

One  corner  of  a field  had  been  cleared  of  vegetation  for  an  unrelated  experi- 
ment and  NBS  personnel  laid  out  the  test  lane  on  a small  section  of  this 
corner. 

The  same  test  lane  layout  was  used  at  both  sites.  As  shown  in  the 
dir.gram  of  Figure  2,  seven  targets  consisting  of  both  natural  and 
man-made  objects  were  buried  60  inches  apart  along  a straight  line.  Pre- 
liminary background  runs  were  conducted  to  collect  reference  data  at  closely 
spaced  stations  along  two  separate  five-foot  stretches  of  the  test  lanes 
before  the  targets  were  burled.  After  the  measurement  program  got  underway, 
the  cart,  provided  with  electric  stepping  motors  controlled  by  the  minicomputer 
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Figure  2.  Test  lane  layout. 
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in  the  system,  was  commanded  to  stop  at  given  positions  over  or  near  each  tar- 
get. These  positions  were  -8,  -4,  0,  4,  and  8 inches  from  the  target,  and 
a sixth  position  midway  between  targets  (i.e,,  30  inches)  was  included  for 
the  purpose  of  acquiring  "background"  data.  At  each  position,  the  frequency 
was  indexed  through  123  discrete  values  at  which  the  S-parameters  were 
measured  and  recorded  automatically  on  magentic  tape.  A boardwalk  was  laid 
out  and  leveled  alongside  the  test  lane  to  provide  a smooth  runway  for  the 
antenna  cart.  The  cart  was  guided  along  its  path  by  an  aluminum  rail  pinned 
to  the  ground  on  one  side  of  the  boardwalk.  The  test  target  positions  were 
carefully  located  by  strings  stretched  along  and  across  the  target  center- 
line  track.  These  and  other  details  of  the  test  lanes  can  be  seen  in  the 
photographs  of  Figures  3 through  8. 

The  moisture  content  of  the  soil  was  usually,  but  not  always,  deter- 
mined by  two  different  means.  The  most  convenient  method  utilized  a water- 
carbide  reaction  inside  a sealed  pressure  vessel.  Known  amounts  of  soil  and 
the  carbide  reagent  are  deposited  in  the  container  but  not  allowed  to  contact 
each  other  until  the  vessel  has  been  sealed.  Then  the  two  are  mixed  by 
external  agitation  and  the  resulting  chemical  reaction  produces  acetylene, 
generating  heat  and  a rise  in  internal  pressure.  A pressure  gauge  calibrated 
in  per  cent  moisture  content  converts  the  pressure  reading  to  moisture  con- 
tent. The  "bomb"  technique  is  within  a percentage  point  in  accuracy,  it  is 
claimed,  and  oven-drying  techniques  tend  to  support  this  claim,  at  least 
for  soils  with  modest  surface  area  to  volume  ratios.  The  oven  method  relies 
on  a comparison  of  sample  weight  before  and  after  several  hours  of  baking 
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near  the  boiling  point  of  water.  It  is  considered  more  accurate  than  the 

bomb  technique.  | 

I 

Moisture  content  was  not  always  measured  as  a function  of  depth,  but  j 

in  those  cases  where  it  was,  it  tended  to  be  more  uniform  at  the  NBS  site  ; 

than  at  the  Gun  Barrel  Hill  site.  The  soil  at  the  NBS  site  was  lighter  than  | 

! 

at  Gun  Barrel  Hill  and  was  uniform  to  a depth  of  4 feet  or  more.  As  seen  | 

J 

in  Figure  7,  Gun  Barrel  Hill  soil  appears  to  be  thoroughly  laced  with  small  j 

rocks  and  pebbles  ranging  up  to  a few  centimeters  in  size;  however,  the  j 

last  15  feet  of  the  NBS  site  are  also  impregnated  with  assorted  small  rocks  | 

and  pebbles.  j 

Preliminary  Data  Measurements  1 

! 

When  the  contract  at  Georgia  Tech  began,  the  measurement  facility  at 
the  National  Bureau  of  Standards  had  already  been  developed  and  was  being 
used  for  similar  measurements.  Some  preliminary  data  collection  procedures 
had  already  been  established  and  GIT  was  immediately  in  a position  to  review 
some  of  the  data.  ’ 

The  very  first  data  consisted  of  only  23  samples,  starting  at  240  MHz 
and  going  up  to  2000  MHz,  with  80  MHz  increments.  It  was  apparent  that  this 
resolution  in  the  frequency  domain  was  too  coarse  to  derive  a good  data  base. 

We  therefore  suggested  that  the  frequency  resolution  be  improved.  The  next 
data  received  contained  48  points  from  120  MHz  to  2000  MHz,  giving  a 40  MHz  ■ 

resolution.  Considerably  more  detail  was  seen  in  this  spectral  data,  suggest-  j 

ing  that  we  improve  the  resolution  even  more  until  it  was  apparent  that  addl-  | 

tional  resolution  would  create  a diminishing  return.  Initially  there  appeared  j 


UNCLASSIFIED 

15 


UNCUSSiFIED 


to  be  a software  problem  in  storing  more  than  53  points;  however,  the  memory 
space  problem  was  resolved  and  a system  capability  of  128  points  was  achieved. 
The  next  set  of  data  was  taken  using  two  frequency  bands  per  target.  The 
lower  band  was  105  to  1995  MHz.  The  higher  band  was  2010  to  3915  MHz.  Both 
bands  had  a resolution  of  15  MHz.  These  data  gave  us  the  capability  to 
evaluate  the  deleterious  effect  of  degraded  resolution  by  eliminating  data 
points  during  computer  processing.  It  was  apparent  that  the  15  MHz  resolution 
was  fully  sufficient  and  in  fact  using  every  other  point  (simulating  30  MHz 
spacing)  provided  sufficient  resolution.  At  that  time,  it  was  determined 
that  data  could  be  recorded  from  about  100  MHz  to  about  3 GHz  without 
requiring  recalibration  of  the  measurement  system.  These  observations 
led  to  the  final  choice  of  measurement  range  from  120  MHz  to  3048  MHz  in 
24  MHz  steps. 

In  addition  to  conventional  mine  data,  it  was  suggested  that  some 
other  special  situations  would  be  required.  One  set  of  data  was  measured 
for  the  case  in  which  a mine  was  unearthed,  and  the  hole  was  refilled  with 
dirt.  This  data  is  useful  in  determining  the  difference  between  the  undis- 
turbed dirt  and  the  dirt  packed  on  top  of  the  various  targets. 

Another  special  run  was  made  with  the  bistatic  pair  of  antenna  aimed 
at  the  sky,  rather  than  the  ground,  to  obtain  some  information  on  the  undis- 
turbed coupling  between  the  antennas.  While  the  antennas  were  elevated,  a 
mine  was  positioned  8 inches  above  them  for  one  run  to  determine  the  character- 
istics of  one  of  the  targets  of  interest  without  the  effect  of  clutter. 
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Selection  of  the  Test  Matrix 

Early  in  the  project  it  was  determined  that  a test  matrix  had  to  be 
developed,  representing  the  method  for  collection  of  all  subsequent  data. 

Early  in  May,  some  preliminary  data  was  received,  having  only  23  points  with 
a frequency  spacing  of  80  MHz.  It  was  soon  decided  that  the  spectral  reso- 
lution was  not  adequate  for  picking  up  the  high  rate  of  oscillation  in  the 
spectrum.  Subsequently,  data  was  received  containing  measurements  from 
105  MHz  to  3915  MHz,  in  15  MHz  increments.  This  data  certainly  had  enough 
resolution,  as  determined  by  evaluating  the  degradation  when  every  other 
data  point  was  deleted,  providing  30  MHz  resolution.  This  data  also  demon- 
strated that  there  was  very  little  information  beyond  3000  MHz  because  of 
the  loss  in  the  antennas  and  the  loss  in  the  ground  at  the  higher  frequencies . 
Consequently  a measurement  program  was  defined  in  which  the  upper  frequency 
was  held  at  3048  MHz  and  the  lower  at  120  MHz,  with  the  measurements  made 
every  24  MHz.  This  provided  precisely  128  frequencies,  including  those 
below  120  MHz  which  could  not  be  measured  and  had  to  be  artificially  assigned 
some  value,  usually  zero.  The  measurements  could  be  carried  out  using  a 
single  calibration  run,  thereby  saving  time,  since  frequencies  above  (about) 
3000  MHz  would  have  required  a second  calibration. 

There  remained  the  problem  of  selecting  an  optimum  antenna  height  for 
the  measurements.  It  had  been  pointed  out  previously  that,  due  to  refractive 
effects  at  the  ground  surface,  it  may  be  possible  to  choose  the  antenna 
height  so  as  to  maximize  the  signal  return  due  to  a specular  reflection  from 
the  target.  Since  the  refraction  depends  on  soil  moisture  content,  the 
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optimum  antenna  height  would  be  different  for  different  moisture  conditions. 
NBS  personnel  therefore  carried  out  a sequence  of  measurements  in  which  the 
antenna  was  placed  2,  4,  6,  and  12  inches  above  the  ground  for  both  dry  and 
damp  conditions. 

The  "dry"  soil  moisture  content  was  reported  to  be  of  the  order  of 
2%  and  the  "damp"  about  13%.  The  measurements  were  conducted  in  the  NBS 
indoor  test  lane  and  the  damp  condition  was  created  by  flooding  the  test 
lane  with  several  inches  of  standing  water.  The  water  was  permitted  to 
percolate  through  the  ground  for  several  days  before  the  measurements  were 
made.  The  dry  soil  tests  showed  that  a 12-inch  antenna  height  was  not  as 
desirable  as  lower  heights,  but  that  the  differences  between  the  2-,  4-, 
and  6-inch  heights  were  not  pronounced.  Nevertheless,  the  target  seemed  to 
stand  out  a little  more  clearly  for  the  4-inch  height  than  the  others,  and 
this  was  selected  as  the  optimum.  The  damp  soil  tests  were  inconclusive, 
with  no  height  appearing  to  be  better  than  any  other.  Consequently,  GIT 
recommended  to  NBS  personnel  that  the  antenna  be  fixed  at  4 inches  above 
the  ground  for  all  tests. 

Upon  completion  of  this  evaluation.  Table  1 was  generated  as  a guide 
to  direct  the  measurements  effort  at  both  the  test  sites.  The  selection  of 
these  parameters  has  proven  to  be  acceptable  in  the  analysis  of  the 
data.  Table  2 shows  the  matrix  of  data  actually  collected. 
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Antenna  polarization: 
Antenna  orientation: 
Antenna  height : 
Targets : 

Positions  per  target: 
Target  spacing 
Depths : 

Moisture  contents: 
Frequencies : 


Table  1 

MEASUREMENT  MATRIX 


parallel  to  direction  of  travel 

bistatic 

4" 

3 mines,  3 false  alarms,  1 metal  plate 
-8",  -4'^0",  4",  8^t  30#  ♦ 

60" 

3",  6" 

low,  medium,  high 

120  MHz  to  3048  MHz,  24  MHz  steps 


Table  2 

ACTUAL  DATA  MATRIX 

Antenna  polarization: 

parallel  to  direction  of  travel 

Antenna  orientation: 

bistatic 

Antenna  height : 

4"  nominal 

Targets : 

mine  types  A,  B,  and  C;  (2) 
(1)  root,  (1)  metal  plate 

rocks ; 

1 Positions: 

-8",  -4",  0",  4",  8",  30" 

1 Target  spacing: 

60" 

Depths: 

3",  6" 

Soil  types/moisture  content 

: Loam:  7%,  17%,  18% 

Bentonite:  12-20%,  13-16%, 

26-30%,  30% 

i Frequencies: 

120  MHz  to  3048  MHz,  24  MHz 

steps 

Results  of  a typical  measurement  are  presented  for  a type  A mine,  a 
metal  plate  and  a rock  of  about  the  same  size  as  the  mine.  At  each  of  the 
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123  frequencies,  the  real  and  Imaginary  parts  of  the  parameter  (the 
coupling  coefficient  from  Port  1 to  Port  2)  were  recorded,  and  the  magni- 
tude of  displayed.  The  phase  angle  at  each  frequency  was  also  cal- 

culated from  the  real  and  imaginary  components  and  used  in  the  FFT  computa- 
tions. A casual  examination  of  the  following  figures  shows  that  the  spectra 
are  complicated. 

Before  the  line  spectra  were  transformed  to  the  time  domain,  the 
spectral  data  were  weighted  with  a Chebychev  weighting  function  designed  to 
yield  low  (40  dB)  time  sidelobes.  The  weighting  function  is  plotted  in 
Figure  9 where  it  is  normalized  to  the  peak  amplitude  occurring  at  mid- 
frequency,  f^.  After  the  spectral  data  are  multiplied  by  this  weighting 
function,  the  data  are  ready  for  transformation  to  the  time  domain.  The 
standard  FFT  requires  that  the  number  of  samples  be  an  integral  power  of  2; 
samples  were  selected  for  128  frequencies,  including  zero  frequency  (DC). 

Since  the  five  lowest  frequency  samples  lie  below  the  capability  of  the 
equipment  and  were  not  measured,  zeroes  were  arbitrarily  inserted  to  represent 
the  unmeasured  values. 

The  FFT  algorithm  used  in  the  generation  of  the  time  domain  signals 
creates  an  array  of  256  datum  points,  128  of  these  being  the  measured  data 
for  positive  frequencies  and  the  other  128  being  their  complex  conjugates 
for  the  corresponding  negative  frequencies.  The  total  time  record  length 
is  the  reciprocal  of  the  frequency  interval;  hence,  a 24  MHz  sampling  interval 
generates  a time  record  41-2/3  nsec  long.  Virtually  all  the  signals  of  inter- 
est occur  in  the  first  half  of  this  record,  however;  hence,  only  the  data  for 
the  first  20-5/6  nsec  are  presented. 
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Line  spectra  data  such  as  those  in  Figures  10  through  12  were  used 
to  synthesize  the  performance  of  a short  pulse  system  by  means  of  the  FFT.  » 

After  weighting  the  data  with  a Dolph-Chebychev  function  centered  on  the  mid- 
band frequencies  to  simulate  the  spectrum  of  a short  transmitted  pulse,  time 
histories  such  as  those  in  Figures  13  and  14  are  generated.  These  two 
figures  represent  the  short  pulse  response  of  a metal  plate  buried  6 inches 
deep  in  relatively  damp  soil  and  of  a type  A mine  in  relatively  dry  soil; 
the  amplitude  scale  of  both  time  histories  has  been  normalized  with  respect 
to  the  largest  amplitude  occurring  in  the  time  record. 

The  first  small  anomaly  near  the  left  side  of  the  time  signature  is 
apparently  a direct  leakage  of  energy  between  the  two  antenna  connectors, 
which  are  relatively  close  to  each  other.  The  large  return  near  the  center 
of  the  diagrams  is  predomlnar^^  the  ground  return,  however  there  is  an 

% 

additional  contribution  due  to  coupling  between  the  antenna  apertures.  The 
target  return  appears  a short  time  later,  after  the  "ringing"  from  the  ground 
return  has  subsided. 

A casual  inspection  of  the  target  returns  in  Figures  13  and  14  shows 
that  the  two  are  quite  different.  In  fact,  the  three-peak  characteristic  in 
Figure  14  was  recognizable  in  most  of  the  time  histories  obtained  for 
the  type  A mine.  The  spectral  characteristics  of  the  target  alone  can  be 
extracted  from  the  time  histories  by  applying  a range  gate  or  time  "window" 
that  deletes  all  the  returns  except  those  from  the  target.  Typically  the 
window  need  be  only  a few  nanoseconds  wide,  and  the  time  data  can  be  trans- 


! 

formed  back  to  the  frequency  domain  using  an  inverse  FFT.  This  process  gene- 
rates potentially  useful  spectral  data,  as  shown  in  Figures  15  and  16. 
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The  metal  plate  response  is  shown  in  Figure  15  and  that  of  a type  A 
mine  in  Figure  16,  and  the  difference  between  the  two  is  far  more  pronounced 
in  the  frequency  domain  than  in  the  time  domain.  Whereas  the  plate  exhibits 
a single  peak  (near  700  MHz),  the  type  A mine  response  contains  a sequence 
of  peaks  and  nulls  whose  frequency  separation  is  remarkably  constant  (800 
MHz  per  cycle  of  the  pattern) . This  Interference  pattern  is  attributable 
to  two  discrete  echo  sources  of  comparable  magnitude  going  in  and  out  of 
phase  as  the  frequency  is  swept.  The  target  is  transparent,  being  made  of 
essentially  lossless  material,  and  one  reflection  occurs  at  its  bottom 
surface  while  the  other  occurs  at  the  top. 

Spectra  such  as  those  in  Figures  15  and  16  were  generated  for  all 
seven  burled  objects  and  were  subjected  to  analysis  in  a search  for  discri- 
mination  features.  Although  this  procedure  is  straightforward  when  the 
target  return  is  separable  from  the  ground  reflection  in  the  time  domain, 
some  difficulty  was  encountered  for  shallow  targets  buried  in  relatively 
dry  soils.  This  is  illustrated  in  Figures  17  and  18.  The  time  signature 
of  a type  A mine  burled  chose  to  the  ground  surface  (3  inches)  in  dry  soil 
is  shown  in  Figure  17.  The  familiar  three- lobed  response,  as  in  Figure  14, 
is  still  present  but  the  initial  portion  of  its  pattern  is  altered  by  the 
trailing  edge  of  the  ground  return.  Since  the  ground  contribution  "spills 
over"  into  the  target  return,  a window  bracketing  the  target  now  removes 
only  part  of  the  ground  return.  Consequently  the  deconvolved  target  spectrum 
shown  in  Figure  18  is  significantly  different  from  the  type  A mine  character- 
istics shown  in  Figure  16.  Because  the  locations  of  peaks  and  nulls  of  these 
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Figure  17.  Time  Domain  Response  of  S 
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:ruin  of  Shallow  Target  in  Dry  Soil. 
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spectra  are  different,  a discrimination  scheme  based  on  Figure  16  would  not 
necessarily  classify  the  spectrum  of  Figure  18  as  being  that  of  a type  A 
mine. 

SHORT-PULSE  RADAR  MEASUREMENTS 

The  short-pulse  radar  system  used  for  collecting  the  short-pulse  data 
is  briefly  described  in  the  following  paragraphs.  In  addition,  the  condition 
of  the  site  where  the  data  was  taken  is  reviewed.  Then,  the  test  matrix  uti- 
lized is  explained  and  some  typical  data  collected  are  presented. 

Measurement  System 

The  experimental  vehicular-mounted  mine-detector  radar  system  supplied 
by  the  U.S.  Government  and  shown  in  Figures  19  and  20  primarily  comprises 
three  short-pulse  receiver/transmitters , three  wideband  antennas,  and  a 
storage  oscilloscope  mounted  on  a GFE  General  Electric  tractor  (Model  E-15) . 
Some  additional  equipment  was  supplied  but  was  not  used. 

The  radar  operates  , at  a pulse  repetition  frequency  of  5 MHz.  This 
pulse  repetition  rate  is  derived  by  dividing  down  a 15  MHz  crystal  controlled 
oscillator.  To  avoid  mutual  interference,  the  output  of  the  transmitter  pulses 
are  staggered  so  that  each  radar  operates  at  a different  time  interval.  This 
staggering  is  accomplished  by  shifting  the  phase  of  the  5 MHz  signal  to  provide 
0°,  120°,  and  240°  outputs.  A fourth  5 MHz  output  is  provided  which  is  divided 
in  frequency  by  a factor  of  500  and  used  for  synchronization  of  the  range 
gates  and  staircase  sweep  generator. 

Synchronization  of  the  radar  is  provided  by  a 5 MHz  sine-wave  wnich  is 
divided  down  by  a factor  of  500  to  produce  a 10  kHz  clock  signal.  This  clock 
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(1)  three  short-pulse  receiver /transmitter  modules 

(9)  oscilloscope  monitor 

(10)  audio  and  video  alarm  unit 

(11)  battery  box 

(12)  three  antenna  units 

(15)  recorder  remote  control  box 

Figure  19.  Vehicular-mounted  mine  detector  radar  (left  side  view). 
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(2)  radar  control  console 

(3)  recorder  junction  box 

(4)  display  generator 

(5)  magnetic  tape  recorder 

(6)  discriminator  unit 

(7)  storage  oscilloscope 

(8)  power  junction  box 
(14)  position  potentiometer 

Figure  20.  Vehicular-mounted  mine  detector  radar  (right  side  view). 


UNCLASSIFIED 

35 


UNCUSSIFIED 


signal  is  used  to  step  the  staircase  sweep  generator  which  moves  the  saippling 
probe  progressively  in  time.  The  received  signal  is  sampled  approximately 
50,000  times  for  each  10  ms  sweep.  By  using  the  sampling  receiver  technique, 
a real-time  waveform  of  approximately  10  nsec  duration  is  displayed  on  the 
monitor  scope  in  approximately  10  msec. 

Measurement  Sites 

All  of  the  sites  for  the  short-pulse  radar  measurements  were  selected 
at  random  in  the  grounds  surrounding  the  Georgia  Tech  Baker  Building.  The 
soil  was  dry  red  clay.  All  the  mines,  metal  plates,  and  rocks  were  buried 
with  their  top  surfaces  approximately  six  inches  beneath  the  soil  surface.  All 
large  rocks  and  other  debris  were  removed  during  the  burying  process.  In  most 
cases,  the  surface  grass  was  replaced  for  the  measurement.  Only  one  of  the 
three  radar  modules  was  used  in  data  collection;  the  same  radar  was  used  for 
all  measurements . The  radar  antenna  was  placed  approximately  four  Inches  above 
the  soil  surface  in  each  case. 

Selection  of  Test  Matrix 

A matrix  of  measurements  to  be  made  with  the  short-pulse  radar  was 
generated  and  is  given  in  Table  3.  Since  the  major  purpose  of  the  short-pulse 
measurements  was  to  verify  that  this  measurement  system  is  effective  for 
target  discrimination/classification,  all  of  the  measurements  were  made  in 
dry  soil . All  objects  were  buried  6 inches  below  the  soil  surface  in  a red 
clay  soil.  The  objects  buried  were  two  types  of  non-metallic  mines,  a metal 
plate  and  a rock  about  the  same  size  as  the  mines.  In  addition,  some  data 
were  taken  in  ground  without  buried  objects  and  is  labeled  in  Table  3 as 
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j background  data.  Each  object  was  buried  in  4 different  locations,  and  for 

i ( 

I each  location  three  radar  measurements  were  taken.  One  measurement  was 

taken  with  the  radar  antenna  centered  over  the  object,  and  one  measurement 
4 inches  on  either  side  of  the  center.  Thus,  a total  of  12  dry  measurements 
per  object  were  taken. 


Table  3 

SHORT-PULSE  MEASUREMENTS 


Type  A Mine  I 

3 measurements/object 

Type  C Mine  ( 

4 locations 

Metal  Plate  1 

Rock  1 

BACKGROUNDS 

3 measurements /location 

4 locations 

Typical  Data 

Some  typical  temporal  and  spectral  plots  of  the  data  taken  with  the 
short-pulse  radar  are  shown  in  Figure  21  and  Figures  22  and  23.  In  each 
figure,  six  plots  are  shown.  The  center  pair  of  plots  represent  the  time 
and  spectral  responses  of  measurements  taken  when  the  antenna  was  centered 
over  the  burled  object.  The  pair  of  plots  at  the  top  and  bottom  of  the 
figures  are  the  data  taken  when  the  antenna  was  centered  4"  on  either  side 
of  the  object  center.  The  time  plots  on  the  left  span  a time  of  zero  to 
seven  nanoseconds  and  the  spectral  plots  on  the  right  represent  zero  to  3.33  Ghz. 

i 

ft 
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0 Time  7ns  0 Frequency  3.3  GHz 

Figure  21.  MERADCOM  short-pulse  radar  returns  from  metal  plate  buried  6" 


below  ground  surface  in  dry  clay  soil. 
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The  plots  in  Figure  21  are  radar  returns  from  a metal  plate  buried  6" 
below  the  ground  surface  in  dry  soil.  The  signature  between  the  vertical 
lines  shown  in  the  time  plots  (outer  boundaries  of  the  range  gate)  approxi- 
mates the  transmitted  waveform  (inverted)  and  is  about  one  period  of  a sine- 
wave  with  center  frequency  of  approximately  1 GHz.  Thus,  the  period  of  the 
signature  is  about  one  nanosecond.  The  spectral  response  obtained  for  the 
gated  temporal  response  is  consistent  with  the  time  signature.  A calculated 
spectral  response  for  one  period  of  a perfect  sine  waveform  is  shown  is 
Figure  24.  Note  the  close  resemblance  of  the  measured  data  spectral  response 
to  the  calculated  one. 

The  time  response  preceedlng  the  range  gate  is  part  of  the  ground 
return.  The  sharp  peak  following  the  gate  is  the  edge  of  a 6.4  nanosecond 
gate  built  into  the  short-pulse  radar.  The  beginning  of  this  gate  is  at  the 
left  edge  of  each  figure.  It  shows  in  the  figures  because  the  plots  exceed 
6.4  nanoseconds  in  time  extent  (i.e.,  7.0  nanoseconds). 

Returns  from  a type  A non-metallic  mine  burled  6 Inches  below  the  soil 
are  shown  in  Figure  22.  Likewise,  mine  type  C returns  are  shown  in  Figure  23. 
In  each  case,  the  temporal  signature  of  interest  is  located  between  the  verti- 
cal lines  shown  in  the  figures.  The  corresponding  spectral  responses  are  also 
shown.  The  position  of  the  null  in  the  spectral  response  is  characteristic 
of  the  mines  being  sensed.  The  position  of  the  null  in  hertz  is  equal  to  the 
reciprocal  of  the  electrical  thickness  of  the  mine.  The  electrical  thickness 
is  related  to  physical  thickness  as  follows: 

t = ^ y/T-  (3) 

c ''  r 
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Figure  24.  Calculated  spectral  response  for  one  period  of  sine  waveform. 
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where 

d = physical  thickness 
c = velocity  of  light  in  free  space 
t = time 

= the  relative  dielectric  constant  of  the  material 
in  the  mine  (for  paraffin  wax,  = 2.25) 

The  type  A mine  is  about  4.5  inches  in  thickness.  Thus,  its  time 
thickness,  t,  is  about  1.14  nanoseconds.  The  null  should  be  at  1/t  or  0.88 
GHz.  Careful  measurement  of  the  spectral  plots  will  indicate  close  agreement. 
The  type  C mine  is  about  3.2  inches  thick;  therefore,  t is  0.81  nanoseconds. 
The  null  is  located  at  1.2  GHz  which  it  should  be. 
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SECTION  III 
DATA  REDUCTION 

One  of  the  objectives  of  the  mine  detection  program  is  to  develop 
algorithms  to  automatically  process  real-time  data  in  order  to  detect  and 
identify  non-metallic  mines  buried  in  soil.  However,  before  selecting  and 
testing  candidate  detection  and  classification  schemes,  it  was  necessary 
to  visually  observe  and  compare  the  temporal  and  spectra  responses  for 
the  various  buried  objects  sensed  with  the  measurement  system.  To  provide 
plots  of  the  time  and  frequency  responses,  the  data  were  processed  off-line 
with  Georgia  Tech's  CYBER  74  computer  as  indicated  in  Figures  25  and  26. 
Figure  25  outlines  the  processing  perforuicd  on  the  NBS  data  to  provide  the 
necessary  plots,  and  Figure  26  show'j  the  processing  necessary  to  generate 
plots  of  the  short-pulse  radar  data. 

As  indicated  previously,  the  NBS  data  were  collected  with  a network 
analyzer  which  is  a CW  device.  Each  CW  measurement  (amplitude  and  phase) 
represented  one  line  of  a line  spectrum.  To  simulate  a short-pulse  return, 
the  J28  lines  representing  frequencies  from  0 MHz  to  3.048  GHz  were 
transformed  from  the  frequency  domain  to  the  time  domain.  In  addition, 
amplitude  weighting  was  used  to  control  the  time  sldelobes , the  ground 
return  was  removed,  and  the  remaining  data  were  searched  for  the  proper 
range  gate  position.  An  inverse  Fast  Fourier  Transform  was  used  to  generate 
a "gated  spectrum" — the  spectrum  of  the  suspected  temporal  signature.  Plots 
of  the  time  and  frequency  signatures  were  generated  for  visual  inspection 
and  comparison. 
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Figure  25.  Processing  of  NBS  data  to  provide  temporal  and 
spectral  plots. 
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Figure  26. 


Processing  of  short-pulse  radar  data  to  provide 
temporal  and  spectral  plots . 
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The  short-pulse  radar  data  were  processed  in  a similar  manner  as  shown 
in  Figure  26  to  provide  the  necessary  plots;  however,  these  data  were  recorded 
as  a temporal  signature.  Therefore,  the  "forward"  Fourier  transform  was 
not  necessary. 

The  following  paragraphs  describe  in  detail  the  processing  required  to 
present  the  measured  data  in  a proper  visual  format  for  analysis . 

NBS  DATA 

Preliminary  data  were  recorded  and  stored  on  magnetic  tape  at  the  NBS 
test  site  in  Boulder,  Colorado  and  sent  to  Georgia  Tech.  The  data  on  these 
tapes  were  loaded  into  memory  in  the  Georgia  Tech  CDC  CYBER  74  computer  for 
processing  and  analysis.  Each  tape  consisted  of  several  data  runs,  each  con- 
taining scattering  parameter  frequency  response  measurements  of  a specific 

target.  These  data  tapes  contained  both  monostatic  measurements  and 

J.i  Q 

bistatic  measurements.  Each  data  run  contained  a heading  record  which  gave 
such  information  as  date,  antenna  and  target  position,  number  of  samples, 
and  name  of  the  target  measured.  Following  the  heading  were  the  actual  data 
samples,  keyed  by  frequency  and  with  samples  for  both  and  82^^  data. 
Following  the  samples  was  a delimiting  record  containing  two  values,  -1.0 
and  0,  to  identify  separation  between  data  runs.  In  response  to  the  need 
to  reduce,  process,  and  analyze  the  data  on  these  and  subsequent  tapes,  a 
series  of  analysis  computer  programs  was  developed. 

Computer  Program  Development 

The  processing  requirement  for  the  data  was  to  perform  an  inverse 
Fast  Fourier  Transform  (FFT)  on  a desired  spectral  data  run  and  display  the 
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resulting  time  response  data  in  a manner  that  could  easily  be  interpreted  by 
the  user.  The  FFT  program  read  either  monostatic  or  82^^  bistatic  data, 
as  specified  by  the  program  user,  from  the  data  file  and  data  run  of  the  user's 
choice.  The  program  was  designed  to  read  the  header  for  pertinent  information 
concerning  the  sampling  period  and  adjust  the  data  automatically  to  fit  the 
power-of-two  array  required  by  the  FFT.  The  output  was  in  the  form  of  a 
numerical  print-out  and  a plot  provided  by  the  CALCOMP  digital  plotter. 
Verification  of  the  inverse  FFT's  performance  was  obtained  by  entering  test 
cases  whose  temporal  responses  were  known.  To  provide  a faster  turnaround 
on  plots,  a plotting  routine  employing  a Tektronix  interactive  graphics  termi- 
nal was  added  to  the  program.  This  made  it  possible  for  the  user  to  employ 
the  program  in  an  on-line,  interactive  analysis  mode,  as  graphic  output  was 
now  returned  to  the  user  immediately.  It  was  also  desired  that  frequency 
response  data  be  displayed  graphically  so  a CALCOMP  routine  to  plot  signal 
magnitude  versus  frequency  was  added. 

Casual  observance  of  the  resultant  time  plots  indicated  that  all  useful 
information  involving  the  target  returns  in  the  time  response  array  was 
located  within  the  first  20  nanoseconds  of  time  samples.  To  display  only 
this  portion  of  the  array,  the  plotting  routines  were  modified  to  plot  only 
one-half  or  one-quarter  (depending  on  the  total  time  length)  of  the  time 
response,  rather  than  that  suggested  by  the  sampling  resolution. 

Another  addition  to  the  program  was  the  filtering  subroutine  which 
weighted  the  frequency  response  data  before  the  FFT.  Low  and  high  frequency 
limits  of  the  filter  could  either  be  set  by  the  user  during  program  execution 
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or  by  the  program  as  the  low  and  high  sample  frequencies  of  the  data  run 
being  processed.  Initially,  the  main  lobe  of  a sin  x/x  weighting  function 
was  used,  with  subsequent  addition  of  the  Dolph-Tchebychev  weighting  for 
sidelobe  control. 

At  this  point,  the  analysis  program  had  a variety  of  optional  opera- 
tional modes  based  on  selection  of  the  filtering  mode,  the  output  format, 
the  selection  of  appropriate  data,  and  others.  To  allow  selection  of  these 
various  options  in  an  efficient  manner,  a system  of  control  variables  was 
implemented.  Each  option  was  controlled  by  the  value  of  an  associated  inte- 
ger variable.  Values  were  assigned  to  these  variables  at  the  beginning  of 
the  program  in  one  of  two  ways.  In  the  first,  the  program  interrogated 
the  user  as  to  the  value  of  each  control  variable.  A message  explaining 
the  options  for  each  control  variable  was  printed  before  each.  In  the 
second,  predetermined  values  were  assigned  to  each  variable  as  part  of  the 
program  flow,  and  the  user  was  not  interrogated. 

Study  of  time  response  plots  for  data  revealed  a "ringing"  oscilla- 
tion. It  was  concluded  that  this  was  brought  about  by  a "spike"  at  105  MHz, 
detected  on  the  frequency  response  magnitude  plot  of  the  data.  To  suppress 
this  undesirable  ringing,  the  option  of  deleting  data  samples  below  300  MHz 
was  added  to  the  system,  made  possible  by  the  control  variable  system.  The 
low  frequency  spike  appeared  to  be  peculiar  to  the  data  and  the  low  fre- 
quency deletion  was  not  required  for  processing  of  82^^  data  associated  with 
the  bistatic  arrangement  of  the  measurement  system. 

Other  options  were  added  to  the  program,  mostly  pertaining  to  program 
control.  Normally,  the  program  would  read  and  process  a data  run  once,  but 
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when  specified  by  the  user,  the  program  could  reprocess  a previous  data  run, 
making  it  possible  to  study  the  effects  of  filtering  and  other  factors  in  an 
interactive  manner.  In  order  to  test  transform  and  data  accuracy,  the  option 
of  varying  the  size  of  the  array  used  in  the  FFT  routine  was  added. 

Data  collection  conforming  to  the  suggested  measurement  matrix  began 
and  more  comprehensive  data  tapes  began  arriving  at  Georgia  Tech.  Analysis 
of  this  data  suggested  new  requirements,  and  further  modifications  of  the 
program  became  desirable.  For  better  reference  and  organization,  the  two 
CALCOMP  plotting  routines  were  consolidated  to  plot  frequency  response  and 
time  response  data  side-by-side  on  a single  plot. 

The  range  gate  and  forward  FFT  option  was  implemented  to  aid  in  analy- 
sis of  data.  This  option  allowed  the  user  to  specify  two  "range  gate  limits" 
in  the  temporal  domain.  These  gates  were  entered  by  the  user  in  units  of 
nanoseconds,  which  the  program  converted  to  the  proper  array  Indices  for  the 
time  response  array.  The  program  then  generated  a gated  time  response  array, 
containing  only  the  time  response  data  located  within  the  range  gates.  Ele- 
ments of  the  array  outside  the  range  gates  were  assigned  the  value  of  zero. 

A forward  FFT  was  then  performed  on  the  gated  time  array  to  produce  a range- 
gated frequency  response  array.  This  modification  made  it  possible  for  the 
user  to  specify  a particular  portion  of  the  time  array,  such  as  that  associated 
with  a target,  and  exclude  all  other  Information  thus  producing  a frequency 
response  array  containing  only  the  desired  target  information.  To  display 
this  new  information,  the  CALCOMP  plotting  routine  was  designed  to  plot  fre- 
quency response,  time  response,  gated  time  response,  and  gated  frequency 
response  on  one  page. 
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As  data  from  the  NBS  measurements  site  arrived,  it  was  placed  in  mass 
storage  on  the  CYBER  74.  Although  the  data  contained  both  and  82^^  data, 
only  82^^  data  were  analyzed.  Each  data  run  contained  123  samples  between 
and  including  120  and  3048  MHz.  Because  of  the  accumulation  of  a large 
volume  of  data  and  this  method  of  storage,  disc  storage  costs  rose  consider- 
ably during  this  time.  Also,  program  input  was  impaired  by  the  large  volume 
of  data  and  the  necessary  A8CII  to  binary  conversion  during  data  read.  To 
cut  down  on  storage  and  processing  costs,  the  data  was  reformatted  and  con- 
verted to  binary  code.  All  redundant  and  unnecessary  information,  such  as 
the  8^j^  data  and  the  frequency,  was  removed.  The  FFT  program  was  modified 
to  read  this  reformatted  data,  and  all  files  in  use  were  reformatted. 

Analysis  of  the  Dolph-Tchebychev  filter  was  performed.  Time  side- 
lobes  of  -20  dB  to  -80  dB  were  tried  and  target  detectability  and  separation 
were  evaluated.  The  filter  parameters  were  subsequently  set  to  a -40  dB 
sidelobe  level.  This  filter  was  used  for  essentially  all  subsequent  pro- 
cessing and  analysis. 

To  augment  the  analysis  of  the  time  response  data  further,  the  Tektronix 
time  response  plotting  routine  was  modified  to  implement  the  terminal's  joy- 
stick cursor  as  a range  finding  device.  When  interrogated  by  the  computer, 
the  joystick  cursor  displays  a set  of  crosshairs  on  the  display  of  the  ter- 
minal, controlled  by  the  attached  joystick.  When  the  user  has  positioned 
the  crosshairs  in  the  desired  position,  he  returns  control  to  the  program  by 
typing  a character  on  the  keyboard.  This  returns  the  coordinates  of  the  cross- 
hairs to  the  program.  The  time  response  plotting  routine  employed  the  cursor 
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in  this  fashion.  After  the  time  response  was  plotted,  the  program  would  go 
into  the  cursor  mode.  To  determine  the  range  of  a given  feature  on  the  time 
response  plot,  the  user  would  position  the  crosshairs  over  that  feature  and 
request  the  time  associated  with  the  feature  and  the  cursor.  When  the  appro- 
priate command  was  given,  the  program  would  return  to  normal  processing. 

As  a result  of  analysis  using  the  Tektronix  cursor  routine  to  determine 
range  information,  it  became  desirable  that  certain  features  of  the  time 
response,  such  as  the  location  of  the  ground  return,  be  determined  automatic- 
ally by  the  program.  To  satisfy  this  need,  a subroutine  was  developed  to 
locate  the  ground  return.  Extensive  observation  of  the  temporal  plots  sug- 
gested that  the  nature  of  the  ground  return  was  consistent  and  could  be 
found  by  an  appropriate  algorithm.  The  subroutine  located  the  ground  return 
by  searching  for  the  first  two  positive  and  negative  peaks  associated  with 
the  ground  return,  which  have  consistent  and  recognizable  relative  amplitudes. 
This  marked  the  location  of  the  ground  return.  The  subroutine  then  located 
each  subsequent  crossing  of  the  zero  axis  of  amplitude,  and  stored  the  loca- 
tions of  these  crossings  in  an  array.  This  subroutine  became  the  basis  of 
later  development  of  an  automatic  range  gate  routine  for  location  of  the 
target  information  in  the  time  array. 

Another  CALCOMP  plotting  routine  was  added  for  generating  the  plots  for 
the  first  interim  report.  This  routine  plots  a series  of  5 data  runs  on  an 
8-1/2  by  11  inch  sheet  of  paper.  For  each  data  run,  there  were  three  adjacent 
plots,  each  2 by  1-1/2  inches.  The  first  plot  consisted  of  the  raw  frequency 
response  magnitude,  the  time  response  with  the  range  gates  marked  by  vertical 
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lines,  and  the  gated  frequency  response  magnitude.  Figures  found  In  Appendix 
A are  representative  of  the  format  in  that  report. 

Filter  Selection 

To  reduce  the  number  of  options  from  which  the  user  had  to  select,  it 

was  determined  that  selection  of  an  optimum  filtering  function  had  to  be  made 

Hanning,  Hamming,  and  Kaiser  window  functions  were  briefly  considered  along 
s in  X 

with  — - — and  various  Dolph-Tchebychev  filters;  however,  it  was  decided  that 
the  latter  two  represented  enough  variability  to  select  the  appropriate  func- 
tion. 

sin  X 

The  — — — function  was  used  in  the  initial  analysis  efforts.  The  sub- 
routine which  performed  the  actual  data  weighting  allowed  the  user  to  specify 
the  center  frequency  and  the  bandwidth  of  the  main  lobe.  The  use  of  a 
narrower  bandwidth  permitted  inclusion  of  more  of  the  sldelobes  within  the 
frequency  extent  of  the  data,  producing  lower  time  sldelobes  at  the  expense 
of  time  resolution. 


Later  the  Dolph-Tchebychev  weighting  function  was  implemented.  It  has 
1 ( 

the  form  1 , where  1^^  is  the  modified  Bessel  function  of  order  one  and 

X 

the  argument,  x,  is  defined  as  follows: 


(4) 


The  constant,  K,  is  related  to  the  desired  time  sidelobe  level  by 
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where  R is  the  square  root  of  the  reciprocal  of  the  sldelobe  level. 
Table  4 gives  the  value  of  K for  the  various  sldelobe  levels  studied. 


Table  4 

DOLPH-TCHEBYCHEV 

CONSTANTS 

Sldelobe  Level 
Magnitude 

dB 

R 

K 

.1 

-20 

10 

2.9932 

.01 

-40 

100 

5.2983 

.001 

-60 

1000 

7.601 

.001 

-80 

10,000 

9.9035 

An  example  of  a Dolph-Tchebychev  weighting  function  is  shown  in  Figure  9. 
Range  Gate  Algorithm  Development 

When  the  range  gated  temporal  response  was  initially  generated,  the 
position  of  the  range  gate  was  entered  by  the  computer  user  by  typing  in 
the  range  based  on  a priori  knowledge  of  the  target  position.  Subsequently, 
the  location  of  the  range  gate  was  determined  Interactively  by  using  the 
cursor  mode  in  the  Tektronix  display.  For  large  scale  processing  of  data, 
an  automatic  range  gate  placement  routine  is  desired.  Several  techniques 
for  accomplishing  this  function  were  explored  before  a satisfactory  algorithm 
was  generated. 

Figure  27  illustrates  some  of  the  general  conditions  which  impact 
the  problem  of  automatically  locating  the  area  of  target  return  in  a typical 
temporal  plot.  For  the  most  part,  the  time  resolution,  6t,  of  the  temporal 
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Figure  27.  Typical  temporal  return  from  a buried 
land  mine. 
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data  is  approximately  0.04  nanoseconds,  resulting  in  a time  delay  to  the  n^^ 
data  point  of  n6t  or  0.04  n in  nanoseconds.  The  first  data  of  interest  occur 
at  about  3 nanoseconds,  resulting  from  coupling  of  the  transmitted  signal 
from  the  transmit  antenna  coaxial  connector  to  that  on  the  receive  antenna. 
This  signal  component  does  not  affect  the  target  return  and  can  therefore  be 
disregarded.  The  next  significant  component  occurs  at  approximately  9 nano- 
seconds and  is  actually  a combination  of  the  coupling  from  the  transmit 
antenna  to  the  receive  antenna  and  the  reflection  from  the  surface  of  the 
ground.  This  region  exhibits  the  characteristic  positive-negative-positive 
sequence  which  is  located  by  the  analysis  program  and  identified  as  the 
ground  return. 

Since  the  ground  return  necessarily  precedes  the  return  from  a target 
buried  below  ground  level,  the  time  associated  with  the  ground  return  provides 
a lower  limit  on  the  position  of  target  information.  In  particular  the  ground 
return  can,  with  some  confidence,  be  associated  with  four  amplitude  peaks, 
alternating  in  sign,  or  equivalently,  five  zero-crossing  points.  The  end 
of  the  ground  return  can  thus  be  defined  to  be  the  time  t^  = 0.04  n^,  of  the 
first  data  point  occurring  after  the  fifth  zero-crossing  of  the  ground  return. 
The  time  t^  can  then  be  used  as  a lower  bound  when  determining  the  target 
near-range  gate. 

The  last  significant,  and  most  important,  feature  of  the  temporal  data 
shown  in  Figure  27  is  that  resulting  from  reflections  from  a buried  object. 

For  the  data  used  in  these  analyses,  no  target  was  buried  deeper  than  6 Inches 
below  the  surface  of  the  ground,  and  primary  target  reflections  could  not  be 
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seen  in  the  data  past  delays  of  approximately  16  nsec,  the  400th  data  point. 

Therefore,  the  search  for  near  and  far  time  gates  could  be  confined  to  the 

interval  (t  , 16  nsec) . 
z 

The  development  of  an  automatic  range  gate  technique  was  based  on  the 
supposition  that,  If  a target  Is  present  In  this  data.  It  should  manifest  It- 
self as  a recognizable  deviation  from  a background  level.  At  least  in  the 
drier  soils,  the  presence  of  a target  should  cause  the  amplitude  of  the 
reflected  signal  to  be  larger  than  this  background  by  some  amount.  The  first 
attempt  at  an  algorithm  for  the  search  for  a target  return  limit  was  effected 
by  determining  the  first  data  point  which  exceeded  a specified  threshold 
value  as  shown  in  Figure  28.  The  background  level,  y,  , , was  calculated  to  be 

DKgd 

the  mean  value  of  forty  data  points  past  the  time  of  any  known  target  return: 

, 439 

Joo 

where  |y^|  is  the  absolute  amplitude  of  the  n^^  data  point.  The  time  of  the 
far  time  gate,  was  calculated  as  the  time  of  that  data  point  y^  first 

encountered  such  that 


<>"„>  i '■J'bkgd' 


n = 400,  399,  398,  .. .,  n 


Likewise  the  time  of  the  near  time  gate,  t , was  calculated  as  the  time 

near 

of  that  data  point  y first  encountered  such  that 

n 
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If  no  such  points  could  be  located  which  satisfied  either  criterion,  default 
time  gate  limits  became  either  t^  or  16  nsec,  as  appropriate. 


Limits  of  Search 


V V 


16  nsec 


400  439 


Figure  28.  Non~adaptive  Threshold  Algorithm. 


When  this  automatic  gating  algorithm  was  implemented,  several  short- 
comings became  apparent.  Although  a multiplicative  factor  of  4 was  a viable 
parameter  in  dry  soil  conditions  most  of  the  time,  targets  did  not  always 
faithfully  exceed  the  background  level  by  the  required  amount,  particularly 
in  moist  soil.  If  this  margin  between  background  level  and  threshold  were 
lowered  to  accommodate  the  moist  soil  conditions,  false  alarms  occur,  causing 
the  placement  of  the  gate  to  stray. 

To  allevalte  this  problem,  the  basic  decision  mechanism  was  modified 
to  Include  provision  for  an  adaptive  threshold  level.  During  the  initial 
attempt  to  locate  the  time  gates,  a threshold  level  equal  to  4 times 
would  be  used,  where  the  background  level  y,  , , was  calculated  as  previously. 

DK^Q 

Rather  than  define  default  gate  values  if  no  target  decision  was  made,  as 
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had  been  done  before,  the  threshold  value  was  lowered  to  80  percent  of  Its 
previous  value  and  the  process  renewed.  The  second  search  criteria  became 

(y^)  > 3.2  n = 400,  399,  398,  . . . , n^  + 24  (9) 

for  the  far  time  gate  and 


(y  ) > 3.2  y,  , , 

■'n  — bkgd 


^^z+l’  ”z+2’  •••’ 


for  the  near  gate. 

If  the  second  search  did  not  result  in  a target  decision,  the  threshold 
was  again  lowered  to  80  percent  of  its  previous  value,  and  the  search  process 
was  repeated.  This  process  was  repeated  until  either  a target  decision 
was  made  or  the  threshold  finally  became  less  than  the  average  background 
level  I**  the  latter  case,  the  appropriate  default  time  gate  would 

be  used.  To  prevent  the  gates  from  becoming  too  close  to  one  another,  the 
far  gate  would  be  defined  first  and  the  near  gate  would  have  to  meet  the 
decision  criterion  and  occur  at  least  24  data  points  prior  to  the  far  gate. 
(Restricting  the  far  gate  to  the  interval  [n^  + 24,  400]  prevents  the  obvious 
contradiction  in  location  requirements.) 

The  implementation  of  such  an  adaptive  threshold  resulted  in  a much 
improved  target  detection  rate;  however,  it  was  noted  that  this  technique 
was  overly  sensitive  to  randomly  occurring  spikes  in  the  data  above  the 
threshold  level.  In  an  effort  to  alleviate  this  sensitivity,  the  adaptive 
threshold  scheme  was  modified  so  that  multiple  detections  above  the  threshold 
would  be  required  to  set  a time  gate. 
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Two  ways  of  accomplishing  this  requirement  were  investigated:  (1)  two 
successive  data  point  detections  for  gate  placement  (n  and  n-1) , and  (2)  two 
detections  out  of  three  successive  data  points  (n  and  n-1,  n and  n-2,  or 
n-1  and  n-2)  would  be  required.  The  performance  resulting  from  the  use  of 
two  successive  target  decisions  was  somewhat  superior  to  that  of  single  detec- 
tion, while  using  two  out  of  three  resulted  in  excessive  time  differences 
between  near  and  far  time  gates. 

Of  all  of  these  techniques  which  employed  straight  amplitude  thres- 
holds for  decisions  in  setting  gates,  that  which  employed  two  successive 
decisions  against  an  adaptive  threshold  was  the  most  effective.  However,  the 
level  of  performance  was  considered  to  be  far  from  that  which  could  be  ob- 
tained by  manually  setting  both  gates  upon  inspection  of  the  temporal  data. 

In  the  hope  that  a more  effective  automatic  technique  could  be  found,  a 
different  attack  on  the  problem  was  sought. 

To  desensitize  the  detection  algorithm  to  random  spikes  in  the  back- 
ground data,  sliding  window  averaging  was  utilized  on  the  absolute  value  of 
the  amplitude  versus  time  data.  This  type  of  averaging  accomplishes  two 
things:  (1)  it  suppresses  those  "spikes"  in  the  data  which  have  short  extent 

(in  time)  and  (2)  it  enhances  the  desired  signal  with  respect  to  an  average 
background  level. 


The  absolute  values  of  the  amplitude  data  (y^)  were  averaged  using 
sliding  window  widths  of  13  and  21  points  wide  in  separate  tests.  A new 
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for  the  21-wide  case  or  the  relationship 
in*n-t-6 
n)*n-6 

for  the  13-wide  case.  Note  that  the  averaging  was  undertaken  only  between 
the  appropriate  end  of  the  ground  return  and  the  maximum  range  of  Interest, 

16  nanoseconds. 

A typical  plot  of  the  smoothed  data  y is  shown  in  Figure  29 . Averag- 

m 

ing  over  21  data  points  proved  to  be  superior  to  that  over  13  in  terms  of 
amplitude  trends  for  target  location  identification.  The  "center"  of  the 
target  is  defined  to  be  that  data  point  at  time  corresponding  to  the 

largest  member  y(max)  of  the  smoothed  array.  Thus  the  largest  reflected  aver- 
age signal  is  construed  to  be  the  target. 
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Figure  29.  Smoothed  Data  Array  y. 
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The  temporal  extent  of  the  target  return  is  determined  by  locating  the 
closest  local  nulls  in  the  smoothed  array  y^  on  either  side  of  the  center 
point  at  time  By  observation,  these  nulls  indicate  the  actually 

observed  target  signal  boundaries  under  the  following  restrictions: 

(1)  The  values  of  the  nulls  of  the  smoothed  data  y(t  ) 

^ near 

and  must  be  less  than  .25y(max). 

(2)  Each  null  must  be  located  at  least  12  time  divisions 
(.48  nanoseconds)  from  the  center  of  the  target  return, 
so  that  the  entire  target  extent  must  be  at  leasr  24  time 
units  (.96  nanoseconds)  long. 

(3)  If  a null  value  of  the  smoothed  array  cannot  be  determined, 
the  time  gates  are  placed  at  the  end  of  the  ground  return 
t^  or  at  16  nanoseconds,  as  appropriate. 

The  operation  of  this  computerized  algorithm  appears  to  be  similar  to 
manually  selected  range  gates,  with  few  exceptions.  It  is  adequate  to 
produce  the  range  gates  for  analyzing  large  quantities  of  data,  but  there 
are  some  times  when  the  algorithm  places  the  range  gate  at  a less  than  optimal 
position. 

Elemoval  of  Ground  Return  from  Processed  Data  i 

The  range  gate  algorithm  just  described  sometimes  falls  to  work  properly 

because  of  the  size  of  the  ground  return,  '.^hen  the  soil  is  wet,  the  return  | 

( 

of  a buried  target  is  greatly  attenuated  as  compared  to  the  ground  return.  If 
a portion  of  the  ground  return  falls  within  the  time  window  being  searched  for 
the  buried  target  signature,  the  range  gate  algorithm  is  adversely  affected. 
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This  problem  is  particularly  severe  with  shallow  targets  burled  at  three- 
inch  depths  in  wet  or  moist  soil.  Therefore,  a technique  to  remove  the 
ground  return  before  searching  for  the  target  signature  was  considered  and 
is  discussed  in  the  following  paragraphs. 

When  targets  of  Interest  are  only  a few  Inches  below  the  soil  surface, 
the  target  return  and  the  ground  return  overlap  in  time.  To  facilitate 
detection  and  classification  of  target  returns,  it  is  desirable  to  remove 
the  overalpplng  ground  return.  One  way  to  accomplish  this  is  to  record  the 
return  from  a location  where  a target  is  not  present  and  subtract  it  from 
the  returns  obtained  in  locations  where  a target  of  Interest  may  be  present. 
Since  measured  data  exists  for  areas  with  and  without  targets,  this  scheme 
was  implemented  in  the  data  processing  program  which  has  been  previously 
described.  Since  the  radar  antenna  height  may  change  a little  with  position, 
the  computer  program  locates  Che  negative  peaks  of  the  two  ground  returns  to 
be  subtracted  and  aligns  them  by  the  proper  Indexing.  In  addition,  the  ampli- 
tude of  the  time  data  is  scaled  so  the  amplitude  of  the  two  negative  peaks 
are  equal.  Thus,  the  subtraction  of  the  two  time  responses  guarantees  a 
zero  value  where  the  largest  negative  peak  of  the  ground  return  existed. 

The  results  are  shown  in  Figures  30  and  31  for  two  selected  cases.  In 
Figure  30,  the  time  response  of  a metal  plate  burled  6 Inches  below  the  soil 
surface  with  and  without  the  ground  return  present  is  shown.  In  Figure  30(b) 
measured  data  GT5-65  is  subtracted  from  GT5-68.  Note  that  the  subtraction 
is  not  perfect.  The  residue  remaining  is  about  12%  of  the  original  peak 
amplitude  shown  in  Figure  30(a). 

UNCLASSIFIED 


64 


Amplitude  Amplitude 


UNCLASSIFIED 


r" 


I 

1 

0 


L. 

0 


! ! I 


Metal  Plate 


Ground  Return  || 

< ^ I'  II 


I 


!'  I 


Time 

(a) 


20.83  ns 


f I 


Metal  Plate 


Residue  of  GrounS^  ' 'j 
Return  Subtraction  | 


Time 

(b) 


20.83  ns 


Figure  30. 


Time  response  of  metal  plate  buried  6"  below  ground,  7% 
moisture  a.)  with  ground  return  present;  b.)  with  ground 
return  subtracted  (GT5-65  subtracted  from  GT5-68). 
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Figure  31.  Time  response  of  type  A mine  buried  6"  below  ground,  7% 
moisture  (a)  with  ground  return  present;  (b)  with  ground 
return  subtracted  (GT5-53  subtracted  from  GT5-56) . 
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In  Figure  31,  the  time  response  of  the  type  A mine  buried  6"  below  the 
ground  surface  with  and  without  the  ground  return  present  is  shown.  The  type 
A target  return  is  somewhat  smaller  than  the  ground  return  in  Figure  31(a). 

When  the  ground  return  is  subtracted  in  Figure  31(b),  the  remaining  time 
response  is  scaled  by  the  computer  program  to  fill  the  figure  window.  Thus, 
the  type  A target  return  appears  larger  in  Figure  31(b)  than  in  31(a).  If 
the  amplitudes  were  on  the  same  scale,  one  would  discover  that  the  peak 
amplitude  of  the  ground  subtraction  residue  in  31(b)  is  about  20%  of  the 
ground  return  peak  amplitude  in  Figure  31(a). 

The  results  obtained  in  these  two  cases  are  not  desirable  because  the 
residue  is  larger  than  most  target  returns  of  interest,  especially  in  wet 
soil.  Hence,  the  effect  of  the  ground  return  on  shallow  targets  may  still 
be  appreciable.  Ways  in  which  to  Improve  the  subtraction  were  sought.  It 
was  discovered  that  one  reason  the  ground  returns  were  subtracting  so  poorly 
was  the  position  of  the  time  samples  relative  to  the  ground  return  envelope. 

The  problem  is  illustrated  in  Figure  32  where  a typical  ground  return  envelope 
is  shown.  The  ground  return  is  an  analog  signal  sampled  in  the  computer  pro- 
gram processing  procedure.  In  the  present  processing  configuration,  the 

-12 

sampling  period  is  approximately  .041  nanoseconds  (41  x 10  ) seconds)  which 

corresponds  to  a one  way  range  of  0.5  inch  in  free  space.  Since  the  time 
position  of  the  samples  are  fixed,  if  the  antenna  height  above  ground  changes 
by  0.25  inches,  the  position  of  the  samples  relative  to  the  analog  envelope  ^ 

changes  by  one-half  sample  period  as  shown  in  Figure  32(b).  If,  for  example,  . 

I 

the  two  ground  returns  which  are  subtracted  are  each  sampled  as  illustrated 
in  Figure  32(a)  and  (b),  the  residue  will  be  sizeable.  Tho  computer  program  i 
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Figure  32,  Sampling  Position  Relative  to  the  Ground  Return  Analog 

Envelope.  (a)  Sampling  on  Peak  Value,  (b)  Sampling  off  * 

Peak  Value. 
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locates  the  largest  sample  value  In  each  case,  normalizes  the  two  sampled  ^ 

waveforms  to  have  the  same  maximum  value,  and  then  Indexes  one  sampled 

return  relative  to  the  other  to  align  the  two  maximum  samples.  Although  • 

I 

I 

I the  maximum  values  will  be  aligned,  the  envelopes  will  not  be.  The  dlffer- 

l 

I ence  in  sampled  values  near  the  zero  crossing  points  can  be  as  large  as  25% 

I 

of  the  peak  value. 

To  alleviate  this  problem  in  subtracting  the  ground  return,  the  measured 
data  with  the  ground  return  only  (no  target  present)  is  modified  to  reposition 
the  analog  envelope  relative  to  the  sampling  points.  This  is  accomplished  by 
introducing  the  appropriate  linear  phase  slope  in  the  measured  frequency  data. 

The  results  of  this  modification  to  the  subtracting  procedure  is  shown  in 
Figure  33  where  measured  data  GT5-71  is  subtracted  from  GT5-74  (type  B mine 
buried  6",  7%  moisture).  The  time  responses  of  these  two  records  are  shown 
in  Figure  34  where  they  have  been  aligned  in  time  and  overlayed  on  the  same 
plot.  Note  that  the  amplitude  scale  is  the  same  for  both  Figures  33  and  34. 

The  subtraction  residue  is  still  quite  large,  about  12%  of  the  ground 
return  peak  amplitude.  A close  observation  of  the  two  overlayed  waveforms 
in  Figure  34  indicates  why.  Note  that  the  trailing  edges  of  the  two  ground 
returns  are  significantly  different  and  are  the  source  of  the  largest  part 
of  the  residue  in  Figure  33.  Since  the  type  A mine  is  6"  below  the  ground 
surface  and  is  well  separated  in  time  from  the  ground  return,  this  difference 

* 

cannot  be  attributed  to  the  type  A target  time  response.  It  was  thought  that 
perhaps  the  difference  was  due  to  the  presence  of  a small  rock  near  the  sur- 
face of  measurement  GT5-71.  To  explore  this  possibility,  another  nearby 
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Figure  33.  Measured  Data.  Time  response  of  type  A mine  burled 
6"  below  ground,  7%  moisture — ^wlth  ground  return 
subtracted — GT5-71  subtracted  from  GT5-74. 
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ground  return  only  measurement,  GT5-77,  was  compared  to  GT5-71  by  plotting  ^ 

i 

both  curves  on  the  same  plot  as  given  in  Figure  35.  As  will  be  noted,  these  • j 

two  time  responses  track  each  other  very  closely.  The  only  remaining  explana- 
tion is  that  there  was  some  irregularity  near  the  soil  surface  in  measurement 
GT5-74. 

Although  the  reason  for  the  difference  in  the  two  waveforms  cannot 
conclusively  be  determined,  it  is  obvious  that  Irregularities  in  either  the 
ground  return  only  measurement  or  the  ground  return  and  target  measurement 
can  seriously  interfer  with  the  subtraction  of  the  ground  return.  Thus,  it 
was  concluded  that  the  ground  subtraction  approach  cannot  be  used  for  the 
consistent  removal  of  the  unwanted  ground  return. 

SHORT-PULSE  RADAR  DATA 

It  was  brought  to  our  attention  that  the  Army  was  in  possession  of  a 
short-pulse  radar,  originally  developed  by  the  CALSPAN  Corporation,  which 
could  possibly  be  utilized  to  test  and  evaluate  discrimination  techniques. 

Subsequently,  it  was  decided  by  the  Army  that  the  short-pulse  system  be  used 
as  a tool  by  Georgia  Tech/EES  to  obtain  wide-band  target  signatures  for  use 
in  verifying  the  utility  of  the  processing  algorithms  which  have  been  developed 
by  Georgia  Tech  with  the  Hewlett-Packard  network  analyzer.  The  availability 
of  this  system,  which  develops  a wide-band  transmitted  pulse  (approximately 
1.1  ns)  was  valuable  in  the  natural  follow-on  efforts  Involved  with  developing 
actual  signal  processing  algorithms  for  target  discrimination.  By  virtue  of 
the  wide  spectrum  of  the  transmitted  pulse,  this  system  represented  the  analog 
hardware  Implementation  of  the  FFT  technique  used  to  synthesize  time  domain 
data  in  the  previous  work.  With  this  system,  the  hardware  exists  which  allowed 
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Figure  35.  Measured  Data.  Time  responses  of  ground  without 
target — GT5-71  and  GT5-77  (ground  areas  adjacent 
to  type  C mine) . 
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us  to  obtain  time  domain  data  directly.  The  radar  was  used  to  derive  the 
signatures  which  were  recorded  on  magnetic  tape  and  read  into  the  CYBER  74 
computing  system  at  Georgia  Tech  for  synthetic  processing  as  illustrated  in 
Figure  26.  " The  outputs  from  the  radar  were  made  compatible  with  the  Georgia 
Tech  analysis  tools  and  computer  programs.  The  primary  objective  and  emphasis 
of  these  radar  measurements  were  to  verify  the  computerized  analysis  which 
was  previously  based  on  CW  measurements  over  a broad  range  of  frequencies. 

With  the  NBS  data,  a simulated  return  pulse  is  obtained  using  linear  Fourier 
analysis  techniques.  The  ground  and  buried  object,  however,  do  not  necessarily 
represent  a linear  system  and  realistic  measurements  are  required  for  compari- 
son purposes.  The  return  signature  was  obtained  by  using  the  GFE  short-pulse 
radar  and  temporal  and  spectral  analysis  was  applied  to  the  receiver  output 
for  comparison. 

The  three  basic  components  of  the  short-pulse  system  are  the  transmitter, 
the  receiver,  and  the  antenna  system.  The  transmitter  generates  a 1.1  nano- 
second monocycle  which  is  radiated  into  the  ground  from  a wide-band  horn 
antenna.  The  wide-band  receiver  signal  is  sampled*,  using  a technique  similar 
to  that  in  a sampling  oscilloscope,  at  a rate  compatible  with  conventional 
oscilloscope  displays,  conventional  laboratory  analyzers,  or  magnetic  tape 
recording  systems. 

In  order  to  collect  the  data  for  use  on  the  CYBER  74  computer,  the 
low  frequency  analog  signals  had  to  be  sampled,  stored  in  a random  access 

*The  effect  of  this  sampling  is  to  expand  the  time  response  by  10^:1.  Thus, 
the  1.1  nanosecond  pulse  appears  on  the  radar  output  as  a 1.1  millisecond 
analog  pulse. 
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memory  (RAM)  and  transported  from  out-of-doors  where  the  measurements  were 
taken  to  a magnetic  tape  recorder  Inside  an  adjacent  building.  The  digital 
data  was  then  placed  on  magnetic  tape  in  a format  so  as  to  be  compatible 
with  the  CYBER  74.  The  following  paragraphs  describe  some  of  the  hardware 
and  software  that  was  developed  to  provide  data  compatible  with  the  computer 
programs  previously  developed  for  the  NBS  data. 

Design  and  Construction  of  the  Event  Recorder 

The  design  and  purpose  of  the  portable  event  recorder  is  to  allow  the 
sampling,  digitizing,  and  storing  of  an  electrical  event  or  return,  usually 
measured  from  some  reference  time  or  marker.  In  our  application,  the 
reference  mark  was  the  range  gate  signal,  available  from  the  radar  control 
console.  The  event  recorder  sampled  the  gated  video  beginning  at  the  range 
gate  mark  at  a 75  KHz  rate,  converted  the  analog  video  signal  to  eight  bits 
to  form  a sign  plus  magnitude  binary  word.  This  sampling/digitizing  process 
continued  for  approximately  7.0  msec,  or  until  512  words  had  been  recorded 
in  the  RAM.  This  comprises  the  "WRITE"  mode  of  operation.  The  event  recorder 
may  then  be  switched  to  battery  power,  unplugged  and  transported  for  up  to 
30  minutes  without  jeopardizing  the  recorded  data. 

This  capability  greatly  simplifies  the  data-taking  process.  The  "READ" 
mode  is  used  to  transfer  the  data  from  memory  to  a digital  tape  recorder. 
Either  Internal  or  external  clock  and  TTL  or  CMOS  levels  up  to  15  VDC  may 
be  used  to  read  the  data. 

The  block  diagram  of  the  event  recorder  is  shown  in  Figure  36.  As 
indicated  earlier,  it  is  a digital  data  recorder  with  associated  clocking 
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Generatlnfi  Computer  Compatible  Data  Tapes 

After  data  is  stored  In  the  event  recorder,  there  are  three  steps 
required  to  convert  it  to  a CYBER  74  computer  file: 

(1)  Data  is  read  from  event  recorder  onto  a magnetic  tape. 

(2)  The  tape  is  taken  to  the  computer  center. 

(3)  The  tape  is  mounted  and  read  and  then  converted  to  a 
useable  form. 

A special  program  named  "INTAPE"  reads  a binary  file  and  separates 
the  data  in  that  file  into  8-bit  integers.  "INTAPE"  is  designed  to  read  a 
file  assigned  to  a 9-track  magnetic  tape.  "INTAPE"  will  read  the  entire 
input  file  and  determine  how  many  8-bit  integers  could  be  packed  into  that 
file.  The  program  then  calculates  the  greatest  power  of  2 that  is  less  than 
the  maximum  number  of  8-bit  integers.  The  resulting  power  of  2 is  assumed 
to  be  the  desired  number  of  8-bit  integers. 

"INTAPE"  is  used  to  read  a tape  which  contains  8-bit  bytes  recorded 
across  the  9 tracks  (one  track  is  for  parity) . As  the  computer  reads 
the  tape,  the  bytes  are  packed  into  central  memory  (CM)  with  60  bit  words 
as  shown  in  Figure  37 . 
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The  8-blt  bytes  fit  evenly  Into  consecutive  pairs  of  CM  words,  15 
bytes  per  word.  "INTAPE"  takes  a pair  of  CM  words  and  converts  them  to 
the  15  8-bit  bytes.  The  program  looks  at  the  first  word  and  extracts  the 
first  7 bytes  by  shifting  the  word  unitl  the  desired  byte  is  right-justified. 
The  CM  word  is  then  masked  or  "And"  with  a word  which  has  only  the  right 
8-bits  set.  To  obtain  the  8th  byte,  the  program  takes  the  last  4 bits  of 
the  first  CM  word  and  the  first  4 bits  of  the  second  CM  word.  These  two 
4-bit  words  are  "Or'ed",  then  right-justified.  Bytes  9 through  15  are  ob- 
tained by  a process  similar  to  the  one  used  to  obtain  bytes  1 through  7 . 

These  8-bit  bytes  are  stored  in  an  Integer  array  with  values  of  0 to  255. 
Computer  Processing  Development 

Once  the  short-pulse  radar  data  were  in  a format  compatible  with  the 
CYBER  74,  then  the  data  files  were  operated  upon  by  a program  called  "FREQ" 
to  plot  the  corresponding  temporal  and  spectral  responses.  Simply  plotting 
the  data  in  each  data  file  will  yield  the  time  plot  since  the  data  out  of 
the  radar  is  already  in  this  form.  However,  the  spectral  response  is  ob- 
tained by  range  gating  the  temporal  data  to  include  only  the  signature  of 
interest,  removing  a constant  amplitude  bias  introduced  by  the  A/D  conversion 
processing  and  then  performing  a Fast  Fourier  Transform  (FFT)  on  the  gated 
temporal  data. 

As  indicated  in  Figure  26,  the  range  gating  is  performed  manually  rather 
than  automatically  as  with  the  NBS  data.  Since  the  data  was  for  dry  soil 
only,  it  was  easy  for  a human  operator  to  locate  the  signature  of  Interest. 

Of  course,  the  automated  placement  of  the  range  gate  is  preferred;  however. 
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the  range  gating  algorithm  developed  for  the  NBS  data  would  not  work  for  the 
short-pulse  data  without  modification  because  the  transmitted  waveform  was 
• different.  Time  did  not  allow  such  a modification;  thus,  manual  range  gating 

was  used. 

^ In  the  A/D  conversion  process,  the  Input  analog  signal  was  represented 

on  the  A/D  converter  output  by  8-blts — an  output  scale  ranging  from  0 to  255. 
Thus,  negative  Input  voltages  as  well  as  positive  Input  voltages  are  repre- 
sented by  positive  numbers  on  the  A/D  converter  output.  Effectively,  a DC 
bias  of  128  has  been  added  by  this  conversion  process.  This  offset  was 
removed  before  plotting  the  temporal  response  and  before  taking  the  Fourier 
Trans  form. 

The  gated  spectral  responses  so  obtained  are  utilized  by  the  classifi- 
cation programs  to  be  described  later.  It  was  important  to  generate  a spectral 
response  in  a format  compatible  with  the  classification  programs  already 
developed  for  the  NBS  data.  The  NBS  frequency  response  consisted  of  128 
lines  which  were  separated  by  multiples  of  24  MHz.  Thus,  a frequency  span 
of  0 to  3.048  GHz  was  covered  by  the  128  lines.  However,  only  every  forth 
line  was  used  for  classification.  Thus,  tue  frequency  separation  between 
lines  used  was  96  MHz.  To  make  the  short-pulse  data  yield  a comparable  line 
spectra,  only  every  third  temporal  data  point  was  used.  This  provided  a line 
spectra  with  line  separations  of  about  104  MHz.  The  resulting  spectrum 
spanned  frequencies  from  0 to  3.33  GHz.* 

As  explained  earlier,  the  radar  receiver  output  signals  are  expanded  by  a 
ratio  of  10^:1  by  a special  sampling  technique  within  the  short-pulse  radar. 
Thus,  the  output  signals  are  actually  at  millisecond  and  KHz  rates.  However, 
for  comparison  purposes,  the  unexpanded  frequencies  and  times  are  used  in 
this  discussion  (nanoseconds  and  gigahertz). 
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SECTION  IV 

MATHEMATICAL  MODEL  OF  BURIED  TARGET  AND  ENVIRONMENT 

When  the  measured  data  was  processed  and  examined , it  was  noted  that 
radar  returns  from  the  same  target  at  the  same  depth  varied  with  soli  mois- 
ture. Since  consistency  In  target  return  signatures  seemed  lacking  In  the 
measured  data,  a generalized  mathematical  model  was  developed  for  the  metal 
plate  and  the  type  A,  B,  and  C targets  burled  In  soil.  This  model  provided 
an  understanding  of  the  target  and  the  environment.  The  mathematical  model 
was  programmed  on  a general  purpose  computer  to  yield  plots  of  the  temporal 
and  spectral  characteristics  of  various  type  target  returns.  Parameters 
such  as  target  depth  and  soil  moisture  were  varied,  and  changes  In  target 
return  characteristics  were  observed  and  compared.  A discussion  of  the 
mathematical  model  and  Its  describing  equations  follows. 

DEVELOPMENT  OF  MATHEMATICAL  MODEL 

The  target  and  Its  surrounding  environment  Is  modeled  as  shown  In 
Figure  38.  The  source  and  receiver  of  the  electromagnetic  (EM)  waves  are 
located  above  the  ground,  and  the  waves  are  directed  Into  the  ground  by  an 
appropriate  antenna.  The  EM  waves  are  reflected  by  the  ground  surface  and 
dielectric  discontinuities  in  the  soil.  A target  may  consist  of  several 
different  layers  of  material,  each  of  a different  dielectric  constant,  as 
illustrated  in  Figure  38.  The  model  shown  is  deliberately  simplistic  so 
the  describing  mathematics  are  tractable.  The  following  simplifying  assump 
tlons  are  made  In  the  model  development: 
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Figure  38.  Target  and  Environment  Model, 
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(1)  The  incident  waveform  consists  of  many  monochromatic  plane  waves. 

(2)  The  ground  surface  is  flat  and  perpendicular  to  the  incident 
wavefront. 

(3)  The  target  is  totally  represented  by  layers  of  material  each 
with  a different  dielectric  constant. 

(4)  The  boundary  between  each  layer  in  infinite  in  extent  and  is 
perpendicular  to  the  incident  EM  wavefront. 

The  model  can  be  mathematically  equated  to  an  N-section  transmission 
line.  There  is  a one-to-one  correspondence  between  each  section  of  trans- 
mission line  and  each  dielectric  region  layer  in  Figure  38.  As  each  section 
of  the  transmission  line  has  its  own  characteristic  Impedance,  so  does  each 
region  or  medium  through  which  the  EM  wave  passes  in  the  target  model.  The 
characteristic  Impedance  of  the  ^th  layer  is  given  by 


Z (i)  = 376^ 


ohms ;i*=l,  2,  ...,N 


(13) 


where  is  the  relative  dielectric  constant  of  the  ^th  layer.  The  input 
impedance  of  the  l^th  layer  (at  the  upper  surface  of  the  layer  and  looking 
down  into  the  layer)  is 


Zin(l)  = Z^(i) 


Z (i)  cosh 

Li 


(Yi^i)  + (i)  sinh 


Z^(i)  cosh 


(Yi^i)  + Zj^(i)  sinh 


iliVl 

(Yi^i)J 


(14) 


where 


Z (1)  is  the  impedance  seen  at  the  bottom  of  the  d^th  layer  (looking  down) , 

Li 
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Is  the  thickness  of  the  ^th  layer,  and 
is  the  propagation  constant  of  the  ^th  layer. 


The  propagation  constant  is  related  to  free  space  wavelength  (X  ) and  rela- 

o 

tive  dielectric  constant  (e)  as  follows; 


2tt  . 


Yi  = j ; i = 1,  2 N 


In  all  of  the  equations  describing  the  model,  the  first  layer  is  assumed  to 

be  the  bottom  layer.  The  second  layer  is  the  next  highest  region,  etc. 

In  Equation  (14)  given  above,  the  load  impedance  of  the  ^th  layer  is 

equal  to  the  input  impedance  of  the  (i  - 1)^^  layer.  Since  the  first  or 

bottom-most  layer  is  assumed  to  be  infinite  in  depth  or  thickness,  its  input 

impedance  is  equal  to  its  characteristic  impedance,  (1) . Thus,  the  load 

Impedance  of  the  second  layer  is  the  input  impedance  of  the  first  layer,  or 

Zq(1).  Also,  the  input  impedance  of  the  second  layer  is  a function  of  its 

load  Impedance,  Z (1),  its  own  characteristic  impedance,  Zo(2),  and  its  own 
o 

propagation  constant  and  length,  Yo  The  input  impedance,  Z.  (2), 

L 2 in 

becomes  the  load  impedance  for  the  third  layer,  Z (3) , and  so  forth.  Tlius, 

Li 

the  input  impedance  to  the  or  top-most  layer  can  be  determined  by  use  of 
Equation  14.  This  impedance  characterizes  the  response  of  the  incident  elec- 
tromagnetic wave  to  the  target  and  its  environment.  The  reflected  EM  spectrum 
is  related  to  the  transmitted  spectrum  and  the  reflection  coefficient  as 
follows ; 


Sj^(f)  = p (f)  • S^(f) 
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where  the  reflection  coefficient  is  related  to  the  layer  input  impedance 
in  the  following  manner: 


P (f)  “ 


-1„  <■*>  - "s 
"in  '»  + "s 


The  term,  Zg,  is  the  source  impedance  and  is  assumed  to  be  equal  to  the 
characteristic  impedance  of  the  layer. 

Note  that  all  of  the  equations  derived  thus  far  are  for  a monochromatic 
or  single  frequency  wavefront.  For  each  frequency  or  spectral  line  of  the 
transmitted  spectrum,  a particular  input  Impedance  and  complex  reflection 
coefficient  is  obtained. 

If  the  medium  through  which  the  EM  wavefront  passes  is  lossy,  the 
dielectric  constant  will  be  complex  valued  as  given  below: 


e = - 3 


L "rJ 


where  e and  e denote  the  real  and  imaginary  components,  respectively.  The 
ratio  of  the  imaginary  to  real  component  is  a measure  of  the  attenuation  of 
EM  waves  in  the  medium  and  is  called  the  "loss  tangent"  of  the  dielectric 
material.  The  value  of  both  components  of  a dielectric  constant  is  a function 
of  both  frequency  and  temperature.  After  a brief  search,  very  little  informa- 
tion on  the  dielectric  constant  of  dry  and  moist  soil  as  a function  of  fre- 
quency and  temperature  was  obtained.  The  information  given  in  Table  5 was 
obtained  from  the  sources  indicated. 
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Table 

5 

DIELECTRIC 

CONSTANT  OF 

VARIOUS  TYPES  OF 

SOIL 

Soil  Type 

Degree  of 
Moisture 

Frequency 

Temperature 

Dielectric 

Constant 

Sand  [3] 

Dry 

1 MHz 

20“C 

2.5 

Sand  [3] 

15%  Water 

1 MHz 

ZO^C 

9 

"Soil"  [3] 

Dry 

1 MHz 

20°C 

3 

"Soil"  [3] 

Moist 

1 MHz 

20°C 

10 

Earth  [4] 

Dry 

— 

— 

5 

Earth  [4] 

Wet 

10 

Table  5 leads  one  to  the  conclusion  that  the  dielectric  constant  of  dry  soil 
at  1 MHz  and  ZO^C  varies  from  2.5  to  5 depending  on  the  soil  type.  Also, 
the  dielectric  constant  for  wet  soil  at  the  frequency  and  temperature  indi- 
cated is  somewhat  independent  of  soil  type  and  has  a value  of  9 or  10.  Un- 
fortunately, the  cited  references  did  not  give  any  information  on  the  loss 
tangent  for  soil. 

There  is,  however,  a lot  of  information  on  the  dielectric  constant 
and  loss  tangent  versus  frequency  and  temperature  for  pure  water.  [5]  A 
summary  of  this  information  is  given  in  Figure  39  and  by  the  statement, 

"For  frequencies  between  300  and  3,000  me,  the  power  factor  (loss  tangent) 
at  any  temperature  is  linearly  proportional  to  the  frequency,  while  the 
dielectric  constant  is  slightly  higher  at  the  lower  frequencies."  [6] 
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Figure  39.  The  dielectric  constant  and  loss  tangent  of  water  at  a 
frequency  of  3,000  me. 


To  exercise  the  mathematical  model,  assumptions  concerning  the  dielec- 
tric constant  and  loss  tangent  of  soil  of  various  moisture  contents  were  made 
based  on  the  meager  information  given  above.  The  soil  was  modeled  as  follows 

(1)  The  real  part  of  the  dielectric  constant  is  independent  of 
frequency  for  all  soil  moisture  levels. 

(2)  The  real  portion  of  the  dielectric  constant  varies  linearly 
from  3 to  10  with  moisture  content.  Thus,  a dry  soil  would 
have  a dielectric  content  of  3 while  a very  wet  soil  (30%) 
would  have  a dielectric  constant  of  10. 

(3)  The  loss  tangent  varies  linearly  with  frequency  and  moisture 
content.  For  a very  dry  soil,  the  loss  tangent  is  assumed 
to  be  zero.  For  a very  wet  soil  (30%),  the  loss  tangent  is 
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I 

f 


# 

taken  as  that  of  pure  water;  i.e., 
tan  6=0. 06  X 10 
where  t is  in  hertz. 

(4)  The  soil  was  assumed  to  be  at  a temperature  of  20°C. 

These  characteristics  for  the  soil  may  not  be  entirely  accurate  but  will  serve 

to  provide  some  insight  into  the  effect  of  soil  moisture  on  target  returns. 

Type  A and  C mines  are  assumed  to  be  filled  with  paraffin  wax  with  a 

-4 

dielectric  constant  of  2.2  and  a loss  tangent  of  2 x 10  , both  of  which 

are  approximately  independent  of  frequency.  [5]  Type  A and  C targets  are 
assumed  to  be  5 inches  and  3 inches  thick,  respectively. 

The  model  shown  in  Figure  38  will  yield  a ground  reflection  as  well  as 

* 

a target  reflection.  In  some  instances  it  is  desirable  to  remove  or  eliminate 
the  ground  reflection  so  the  characteristics  of  the  target  can  be  more  clearly 
examined.  This  can  be  accomplished  with  the  model  by  assuming  the  source  and 
receiver  of  electromagentic  waves  are  immersed  in  a medium  of  the  same  dielec- 
tric constant  and  loss  tangent  as  the  subsurface  soil.  Of  course,  the  appro- 
priate electric  length,  !i\,  will  have  to  be  used  to  simulate  the  actual  elec- 

o 

trical  length,  In  addition,  the  source  impedance  is  assumed  to  be  matched 

to  the  environment  impedance. 

EFFECT  OF  WET  SOIL  ON  TARGET  RETURNS 

The  effect  of  wet  soil  on  target  returns  can  be  determined  qualitatively 
by  considering  a simple  target  and  its  environment;  e.g.,  a metal  plate  in  wet 
soil.  Also  assume  that  the  electromagnetic  source  and  receiver  are  also 
Immersed  in  a medium  of  the  same  complex  dielectric  constant  as  the  soil 
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surrounding  the  metal  plate.  In  this  sltuacion,  the  monochromatic  plane  wave 
impinging  on  the  metal  plnte  1'^  described  by 


(19) 


where  x is  the  distance  the  wave  travels, 

Y is  the  complex  propagation  constant, 

t is  the  time, 

w is  the  radian  frequency,  and 

is  the  magnitude  of  the  electric  field. 

Assuming  the  metal  plate  has  zero  impedance,  the  reflected  wave  is  given  by 


E 


E e^ 


(20) 


The  reflection  coefficient,  or  tne  position  of  the  incident  wave  that  is 
reflected,  is  defined  as 


P = ^ (21) 

E 

Substituting  Equations  (19)  and  (20)  into  the  above  equation  gives  the  follow- 
ing expression  for  the  monochromatic  reflection  coefficient; 


P 


-2yx 
= -e  ' 


(22) 


where  the  propagation  constant  is  defined  by  Equation  (15)  as  was  discussed 
earlier.  The  reflected  voltage  spectrum  is  equal  to  the  product  of  the  trans- 
mitted spectrum  and  the  reflection  coefficient 
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Sj^  (f)  = P (f)  * (f) 


-2Y(f)  • S„  (f) 
-e  1 


(23) 


Since  the  propagation  constant  Is  complex  valued.  It  can  be  represented  by 


Y (f)  = a(f)  + jB  (f) 


and  Equation  (23)  becomes 


(24) 


Sr  (f)  = -e 


-2a(f)x 


e-32B(f)x  _ 


(25) 


Note  that  the  first  term  In  the  above  equation  represents  the  attenuation 
of  the  reflected  spectrum  whereas  the  second  term  only  represents  a phase 
shift  that  varies  with  distance  and  frequency.  The  attenuation  factor  varies 
with  frequency  and  is  given  by 


2nf 


« (O  - \ {J  e 


■I 


'soil J 

Using  the  model  for  wet  soil  described  earlier, 

e = Et,  [1  -J  *06  X 10"^  Kf] 
soil  R 


(26) 


(27) 


where  K equals  zero  for  dry  soil  and  unity  for  very  wet  soil,  the  attenuation 
factor  becomes 


a (f)  - Rg  I j ^ ^ *0^  X 10  ^ Kfj 

a (f ) — v/eT  ( .03  x lO"^  Kf ) ; f < 3 x 10^  Hz 

c R 


(28) 


UNCUSSIFIED 

90 


• \ 


UNCLASSIFIED 

Note  that  a(f)  varies  as  the  square  of  the  frequency  In  this  equation.  The 
attenuation  versus  distance  into  the  soil  is  shown  in  Figure  40  for  a medium 
wet  soil  (K  “ y)  as  a function  of  frequency.  The  metal  plate  is  assumed  to 
be  buried  6 inches  be’ow  the  soil  surface  in  this  example.  The  effect  of 
soil  moisture  is  to  attenuate  the  hi^  frequencies  more  than  the  lower  fre- 
quencies. In  this  example,  if  the  soil  were  assumed  to  be  very  wet  (K  = 1, 

e = 10),  the  3 GHz  attenuation  would  be  about  50  dB.  Thus,  the  reflected 
K 

spectrum  is  not  only  an  attenuated  version  of  the  transmitted  spectrum  but 
is  also  distorted. 

APPLICATION  TO  NBS  MEASUREMENTS 


The  computer  program  which  calculates  the  reflected  spectrum  from  the 
transmitted  spectrum  and  other  parameters  was  exercised  for  three  specific 
targets — the  metal  plate,  type  A mine,  and  type  B mine — assuming  a transmitted 
spectrum  corresponding  to  the  processed  NBS  data.  The  results  are  presented 
in  Figures  41  through  49  and  discussed  in  the  following  paragraphs.  In  each 
figure  the  magnitude  of  the  reflected  spectrum  is  given  assuming  the  trans- 
mitted spectrum  is  amplitude  weighted  to  provide  40  dB  Tchebychev  temporal 
sldelobes.  The  frequency  scale  runs  from  zero  to  3.048  GHz.  Also  included 
in  each  figure  is  the  reflected  temporal  waveform  obtained  by  taking  the 
Fourier  transform  of  the  complex  valued  spectrum.  In  each  case,  the  time 
scale  runs  from  zero  to  20-5/6  nanoseconds. 

For  each  type  of  target,  three  soil  moisture  values  were  assumed, 
yielding  three  figures  for  each  target.  The  soil  loss  factor,  K,  and  the 
dielectric  constant  (real  portion)  used  are  indicated  in  each  figure. 

»i 
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In  each  case,  the  top  of  the  target  is  assumed  to  be  6 inches  below  the  ground 
surface.  The  number  in  the  upper  right  comer  of  each  figure  is  the  maximum  "* 

magnitude  of  the  waveform  that  is  shown.  It  is  used  primarily  to  determine  » 

attenuation. 

Figures  41  through  43  are  spectrum  and  temporal  plots  for  a metal 
plate  return  with  loss  factors,  K,  of  0,  .5,  and  1,  respectively.  In  Figure  41 
where  there  is  no  loss  assumed  and  the  magnitude  of  the  reflection  constant 
is  unity,  the  spectrum  shown  is  also  the  spectrum  of  the  "transmitted"  spec- 
trum. The  spectrum  has  a Tchebychev  weighting  to  provide  40  dB  temporal 
sidelobes.  The  temporal  radar  return  is  shown  in  Figure  41(b)  and,  except 
for  an  inversion,  is  also  the  "transmitted"  time  waveform.  Figure  42  indicates  » 

the  effect  of  a moderate  loss  factor,  0.5.  Note  that  the  spectrum  has  been 
attenuated  at  the  high  frequencies  more  than  at  the  lower  frequencies.  Thus, 

♦ 

the  spectrum  is  skewed  towards  the  lower  frequencies.  The  time  waveform  in 
42(b)  also  reflects  this  effect  by  having  a larger  time  spread.  Also,  the 
period  of  the  time  response  oscillation  is  larger.  The  position  of  the 
reflected  waveform  relative  to  the  origin  has  been  changed  due  to  the  larger 
dielectric  constant  used;  e.g.,  6.5  Instead  of  3.  Figure  43  consists  of 
spectral  and  temporal  plots  for  a large  loss  factor  of  unity.  The  character- 
istics outlined  for  Figure  42  are  now  more  pronounced.  The  magnitude  of  the 
temporal  responses  in  Figures  42  and  43  are  attenuated  relative  to  the  loss- 
less  case  by  5.2  and  10.4  dB,  respectively 

The  spectral  and  temporal  plots  for  the  type  A mine  model  are  shown 
in  Figures  44  through  46  for  loss  factors  of  0,  .5,  and  1,  respectively. 
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Attenuation  Versus  Depth  for  a Metal  Plate 
Target  in  Moderately  Wet  Soil  (K  = 0.5, 

= 6.5). 


Figure  40 
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Figure  41.  Math  model  of  radar  return  from  metal  plate 
burled  in  a lossless  medium. 


Math  model  of  radar  return  from  a metal 
plate  burled  6"  in  a lossy  medium  (K“0.5) 
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Figure  43.  Math  model  of  radar  return  from  metal 

plate  buried  6"  in  lossy  soil  (loss  equal 
to  that  of  pure  water  at  20®  C,  K=1.0). 


UNCLASSIFIED 

96 


9 


AMPLITUDE  AMPLl 


UNCLASSIFIED 


Frequency  (GHz) 
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Math  model  of  radar  return  from  non-metalllc 
mine  type  A burled  6"  In  a lossy  medium  (loss 
equal  to  that  of  pure  water  at  20'’C,  K-1.0) 


UNCLASSIFIED 


The  type  A mine  model  consists  of  two  dielectric  discontinuities  and,  therefore, 
two  reflected  pulses  are  shown  in  the  time  domain  plots.  The  distance  between 
pulses  is  determined  by  the  thickness  of  the  mine  and  the  dielectric  material 
inside  the  mine.  Note  that  for  the  type  A mine,  the  reflected  pulses  are 
close  enough  together  to  interfere  with  one  another.  Thus,  a particular  com- 
posite pulse  is  formed.  The  corresponding  spectrum  is  also  shown.  Note 
that  several  nulls  are  present  in  the  spectrum  because  the  two  reflections 
are  equal  in  magnitude  and  opposite  in  sense  (180°  out  of  phase)  at  certain 
frequencies.  As  the  attenuation  factor  and  dielectric  constant  are  increased, 
the  spectrum  has  a greater  skew  toward  the  low  frequencies.  The  spectrum  is 
thereby  changed  from  that  of  the  lossless  case.  Also,  the  composite  temporal 
waveform  changes  as  shown  in  Figures  45(b)  and  46(b). 

Figures  47  through  49  contain  plots  of  the  spectral  and  temporal  returns 
from  the  type  C target  model.  Loss  factors  of  0,  .5,  and  1 are  also  assumed 
for  these  figures.  Since  the  type  C mine  is  modeled  with  a smaller  thickness 
of  3 inches,  the  spectrum  nulls  occur  at  frequencies  different  from  those  for 
the  type  A mine.  Also,  the  time  responses  interfere  in  a different  manner. 
Increased  attenuation  by  the  moisture  in  the  soil  changes  the  spectrum  and 

time  responses  as  indicated  in  Figures  48  and  49. 

APPLICATION  TO  SHORT-PULSE  RADAR  MEASUREMENTS 

The  short-pulse  measurements  were  made  late  in  this  program  and  were 
not  modeled  by  the  mathematical  model  previously  described.  The  major  differ- 
ence in  applying  the  model  to  the  NBS  data  and  the  short-pulse  radar  data  is 
in  the  spectrum  of  the  transmitted  pulse.  The  spectrum  of  the  NBS  data  spans 
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Figure  47.  Math  model  of  radar  return  from  non- 

metalllc  mine  type  C buried  6"  in  loss- 
less medium. 
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Figure  48.  Math  model  of  radar  return  from  non- 
metalllc  mine  type  C buried  6"  in  a 
lossy  medium  (K=0.5). 
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from  120  MHz  to  3.048  GHz,  whereas  the  spectrum  of  the  short-pulse  radar  j 

data  spans  from  approximately  100  MHz  to  2 GHz.  Thus,  the  temporal  waveform  ! 

of  the  short-pulse  radar  is  longer  in  duration  than  the  pulse  of  the  NBS  data.  * > 

In  addition,  the  spectral  magnitude  and  phase  responses  are  different  since 
the  two  temporal  waveforms  differ  in  detail. 

The  theoretical  spectral  and  temporal  responses  of  the  short-pulse 
radar  are  shown  in  Figures  22(a)  and  22(b).  They  could  easily  be  modeled 
by  the  computer  program  just  described;  however,  because  of  lack  of  time 
they  were  not. 

COMPARISON  OF  PROCESSED  NBS  DATA  AND  MATHEMATICAL  MODEL 

As  explained  earlier  in  this  report,  the  NBS  raw  measured  data  was 
weighted  and  transformed  to  the  time  domain  for  range  gating.  The  range  gate 
was  positioned  so  that  the  ground  return  and  other  interferring  returns  were 
partially  or  totally  eliminated.  This  gated  return  was  then  transformed  back 
into  the  frequency  domain.  The  gated  spectrum  and  the  temporal  response 
were  plotted  in  various  formats  for  examination. 

In  this  subsection,  several  of  the  processed  NBS  data  plots  are 
discussed  and  compared  with  the  responses  obtained  with  the  mathematical 
model  presented  in  the  proceeding  subsection.  The  mathematical  model  resoonses 
presented  earlier  were  for  some  representative  parameter  values  and,  therefore, 
did  not  correspond  exactly  to  the  processed  data  with  which  it  was  compared. 

However,  general  characteristics  of  both  the  processed  data  and  the  mathe- 
matical model  responses  were  similar. 
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The  processed  data  and  model  for  the  metal  plate  buried  6 inches  below 
the  ground  surface  were  compared.  Figure  50  is  a plot  of  the  processed 
measured  data.  Figure  50(a)  shows  the  spectrum  of  the  measured  data  without 
weighting,  and  Figure  50(b)  shows  the  corresponding  time  responses  with 
weighting  and  with  the  time  gate  boundaries  marked  with  two  vertical  lines. 

The  temporal  response  between  these  lines  was  taken  as  the  metal  plate  return. 
Figure  50(c)  is  the  spectrum  of  the  gated  portion  of  the  time  response,  and 
Figure  50(d)  is  the  gated  spectrum  phase  response  ranging  between  the  limits 
of  -TT  to  +ir  radians . This  data  is  from  the  measured  data  record  GT5-68  which 
is  a return  from  a metal  plate  buried  6"  deep  in  loam  with  a 7%  soil  moisture 
content . 

Compare  the  gated  spectrum  of  this  data  with  the  spectrum  obtained  with 
the  mathematical  model  as  shown  in  Figure  41.  There  is  a great  deal  of  simi- 
larity between  the  two.  One  notes  that  the  spectrum  in  Figure  50  appears  to 
be  skewed  slightly  toward  the  lower  frequencies.  This  is  because  the  7% 
moist  loam  does  provide  some  attenuation  to  the  incident  EM  waves.  However, 
the  amount  of  attenuation  presented  is  not  as  great  as  that  shown  in  Figure  42 
where  the  loss  factor  is  0.5.  The  spectral  phase  appears  to  be  linear  over 
most  of  the  frequency  range.  Based  on  the  mathematical  model,  one  would 
expect  this  result.  The  gated  temporal  response  of  Figure  50  is  somewhat  simi- 
lar to  that  of  Figure  41,  but  differs  in  detail.  Comparison  of  the  type  A 
and  type  B target  data  with  the  math  model  plots  found  in  Figures  44  through 
49  yield  the  same  result;  i.e.,  the  spectral  plots  were  in  very  close  agree- 
ment, but  the  temporal  signatures  differed. 

At  first,  it  was  thought  the  math  model  was  in  error;  however,  the 
following  simplified  analysis  shows  that  the  results  obtained  with  the  math 
model  were  Indeed  correct  and  the  error  was  in  the  data  measurements. 
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Figure  50.  Measured  data — metal  plate  6"  below  ground,  7%  moisture— 
GT5-68. 
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Consider  for  a moment  the  behavior  of  electromagnetic  waves  at  the 
boundary  between  two  media  as  Illustrated  In  Figure  51.  The  reflection 
coefficient  at  the  boundary  Is  given  by 

^2  - ^1 

^ C2‘ 

where  and  N2  are  the  wave  Impedances  of  media  1 and  2,  respectively.  For 
a perfect  conductor,  the  wave  Impedance  Is  zero;  for  a non-conducting,  non- 
ferrous  medium,  the  wave  Impedance  Is 


where  Is  the  magnetic  permittivity  of  space,  and  e Is  the  dielectric  con- 
stant of  the  medium.  If  the  wave  Impedance  for  free  space  (air)  Is  denoted  by 


then  equations  (29 ) and  (30 ) above  may  be  written  as 


and  p = 


vr 

r 

VSi 


where  e Is  the  relative  dielectric  constant  of  the  medium, 
r 

Note  that  If  medium  1 has  a smaller  relative  dielectric  constant  than 
medium  2,  then  p has  a negative  value.  On  the  other  hand.  If  medium  1 has  a 
larger  relative  dielectric  constant  than  medium  2,  then  p Is  positive.  Also, 


note  that  for  metal,  p Is  also  negative;  l.e.,  p = -1. 
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Suppose  a metal  plate  Is  burled  In  soil  as  Illustrated  In  Figure  32. 

♦ 

Also,  assume  the  transmitted  pulse  Is  of  the  shape  shown  In  the  same  figure. 

What  should  the  radar  return  from  the  alr/soll  and  the  soll/metal  Interfaces 

look,  like?  The  soil  has  a dielectric  constant  greater  than  that  of  air; 

thus,  p at  the  alr/soll  Interface  has  a negative  value.  Likewise,  p at  the  * 

soll/metal  Interface  Is  negative.  Thus,  the  radar  return  should  appear  as 

shown  In  Figure  52.  Assuming  the  metal  plate  Is  deep  enough  below  the  soil 

surface  so  the  two  reflected  pulses  do  not  overlap,  each  of  the  returns 

should  appear  as  the  Inverse  of  the  transmitted  waveform  as  Illustrated. 

Now  refer  to  the  measured  data  for  a metal  plate  burled  In  soil  shown  In 

* 

Figure  50.  Note  that  the  plate  return  Is  not  a replica  of  the  ground  return 
as  theory  Indicates  It  should  be. 

Extending  this  simplified  analysis  to  a mine  burled  In  the  soil  yields  * 

the  expected  returns  illustrated  in  Figure  53.  A comparison  of  the  measured 
mine  temporal  signatures  with  those  of  this  figure  also  Indicates  a differ- 
ence and  that  something  is  wrong. 

After  an  extensive  search  for  the  reason  for  this  discrepancy,  it  was 
found  that  If  the  measured  spectra  lines  were  shifted  down  in  frequency  by 
about  24  MHz,  the  measured  data  and  the  analytically  derived  signatures  agree. 

For  example,  compare  the  measured  temporal  signature  of  the  metal  plate  re- 
turn after  the  spectrum  has  been  down-shifted  about  24  MHz  as  shown  in  ‘ 

Figure  54  to  the  returns  derived  in  Figure  52.  Except  for  amplitude,  there 
is  ^^mplete  agreement!  Figures  55  through  58  show  the  measured  temporal 
responses  of  the  type  A and  type  C mines  before  and  after  the  24  MHz  spectral 
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Figure  51.  Behavior  of  EM  Waves  at  a Boundary. 
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Figure  52.  Metal  Plate  Buried  in  Soil. 
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Figure  53.  Non-metalic  Mine  Buried  in  Soil. 
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Figure  55.  Type  A mine,  6"  depth,  7%  moisture, 
with  time  distortion. 
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shift.  A comparison  of  these  responses  of  the  spectral  shifted  data  with 
those  derived  in  the  simplified  analysis  and  also  those  derived  via  the 
mathematical  model  indicates  almost  complete  agreement.  Apparently,  when 
the  NBS  data  was  collected,  an  error  was  made  in  the  frequency  assignment. 
After  discovering  this  error,  all  of  the  NBS  data  was  corrected  by  a spec- 
tral shift. 
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SECTION  V 

SELECTION  AND  DEVELOPMENT  OF 
DISCRIMINATION / CLASS I FI CATION  ALGORITHMS 

The  process  of  initially  detecting  a potential  target  and  then  deter- 
mining not  only  if  it  is  a mine,  but  its  type,  is  divided  into  three  steps 
as  illustrated  in  Figure  59.  The  first  step  is  to  separate  potential  target 
returns  from  radar  returns  due  to  soil  anomalies.  For  this  study,  a poten- 
tial target  is  considered  to  be  any  buried  object  which  gives  a return  sig- 
nificantly larger  than  from  the  soil  alone.  Some  examples  of  buried  non- 
mine objects  one  might  encounter  in  a typical  field  are  rocks,  roots,  metal 
scraps,  etc.  The  discrimination/classification  steps  which  follow  the 
detection  process  involve  a relatively  complicated  and  time-consuming  process. 
Thus,  not  every  location  in  a field  can  be  subjected  to  all  three  steps  shown 
in  Figure  59.  The  purpose  of  the  detection  processor  is  to  drastically 
reduce  the  number  of  locations  in  a field  that  must  be  scrutinized  by  all 
three  sections  of  the  processor. 

The  second  processor  section  is  referred  to  as  the  discriminator.  The 
purpose  of  this  step  is  to  determine  which  of  the  potential  targets  are  indeed 
targets;  i.e.,  mines.  Later  in  this  report,  the  algorithm  used  to  accomplish 
this  is  described.  The  third  and  final  section  of  the  three-step  processor 
is  designed  to  indicate  to  which  of  several  classes  of  mines  the  detected 
target  belongs. 

The  detection  process  is  a relatively  simple  one.  The  returns  from  the 
upper  one  foot  of  the  soil  are  subjected  to  a threshold  which  is  normalized 
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to  the  RMS  average  return  from  a portion  of  the  soil  below  one  foot  from  the 
soil  surface.  In  the  region  with  depth  greater  than  one  foot,  mines  are  not 
likely  to  be  present.  This  normalization  of  the  radar  returns  should  remove 
the  effects  of  transmitter  power  variations  and  other  unexplained  variables. 
Thus,  the  false  alarm  rate  should  be  fairly  constant  for  each  type  of  soil. 

If  the  radar  returns  from  a particular  location  do  not  pass  the  threshold 
of  the  detection  processor,  they  are  ignored  and  another  location  is  interro- 
gated. If  the  detection  processor  threshold  is  exceeded,  then  the  algorithm 
of  the  discrimination  processor  is  exercised.  Preferably,  the  false  alarm 
rate  for  the  detection  processor  should  be  about  one  false  alarm  rate  per 
1000  meters  of  travel.  As  will  be  shown  later,  this  is  equivalent  to  a 
probability  of  false  alarm  of  about  10 

If  the  radar  returns  from  a particular  location  do  pass  the  detection 
processor  threshold,  the  radar  return  range  profile  is  searched  for  the  range 
ceil  with  the  largest  amplitude.  A target  signature  cell  is  then  defined 
with  the  largest  amplitude  return  as  its  centroid.  The  spectrum  of  the  time 
signature  within  this  cell  is  then  obtained  by  use  of  the  Fourier  transform. 
The  discrimination  algorithm  described  in  the  next  section  is  exercised  uti- 
lizing the  spectral  signature  data.  Based  on  the  output  of  the  algorithm, 
the  signature  is  declared  to  be  that  of  a mine  or  a non-mine,  e.g.,  a rock 
or  a root.  If  it  is  decided  that  the  signature  is  that  of  a non-mine,  the 
processing  is  terminated  and  the  radar  is  moved  to  the  next  location. 

On  the  other  hand,  if  a mine  is  declared,  the  spectral  signature  is 
scrutinized  to  determine  to  which  class  of  mine  it  belongs.  At  this  point, 
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the  processing  Is  complete  and  the  radar  moves  to  a new  location.* 

The  overall  performance  of  the  three-step  processor  and  radar  depends 
upon  the  operational  scenario.  For  the  sake  of  illustration,  the  following 
two  scenarios  will  be  considered: 

Scenario  A 


- The  width  of  ground  to  be  covered  by  the  mine  detection  radar: 

4 meters  ( '■■'  13  feet)  ; 

- The  rate  at  which  ground  is  to  be  covered  by  the  mine  detection 
radar:  10  KM/hr  (^9  ft/sec); 

- The  false  alarm  rate:  less  than  1 per  1000  meters  of  linear 
travel  (1  per/3280  feet); 

- The  ground  is  sampled  by  the  radar  on  a 1 foot  square  grid — i.e., 
each  radar  measurement  is  located  one  foot  from  all  other  measure- 
ments ; 

- The  ground  is  void  of  any  rocks,  roots,  etc.  Only  mines  or  soil 
exist  within  the  path  of  the  mine  detection  radar. 

Scenario  B 


- The  first  four  statements  under  Scenario  A apply; 

- One-half  of  the  locations  where  radar  measurements  are  taken 
contain  non-mines — i.e.,  rocks,  roots,  etc.  The  other  half  con- 
tain soil  only.  In  addition,  a mine  may  be  present  at  any  of 
the  locations. 


The  radar  is  assumed  to  be  continuously  moving.  The  processing  required  must 
be  completed  before  the  radar  has  moved  to  the  next  position  to  be  interro- 
gated. 
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The  number  of  locations  that  are  Interrogated  by  the  radar  per  second 

• can  be  determined  from  the  scenarios  as  follows : 

« 13  locatlons/llnear  foot  of  travel 

9 linear  feet  of  travel/second 
or  13  X 9 = 117  locations /second. 

* 

This  indicates  how  much  processing  has  to  be  accomplished  In  a given  time 
period  and  will  be  used  in  Section  VII  to  determine  the  processor  hardware 
requirements. 

To  convert  the  specified  false  alarm  rate  into  false  alarm  probability, 
the  number  of  locations  interrogated  by  the  radar  per  1000  meters  must  be 
determined.  As  indicated  above,  3280  linear  feet  or  13  x 3280  = 42,640 
locations  must  be  interrogated  per  one  false  alarm.  Thus,  the  overall  prob- 
ability of  false  alarm  must  be 

Pp^<_  1/42,640^^2.4  X lO"^  (34) 

The  that  each  section  of  the  three-step  processor  must  have  to  meet  this 

requirement  depends  on  which  of  the  above  scenarios  is  assumed.  In  Scenario  A 

the  P of  — '10  ^ may  be  divided  between  steps  1 and  2 since  the  overall  P 

FA  f A 

is  given  by  the  product  of  the  individual  P^^'s.  For  example,  the  P^^ 

-4 

assigned  to  the  detection  processor  could  be  10  while  that  assigned  to  the 
discriminator  could  be  10  This  assignment  would  require  the  discriminator 
to  be  exercised  only  about  once  in  10^  interrogations,  or  once  in  about  100 
seconds.  Consequently,  a false  alarm — i.e.,  an  indication  of  the  presence  of 
a mine  when  one  was  not  present — would  occur  only  once  in  about  10^  interro- 
‘ gations  or  once  in  about  1000  seconds.  On  the  other  hand,  the  relatively 
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simple  processor,  the  detection  processor,  would  be  exercised  for  every  loca- 
tion Interrogated.  In  Scenario  B,  the  Pp^asslgned  to  each  processor  Is  quite 
different  from  that  described  for  Scenario  A.  To  Illustrate,  suppose  there 
was  a rock  or  root  at  each  location  which  was  searched.  The  total  or  overall 
Pp^  for  the  three-step  processor  Is  that  of  the  discriminator.  If  the  Pp^ 
assignment  made  for  Scenario  A is  used,  the  total  Pp^  would  be  10  ^ rather 
than  10  ^ as  before. 

In  Scenario  B,  only  one-half  of  the  interrogated  locations  are  assumed 
to  contain  non-mines.  Thus,  the  total  or  overall  probability  of  false  alarm 

is  related  to  the  P . for  each  step  as  follows: 

FA 

^FA^  " ^1  ^FA^  ^FA2  ^2  ^FA2 

(35) 

Wi  + W2  = 1 


where  is  the  fraction  of  locations  without  non-mines  (soil  only) , and 
W2  is  the  fraction  of  locations  with  non-mines  (rock,  roots,  etc.). 
For  the  situation  assumed  In  Scenario  B, 

Wi  = W2  = Y (36) 

-4 

Thus,  with  the  Pp^  assignment  assumed  for  Scenario  A,  e.g. , Pp^  = 10  , 

Pp^  = 10  the  overall  Pp^  is 

Pp  = i X lO"^  X lO"^  + y X 10"^  X 10“^  (37) 

FA  2 2 2 

One  concludes  from  the  foregoing  analysis  that  is  a fair  proportion 

of  the  interrogated  Icoations  are  expected  to  contain  non-mines,  the  P„. 
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assigned  to  the  discriminator  processor  should  be  low.  To  meet  the  of 
10  ^ for  the  overall  processor,  the  for  the  discriminator  should  be 
approximately  10  ^ . Thus,  in  this  situation,  the  detection  processor  does 
little  in  the  way  of  reducing  the  data  the  discriminator  has  to  process. 

On  the  other  hand,  if  most  of  the  locations  only  contain  soil,  i.e. , Sce- 
nario A,  then  the  detection  processor  can  significantly  reduce  the  data  to 
be  handled  by  the  discriminator/classifier. 

The  remainder  of  the  section  is  divided  according  to  the  various  tech- 
niques considered  to  perform  each  of  the  three  processes  Illustrated  in 
Figure  59,  e.g.,  detection,  discrimination,  and  classification.  First,  the 
performance  of  the  CFAR  type  detection  processor  is  evaluated  using  the  NBS 
measured  data.  Next,  the  pattern  recognition  algorithms  considered  for  dis- 
crimination and  classification  are  evaluated  with  the  NBS  data.  In  addi- 
tion, two  types  of  correlation  processing  are  investigated  for  target  dis- 
crimination. 

CFAR  ANALYSIS 

As  described  earlier,  the  detection  processor  subjects  the  sampled 
radar  range  profile  to  a normalized  threshold.  If  any  of  the  range  profile 
exceeds  the  threshold,  a potential  target  is  declared.  The  threshold  is 
normalized  to  the  RMS  average  of  the  time  profile  representing  the  soil  1 foot 
from  the  surface  or  deeper.  However,  in  this  analysis,  it  was  advantageous 
to  fix  the  threshold  and  normalize  the  data.  The  normalization  was  applied 
to  that  portion  of  the  range  profile  of  interest  (i.e.,  that  representing 
the  soil  1 foot  below  the  surface  or  less).  The  normalized  profile  was 
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then  subjected  to  a fixed  threshold.  This  procedure  is  equivalent  to  the 
process  of  normalizing  the  threshold  described  above  and  is  easier  to  display. 

Each  of  the  NBS  data  measurements  were  normalized  in  this  manner  and 

the  peak  value  of  each  profile  was  noted.  In  Figure  60,  a histogram  of  these 

peak  values  is  given  for  each  moisture  category  and  for  either  a soil-only 

(called  background)  or  a mine-only  case.  In  this  figure,  the  profile 

values  which  contained  a mine  signature  are  denoted  by  an  "X"  and  those 

where  only  a soil  signature  existed  are  noted  by  the  symbol  "0."  In  Figure 

60(a)  the  case  of  7%  moisture  with  loam  soil  is  Illustrated.  The  magnitude 

values  are  spread  out  along  the  axis  where  the  mine-only  and  soil-only  values 

are  non-overlapping.  Thus,  a threshold  can  be  placed  at  about  1/5  of  the 

axis  length  and  a probability  of  detection  of  1.0  can  be  obtained  with  an 

accompanying  probability  of  false  alarm  of  0.0.  This  degree  of  performance 

is  shown  in  Figure  61  where  versus  is  plotted  for  this  case.* 

D FA 

In  Figure  60(b),  the  case  of  17-18%  moisture  with  loam  soil  is  shown. 


i 

I 

I 

• ] 


The  normalized  magnitude  values  for  this  case  are  not  spread  out  as  much  as 
for  the  previous  case  but  are  overlapping  as  shown.  Thus,  the  probability 
of  detection  performance  will  not  be  nearly  as  good  as  that  of  the  7%  loam 
case  discussed  previously.  The  corresponding  performance  curves  indicating 
the  probability  of  detection  versus  probability  of  false  alarm  are  shown  in 
Figure  62.  In  a similar  fashion,  a histogram  of  the  profile  peak  values 
for  the  12-20%  moisture  clay  and  the  26-30%  moisture  clay  cases  is  shown  in 


Since  the  number  of  data  points  used  is  somewhat  limited,  one  cannot,  with 
confidence,  state  that  a perfect  detection  performance  will  be  obtained  in 
a practical  application;  however,  very  good  performance  is  highly  likely. 
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Figure  62.  Potential  target  detection  - CFAR. 
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Figures  60(c)  and  60(d).  In  each  of  these  cases,  the  mine-only  and  soil- 
only  values  overlap  as  in  Figure  60(b).  The  corresponding  performance 
curves  are  shovm  in  Figures  63  and  64,  respectively.  The  histogram  shown 
in  Figure  60(e)  Includes  all  of  the  cases  discussed  previously.  Thus, 
the  soil-only  and  the  mine-only  histograms  overlap  and  the  corresponding 
performance  is  rather  poor.  Since  a probability  of  false  alarm  of  approxi- 
mately 10  or  something  of  that  order,  with  the  corresponding  probability 
of  detection  on  the  order  of  0.9  to  1.0  is  necessary  for  this  mine  detec- 
tion application,  it  appears  that  only  the  first  case  where  we  have  a 7% 
moisture  with  loam  soil  yields  the  performance  required.  The  other  cases 
we  have  discussed  and  illustrated  do  not  yield  performances  that  are  accept- 
able for  this  application.  For  example,  the  performance  curves  shown  in 
Figure  62  for  the  17-18%  moisture  loam  soil  case  yield  a probability  of 
detection  of  0.9  with  a corresponding  probability  of  false  alarm  of  about 
0.3.  This  is  clearly  unacceptable. 

It  is  of  interest  to  consider  the  detection  performance  for  the  vari- 
ous moisture  and  soil  types  previously  discussed  where  the  threshold  is  not 
normalized  but  is  subjected  to  a fixed  threshold.  Figure  65  shows  the  his- 
togram for  the  soil-only  and  mine-only  data  and  indicates  the  degree  of 
overlap  that  occurs  in  the  histogram.  The  results  obtained  are  essentially 
the  same  as  for  the  normalized  data.  In  the  case  of  the  7%  moisture  with 
loam  soil,  the  soil-only  and  the  mine-only  histograms  do  not  overlap;  thus, 

it 

a perfect  detection  performance  can  theoretically  be  obtained.  In  the  case 

*Since  the  number  of  data  points  used  is  somewhat  limited,  one  cannot,  with 
confidence,  state  that  a perfect  detection  performance  will  be  obtained. 
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of  the  17-18%  moisture  with  loam  soil,  the  two  histograms  do  overlap.  The 
same  is  true  of  the  wet  clay  cases  as  shown  in  Figures  65(c)  and  65(d).  The 
corresponding  detection  performance  curves  for  the  cases  illustrated  in 
Figure  65  are  given  in  Figures  66  through  69.  When  these  performance  curves 
are  compared  to  those  for  the  CFAR  case,  one  notes  that  there  is  hardly  any 
difference  at  all.  Thus,  we  can  conclude,  based  on  this  limited  set  of  data, 
that  normalizing  the  threshold,  or  normalizing  the  data  which  is  subject  to 
a fixed  threshold,  does  not  improve  the  detection  performance  or  give  any 
apparent  advantage. 

All  of  the  previous  histograms  are  for  the  case  where  the  transmit 
pulse  was  designed  with  a 40  dB  time  sidelobe.  It  was  thought  that  perhaps 
the  time  sidelobes  of  the  large  ground  return  might  be  interfering  with  the 
detection  of  a much  smaller  return  from  a buried  target,  thus  yielding  a 
poor  detection  performance.  To  determine  if  this  were  indeed  true,  histograms 
similar  to  those  shown  in  Figures  60  and  65  for  the  normalized  and  unnorma- 
lized cases  were  generated  where  the  time  sidelobes  of  the  transmitted  pulse 
were  designed  to  be  at  60  dB  below  peak  value  instead  of  40  dB  as  in  the 
previous  examples.  These  histograms  are  shown  in  Figures  70  and  71,  respec- 
tively. One  notes  that  the  degree  of  overlapping  between  the  soil-only  and 
mine-only  histograms  are  about  the  same  as  for  the  40  dB  data;  thus,  one 
concludes  that  with  this  set  of  data,  the  40  dB  time  sidelobes  are  adequate 
such  that  the  detection  performance  is  not  degraded. 

An  explanation  of  why  the  detection  processor  performs  so  poorly  in 
the  wetter  soils  is  needed.  Perhaps,  a reason  for  the  poor  detection 
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Figure  67.  Potential  target  detection  - unnortnalized . 
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Figure  71.  Histograms  of  time  data  subjected  to  a fixed  threshold 
(CFAR) . Sldelobes  of  the  temporal  waveform  is  60  dB. 
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performance  is  that  the  soil  is  composed  of  roughly  parallel  layers  or 
boundaries  which  allow  deposits  of  water.  Since  the  dielectric  constant 
of  the  water  deposits  is  so  much  larger  than  that  of  the  surrounding  soil, 
the  thin  layers  of  water  reflect  a substantial  portion  of  the  incident 
« microwave  energy  and  thus  Interfere  with  the  detection  of  mines  buried  in 

soil.  Another  explanation  is  that  radar  returns  from  close  proximity  ob- 
jects above  the  ground  are  interferring.  For  example,  perhaps  the  radar 
and  radar  vehicle  are  reflecting  the  radar  signals  to  provide  a multiple 
bounce  situation.  The  sourr.  of  this  problem  was  not  determined  in  this 
study. 

PATTERN  RECOGNITION 

The  information  processing  discipline  known  as  pattern  recognition  is 
implemented  whenever  a decision  is  to  be  made  which  involves  large  quanti- 
ties of  data.  The  object  of  pattern  recognition  or  pattern  classification 
is  to  automate  the  decis.ion  process  based  on  the  experience  of  a learning 
set.  A complete  classification  system  encompasses  both  the  design  and  im- 
plementation of  computer  software  and  hardware.  In  this  section,  only  the 
software  system  will  be  described. 

The  general  approach  is  to  divide  Information  processing  for  the 
purpose  of  pattern  recognition  into  three  steps:  (1)  data  collection, 

• (2)  feature  extraction  and  evaluation,  and  (3)  selection  of  the  classifier. 

Figure  72  gives  a block  diagram  of  the  basic  data  processing  flow  for 
designing  various  classifiers.  The  data  on  magnetic  tape  is  in  the  frequency 

* 

domain.  From  this  information,  radar  signatures  are  generated  and  placed 
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Figure  72.  Classification  data  processing  system. 
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into  a mine  or  non-mine  class . The  non-mine  category  consists  of  returns 
from  objects  which  might  typically  be  found  in  a mine  field,  e.g.,  rocks, 

* roots,  and  flat  metal  plates.  The  mine  category  contains  the  signatures 

for  the  type  A,  B,  and  C mines.  The  preprocessing  steps  for  generating  the 
^ radar  signatures  are  summarized  in  the  following  paragraphs. 

The  raw  data  consisted  of  the  parameters  which  were  measured  at  123 
frequencies  over  a range  of  several  GHz.  The  data  was  weighted  to  reduce 
the  time  sidelobes  of  the  ground  return,  then  transformed  into  the  time 
domain.  A range-gated  window  was  set  automatically  to  separate  the  target 
signature  from  the  ground  return,  and  that  portion  of  the  signature  was  then 
transformed  into  the  frequency  domain.  All  classification  work  used  these 
range-gated  frequency  profiles. 

From  each  signature  a characteristic  set  of  numbers,  i.e.,  features, 

were  extracted  according  to  a prescribed  format.  The  ensemble  of  features 

for  one  particular  class  such  as  a rock,  root,  or  a specific  mine  comprised 

a learning  set  for  that  class.  Several  of  these  learning  sets  were  inserted 

into  the  pattern  recognition  programs  (e.g.,  the  Fisher's  Linear  Discriminant 

or  the  k-nearest  neighbor  method)  to  produce  a classification  algorithm.  The 

classifier  was  designed  to  Identify  an  unknown  sample  as  a mine  or  non-mine 

and  to  minimize  the  number  of  mlsclasslfled  samples. 

Data  Selection 
V 

The  majority  of  classification  work  discussed  in  this  report  was  done 
with  the  NBS  data  which  consisted  of  measurements  of  the  Sj^2  parameters  versus 
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frequency  for  two  types  of  land  mines  and  several  non-mine  objects  such  as 
roots  and  rocks.  The  objects  were  burled  under  three  to  six  Inches  of  soil, 
with  soil  moisture  content  varying  between  three  and  twenty  percent.  Bistatic 
measurements  were  made  with  the  antennas  placed  at  five  positions  with  re- 
spect to  the  target  center:  -8",  -A",  0",  4",  and  8". 

The  selection  of  data  for  classification  purposes  is  one  of  the  most 
critical  considerations  in  the  information  processing  chain.  If  the  data 
base  is  too  small  to  include  signatures  taken  for  all  parameters  of  interest 
(e.g.,  aspect  angle,  soil  moisture,  clutter  level,  etc.),  then  the  resulting 
design  may  do  a good  job  in  classifying  under  some  conditions  but  break  down 
when  signatures  resulting  from  the  neglected  conditions  are  examined.  For 
this  work,  most  of  the  available  data  were  used  for  the  classifier,  with  the 
exception  of  the  data  taken  with  the  very  wet  soil  (record  GT-10)  and  the 
measurements  taken  at  -8"  and  8"  from  the  center  of  the  mine.  The  GT-10 
data  was  discarded  because  the  signal-to-interference  ratio  was  very  low 
and  could  not  be  successfully  used  for  classification.  Measurements  at  -8" 
and  8"  were  left  out  because  they  sometimes  yielded  small  signal  levels. 

At  these  antenna  positions,  the  mine  was  only  partially  illuminated.  The 
former  case  limits  the  applicability  of  the  classification  performance  to 
soils  with  moisture  content  less  than  about  20%.  The  latter  case  would  have 
given  over-estimated  false  alarm  rates. 

Feature  Extraction 

Feature  extraction  is  the  process  whereby  a relatively  small  number 
of  parameters  are  computed  from  the  input  measurements.  An  advantage  of 
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this  step  lies  in  the  substantial  savings  of  computer  time,  provided  this 
task  can  be  accomplished  without  loss  of  information  and  without  sacrificing 
performance.  For  this  reason,  the  parameters  should  be  selected  so  that  they 
best  represent  the  data  for  the  purpose  of  making  a decision.  Once  the 
specified  parameters  were  extracted  from  each  radar  signature,  they  were 
treated  as  components  of  a vector  in  N-dimensional  feature  space.  The  size 
of  the  feature  vector  and,  hence,  the  dimensionality  of  this  space  varied 
between  10  and  32  dimensions;  in  other  words,  the  feature  vector  size  varied 
according  to  the  number  of  parameters  selected  per  signature.  Hence,  each 
signature  represents  one  feature  vector,  but  each  vector  can  also  be  thought 
of  as  a point  in  feature  space,  the  location  of  which  is  determined  by  the 
vector's  tip.  The  collection  of  points  in  feature  space  for  a particular 
mine  or  non-mine  is  known  as  a learning  set.  Several  such  learning  sets 
constituted  the  input  to  the  classifier  programs.  Actually,  two  different 
types  of  learning  sets  were  formed  from  the  gated  data  for  the  mine, 
root,  rocks,  and  metal  plate.  For  the  first  set,  no  intuition  was  used  in 
selecting  the  features.  It  contained  the  parameters  of  each  signature 
sampled  at  constant  frequency  increments.  The  real  and  imaginary  components 
were  used  to  obtain  the  absolute  value  so  as  to  reduce  the  vector  dimen- 
sionality. (Another  advantage  of  this  step  will  be  discussed  in  a later 
paragraph  on  feature  size . ) An  example  of  a vector  formed  from  the  sampled 
frequency  representation  of  the  S^2  data  profile  is  shown  in  Figure  73. 

The  other  set  contained  ten  intuitive  features  for  each  signature  which 
were  selectively  chosen  from  an  original  larger  set  to  maximize  the  separation 
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of  the  data  classes.  These  were: 

(1)  The  distance  between  the  highest  and  second  highest  peaks; 

» (2)  The  distance  between  the  highest  and  third  highest  peaks; 

(3)  The  distance  between  the  second  and  the  third  highest  peaks; 

(4)  The  distance  between  the  left  most  and  right  most  peaks; 

* 

(5)  The  ratio  of  the  distance  between  the  highest  and  left  most 
peaks  over  the  distance  between  the  left  most  and  right  most 
peaks ; 

(6)  The  ratio  of  distance  between  the  second  highest  and  the  left 
most  peaks  over  the  distance  between  the  left  most  and  right 

« most  peaks; 

(7)  The  ratio  of  the  amplitude  of  the  second  highest  over  the 
highest  peaks; 

(8)  The  ratio  of  the  amplitude  of  the  third  highest  over  the 
highest  peaks; 

(9)  The  number  of  peaks  above  10%  of  the  highest  peak; 

(10)  The  number  of  peaks  above  50%  of  the  highest  peak. 

The  primary  motive  behind  the  selection  of  these  Intuitive  features 
was  to  reduce  the  dimensionality  of  the  feature  space.  Again,  other  inherent 
advantages  for  this  step  will  be  discussed  later. 

Classifier  Selection 

Once  a data  base  has  been  acquired,  a set  of  potential  features  ex- 
tracted, and  a subset  evaluated  to  be  promising  as  input  for  decision  making, 
* the  actual  decision  rule  development  process  can  begin.  Usually,  the  first 


,1 

i 


UNCLASSIFIED 

145 


UNCLASSIFIED 


step  in  that  procedure  is  to  divide  the  input  data  sets  into  two  equal  parts 
which  should  be  as  identical  as  possible.  One  part  is  used  as  the  learning 
set,  the  other  as  the  test  set.  This  process  is  called  generalization;  how- 
ever, because  of  limited  amounts  of  available  data,  this  step  was  bypassed. 
Thus,  all  input  data  was  used  not  only  as  a learning  but  also  as  a test  set. 

The  second  step  is  to  select  the  classifier.  The  particular  choice 
depends  on  such  variables  as  (1)  the  distribution  of  vectors  in  feature  space, 
(2)  the  amount  of  time  required  to  process  unknown  "radar  signatures"  and, 
most  important,  (3)  the  performance  of  the  classifier.  For  this  reason  we 

investigated  three  basic  types  of  classifiers.  The  simplest  of  these  were 

->■ 

the  linear  processors.  Here  the  components  of  the  sample  vector  X are  multi- 
plied by  the  coefficients  such  that 
N 

S = 2 (38) 

i 

The  unknown  sample  X is  then  characterized  as  a mine  or  non-mine  depending  on 

whether  the  sum  S exceeds  or  falls  below  a threshold  S . The  W.  are  deter- 

o 1 

mined  during  the  training  (i.e.,  during  processing  of  the  learning  set). 

One  of  these  linear  processors  is  the  Fisher's  Linear  Discriminant 

(FLD)  classifier.  The  FLD  projects  the  feature  vectors  in  multi-dimensional 

->■ 

space  onto  a vector  W such  that  the  scatter  between  the  means  of  the  distri- 
butions is  maximized  and  the  scatter  within  each  distribution  is  minimized. 

* 

Each  projected  vector  defines  a point  on  W.  Thus  if  sufficient  data  is 

->• 

available,  a histogram  can  be  constructed  along  W.  The  distributions  of  the 
classes  can  then  be  approximated  by  probability  density  functions.  This 
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process  is  illustrated  in  two  dimensions,  a and  b,  in  Figure  74.  The  circles 

and  squares  represent  the  tips  of  feature  vectors  for  two  different  classes 

-►  -*■ 

such  as  a mine  and  a non-mine.  The  dashed  lines  and  denote  two  differ- 
->  -► 

ent  choices  of  the  W vector.  The  one  marked  as  is  shown  as  the  solution 
vector.  Along  this  direction,  the  projected  points  form  two  density  distri- 
butions which  separate  into  their  respective  classes;  however,  the  vector 
marked  clearly  cannot  be  the  propor  choice  because  the  projected  points 
overlap  to  a considerable  extent.  To  obtain  the  classifier,  the  distributions 
were  separated  at  the  point  S^.  Hence,  to  classify  an  unknown  signature, 
one  must  first  extract  the  feature  vector  X,  then  multiply  it  by  the  vector 

W and  evaluate  the  sum  S in  Equation  (38)  to  check  whether  S is  greater  than 

-> 

or  less  than  S^.  The  multiplication  projects  the  vector  X onto  W.  Any  point 

which  falls  below  S is  classified  as  a mine,  for  example,  whereas  points 
o 

above  S would  be  labeled  as  non-mines. 


o 

A classifier  using  a threshold  such  as  the  point  is  also  known  as 
the  Bayes  classifier.  To  further  illustrate  this  type  of  classifier,  con- 
sider the  two-category  case  shown  in  Figure  75.  The  two  continuous  distribu- 

-► 

tions  are  similar  in  appearance  to  histograms  along  W^,  but  differ  from  his- 
togram representations  by  a scale  factor  such  that 


p(w/xj^)  dw  = 1 


pCw/x^)  dw  = 1 


(39) 


(40) 


In  other  words,  they  are  probability  density  functions.  The  function  p(w/xj^) 
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denotes  the  normalized  density  of  mine  vectors  projected  onto  at  any 
point  w.  The  normalized  density  for  non-mines  at  the  same  point  Is  given 
by  p(w/x2) . 

Consider  the  two  functions  being  divided  Into  two  regions , and  R2 
at  9.  There  are  two  ways  In  which  a classifier  error  can  occur:  either  a 
vector  falls  Into  R^^  and  the  true  state  Is  X2»  or  It  falls  Into  region  R2 
and  the  true  state  Is  x^^.  Therefore,  the  probability  of  error  Is  given  by 


P (error) 


^Xj^)  P(Xj^)  dw  + J*  pCw/x^)  PCx^)  dw  (41) 


where  P(Xj^)  and  P(x2)  are  the  a priori  probabilities  of  the  two  cases.  For 
Instance,  If  the  learning  set  for  the  mines  would  contain  100  signatures  and 
that  for  the  non-mines  only  50  signatures,  then  P(Xj^)  and  P(x2)  would  be 
approximated  by  values  2/3  and  1/3,  respectively. 

The  two  terms  In  Equation  41  are  merely  the  areas  under  the  tails  of 
the  functions  (shaded  portions).  Because  6,  and  hence,  the  regions  R^^  and 
R2,  were  selected  arbitrarily,  the  probability  ot  error  Is  not  as  small  as 
It  might  be.  However,  the  minimum  error  can  be  found  by  differentiating 
P (error)  and  setting  the  result  equal  to  zero  such  that 


— - 0 (^2) 

a 9 

The  solution  to  Equation  (42)  Is  the  necessary  condition  required  to  minimize 
P (error)  and  Is  given  by 

p(9/xj^)  P(Xj^)  = p(9/x2)  P(x2)  (43) 
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Equation  (A3)  shows  that  P(error)  is  a mlnimuTn  when  0 is  at  S , l.e.,  the 

o 

Rayes  position. 

We  used  the  FLD  not  only  as  a classifier  but  also  to  select  those 
intuitive  features  (i.e.,  feature  evaluator)  which  had  the  greatest  influ- 
ence on  the  separation  of  the  projected  distributions.  Clearly,  if  one  of 
the  components  is  zero  or  much  smaller  than  others,  the  intuitive  feature 
weighted  by  that  component  does  not  contribute  significantly  to  the  separa- 
tion of  the  projected  classes  and,  thus,  was  discarded.  In  this  manner  we 
were  able  to  reduce  the  dimensionality  of  our  feature  vectors  yet  retain 
a selected  subset  of  intuitive  features  which  did  not  alter  the  classifier 
performance  significantly.  That  subset  was  shown  previously  in  the  para- 
graph on  feature  selection. 

Unfortunately,  the  FLD  does  not  perform  well  with  just  any  distribu- 
tion of  feature  vectors.  For  example,  if  the  cluster  of  mine  data  were 
contained  within  the  non-mine  data,  or  if  the  two  distributions  were  to 
overlap  to  a significant  extent,  it  is  expected  that  the  1,  3,  5,  or  7- 
nearest  neighbor  (NN)  classifier  will  give  better  performances.  With  the 
nearest  neighbor  classifier,  all  test  data  was  first  stored  in  memory. 

The  vector  distances  were  then  measured  between  a sample  of  unknown  origin 
and  its  neighbors,  and  the  unknown  was  classified  into  that  group  containing 
the  largest  number  of  the  k-nearest  samples. 

The  NN  algorithm  works  best  with  data  in  which  the  same  scale  is  used 
to  measure  the  amplitude  for  all  the  feature  vector  dimensions.  If  the 
various  features  are  measured  on  different  scales,  then  the  definition  of 
distance  becomes  a problem  because  distance  between  points  in  one  feature 
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dimension  may  be  measured  on  a scale  so  radically  different  from  the  others 
that  the  contributions  of  the  other  dimensions  to  the  total  distance  may 
become  negligible.  Thus,  information  from  the  other  features  which  might 
be  useful  In  separating  data  classes  Is  effectively  lost.  The  distance  mea- 
sure Is  not  a problem  with  feature  vectors  which  were  derived  from  raw  data, 
but  may  become  one  when  intuitive  features  are  used  to  form  the  data  vectors. 

Another  disadvantage  with  the  NN  classifier  is  the  fact  that  storage 
of  the  learning  set  requires  a large  memory.  Furthermore,  large  processing 
times  are  needed  since  all  distances  between  the  unknown  point  and  every 
other  point  in  the  learning  set  must  be  calculated.  This  fact  could  influ- 
ence the  search  speed  of  the  mine  detection  vehicle.  Thus,  we  turned  our 
attention  to  a third  classifier,  the  space  partition  (SP)  method. 

Points  belonging  to  different  classes  may  be  separated  directly  in 
feature  space  by  means  of  a decision  or  hyper-surface  (i.e.,  surface  in  n- 
dimensional  space).  A quadratic  surface  is  defined  by  the  quadratic  dis- 
criminant  function  G(X)  = 0 where 


-y  d d d 

G(X)  = w + 2 W X . + 2 2 w_x.x, 

O i jilij 


and  d is  the  number  of  dimensions.  The  X^  are  the  components  of  a feature 
vector  and  W^,  W^,  and  W^.  are  components  of  a weight  vector.  The  quadratic 
discriminant  function  G(X)  is  not  equal  to  zero  if  the  vector  X is  either 
above  or  below  the  hyperplane  as  shown  in  Figure  76.  The  figure  was  simpli- 
fied to  show  the  hyperplane  in  a two-dimensional  space.  In  this  case,  the 


hyperplane  G(x)  = 0 marks  the  boundary  between  two  different  classes  of  data. 
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Figure  76.  Space  partition  classifier. 
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Vectors  (only  the  tip  of  each  vector  is  indicated  as  a circle  or  square) 


which  fall  above  the  hyperplane  will  give  a value  for  G(X)  which  is  greater 


than  zero.  Below  the  hyperplane  G(X)  will  always  be  less  than  zero.  This 


fact  suggests  the  following  decision  rule  for  a two-category  case  such  as 


mine  vs . non-mine : 


(1)  Decide  class  1 if  G(X)  > 0 (l.e.,  points  on  the  positive  side  of 


the  decision  surface  are  labeled  class  1) ; 


(2)  Otherwise  (G(X)  < 0]  decide  class  2 (i.e. , points  on  the  nega- 


tive side  of  the  surface  are  labeled  class  2. 


The  actual  classification  process  is  similar  to  the  FLD  (in  fact. 


the  FLD  is  a special  case  of  the  SP  classifier)  except  that  the  vector  X is 


expanded  to  include  the  cross  products  x^  x^  for  a quadratic  surface,  such 


that  the  complete  vector  X is  given  by 


X = 1 + X, ...X  ...+  Xt^...X  ^...+  X,-...X  , X 

In  In  12  n-1  n 


The  number  1 is  needed  in  case  is  not  equal  to  zero.  All  components  are 


then  multiplied  by  their  respective  weights  W , W . , W. , and  the  product  G(X) 

o i IJ 


tested. 


To  design  a classifier  of  this  type,  the  weights  W^,  W must  be 


determined  to  give  a minimum  number  of  misclasslfied  samples.  Because  W is 


always  normal  to  the  decision  surface,  a reorientation  of  the  vector  W is 


equivalent  to  reorienting  the  decision  surface.  There  are  several  methods 


of  arriving  at  the  correct  W.  For  this  report,  the  perceptron  criterion 


function  was  used.  The  first  step  was  to  generate  a transformation  of  fea- 


ture vectors  which  simplified  the  treatment  of  the  two-category  case,  i.e.. 
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we  set  Y = X If  X is  in  data  class  1 and  Y ■ -X  if  X is  in  data  class  2. 

-►  -v 

This  results  in  G(X)  > 0 for  all  properly  classified  data  points  but  G(X) < 0 

for  mlsclasslfled  points.  Consider  Y' (W)  to  be  the  set  of  samples  mlsclassi- 

-►  -► 

fled  by  weight  vector  W.  The  perceptron  criterion  function  J (W)  which  is 
-►  P 

minimized  when  W is  a solution  vector,  may  then  be  expressed  as 

J (W)  - 2 (Wy")  (46) 

P y. 

where  the  Y™  are  the  feature  vectors  of  misclassified  samples.  Since  WY  < 0 
for  misclassified  samples,  J will  always  be  greater  or  equal  to  zero.  More- 

P 

over,  J will  be  zero  only  if  W perfectly  separates  the  data  into  their  re- 

P ->■ 

spective  groups.  In  order  to  calculate  W which  minimizes  this  function,  one 
must  first  solve  for  the  gradient  of  •Ip(W)  given  by 
-*■ 

VJ  (W)  » 2 (-Y®)  (47) 

P y. 

->■  -*■ 

To  start  the  W search  sequence,  an  arbitrary  is  introduced  initially. 

This  vector  is  then  moved  toward  the  proper  solution  by  the  gradient  descent 

-+•  -*■ 

algorithm  which  states  that  the  next  value  of  W,  is  obtained  by  moving 

some  distance  from  the  initial  starting  vector  Wj^  in  the  direction  of  the 
steepest  descent  along  (i.e.,  along  the  negative  of  the  gradient 


Vi  = \ \ 


\+l  \ ^ ^ ^ 

where  k > 1.  Thus,  each  vector  Y is  first  tested  for  misclassif ication 
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[i.e.,  G(Y)  < 0].  If  a mis class If led  vector  Y™  is  found,  a fraction  of  Y™, 

-*■  -*■ 

Y™,  la  added  to  the  initial  weight  vector  W,  to  reorient  the  decision 
surface.  Thereafter,  the  search  for  mis classified  samples  and  reorientation 

of  the  hyperplane  continues  until  a position  is  found  which  results  in  the 

->■ 

minimum  number  of  misclasslf led  samples.  The  particular  W corresponding  to 

that  position  is  the  solution  vector.  The  variable,  was  usually  set  to 

decrease  as  1/k  for  convergence  toward  a solution. 

There  are  several  disadvantages  associated  with  this  technique.  It 

was  previously  pointed  out  that  the  quadratic  discriminant  function  uses  an 

expanded  vector  which  Includes  the  cross-terms  X.  X . Thus,  if  the  dimension 

->  J -»• 

size  of  the  original  vector  X is  d then  the  expanded  vector  Y will  have  diraen- 

sions  d equal  to  [(d+1)  (d+2)]/2.  For  example,  for  10  (d  = 10)  intuitive 
A 

features,  d will  be  equal  to  66.  It  is  important  that  each  class  contain 
approximately  several  times  more  samples  than  dimensions  (for  quadratic  func- 
tions the  number  of  samples  should  be  » 60  if  d = 10) . Otherwise  the  vector 
.*■ 

W might  converge  to  a solution  which  is  underdetermined.  In  other  words, 
with  that  W vector  it  is  likely  that  "signatures"  which  were  not  used  for 
the  learning  set  would  be  misclassifled.  Of  course,  the  nundier  of  extracted 


intuitive  features  can  be  reduced,  for  example  to  two,  such  that  d = 6. 
However,  the  amount  of  information  contained  in  each  feature  vector  will 
suffer  correspondingly  probably  resulting  in  Increased  overlap  of  the  dis- 
tribution. 
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CORRELATION  PROCESSING 

In  addition  to  conventional  pattern  recognition  techniques  previously 
discussed,  another  technique  was  investigated  for  discrimination/classi- 
fication. This  technique  involved  correlating  a reference  waveform  with 
actual  data  and  submitting  the  results  to  a threshold  for  a decision. 

The  reference  waveform  is  one  that  is  mathematically  derived  and  used  to 
correlate  with  the  measured  data.  This  process  is  more  accurately  a cross- 
correlation of  an  idealistic  waveform  with  actually  measured  waveforms. 

As  discussed  in  Section  IV  of  this  report,  the  mathematically  derived  tem- 
poral responses  of  mines  buried  in  soil  correspond  well  with  the  temporal 
responses  of  actual  measured  data.  Thus,  it  is  expected  that  the  reference 
waveforms  so  generated  will  look  very  much  like  the  waveform  that  is  found 
in  the  actual  measured  data  and  which  represents  the  return  from  a non-metallic 
mine  buried  in  soil. 

The  details  of  the  cross-correlation  scheme  are  shown  in  Figure  77. 

This  technique  was  only  exercised  for  the  data  taken  by  the  National  Bureau 
of  Standards,  referred  to  as  NBS  data.  Thus,  in  Figure  77,  it  is  assumed 
that  the  NBS  data  which  exists  in  the  frequency  domain  is  used  in  this  pro- 
cess. This  data  is  transformed  by  the  Fast  Fourier  transform  into  the  time 
domain.  The  portion  of  the  temporal  data  corresponding  to  the  antenna  coup- 
ling and  the  ground  return  is  gated  out  and  only  that  portion  from  below  the 
ground  surface,  including  the  non-metallic  mine,  is  left.  The  position  of 
the  gate  is  determined  automatically  by  first  locating  the  most  negative  peak 
of  the  ground  return  and  then  setting  the  gate,  or  zeroing  out  the  data,  a 
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Figure  77,  Cross-correlation  processor. 
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certain  number  of  san5)les  beyond  that  negative  peak,  thus  guaranteeing  that 
the  ground  return  is  indeed  removed  from  the  data  before  correlation.  The 
gated  temporal  response  is  then  transformed  back  to  the  frequency  domain  by 
using  the  inverse  Fast  Fourier  transform.  The  resulting  spectral  response 
represents  not  only  the  response  of  the  mine  return,  but  the  response  of  any 
ocher  interference  which  exists  below  the  ground  surface. 

The  reference  spectral  waveform,  which  is  mathematically  derived  using 
the  mathematical  model  described  earlier  in  this  report,  is  multiplied  with 
the  gated  measured  spectral  response.  Both  spectral  responses  are  sampled  in 
the  frequency  domain  and  thus  multiplication  takes  place  samply  by  sample — 
that  is,  the  first  sample  of  the  math  model  spectral  waveform  is  multiplied 
by  the  first  sample  of  the  gated  spectral  response,  and  so  forth,  on  through 
the  point  of  each  waveform.  The  product  of  these  two  responses  is  then 
transformed  by  the  Fast  Fourier  Transform  back  into  the  time  domain.  The 
result  of  this  transformation  represents  the  cross-correlation  of  the  gated 
temporal  response  with  the  mathematically  derived  temporal  waveform. 

Each  sample  in  the  time  domain  of  the  cross-correlation  function  is 
subjected  to  a threshold  which  has  been  normalized  to  background  level.  If 
any  of  the  samples  exceed  this  threshold,  it  is  decided  that  the  measured 
data  response  matches  that  of  the  math  model  and  thus  represents  the  return 
from  a particular  type  of  non-metallic  mine  which  is  buried  in  the  soil. 

Thus,  discrimination/classification  is  accomplished  by  this  thresholding 
process.  If  the  measured  data  does  not  include  returns  from  the  mine  which 
the  math  model  represents,  then  the  cross-correlation  function  will  not 

UNCLASSIFIED 

158 


k 


I 

1 


♦ 


UNCLASSIFIED 


exceed  normalized  threshold  and  the  decision  is  made  that  the  particular 
type  of  mine  does  not  exist  in  that  measured  data. 

If  the  measured  data  exists  as  a temporal  response  rather  than  a 
frequency  response,  the  first  Fast  Fourier  Transform  shown  in  Figure  77  is 
not  necessary;  the  measured  temporal  data  is  simply  gated  and  then  subjected 
to  the  remainder  of  the  process . 

To  illustrate  the  correlation  function  which  is  expected  at  the  out- 
put of  the  processor  when  the  measured  data  and  the  reference  waveform  are 
perfectly  matched,  a mathematically  derived  temporal  signature  of  a mine, 
shown  in  Figure  78(a),  was  correlated  with  itself  to  yield  the  autocorrela- 
tion function  shown  in  Figure  78(b).  Now  consider  the  cross-correlation 
of  this  ideal  mine  signature,  previously  referred  to  as  the  math  model 
reference,  with  actual  measured  data  as  shown  in  Figure  79.  The  measured 
data  is  the  gated  temporal  response  from  a type  A mine  buried  6 inches  deep 
in  soil  with  7%  moisture.  Figure  79(b)  shows  the  result  of  the  cross- 
correlation. Note  that  the  correlation  waveform  in  Figure  79(b)  compares 
very  favorably  with  the  idealistic  waveform  shown  in  Figure  78(b).  As  a 
comparison,  the  cross-correlation  of  the  ideal  mine  signature  with  the 
return  from  a metal  plate  buried  6 inches  deep  in  7%  moisture  soil  is  shown 
in  Figures  80(a)  and  80(b).  Not  only  does  the  cross-correlation  function 
appear  different  in  structure  but  the  amplitude  of  the  cross-correlation 
function  is  significantly  smaller  than  that  in  Figure  79(b). 

When  the  cross-correlation  function  is  normalized  to  the  measured  tem- 
poral response,  the  amplitude  of  the  metal  plate  correlation  function  is 
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Figure  79 . Cross-correlation  of  ideal  mine  signature  with  measured 
data  from  a type  A mine. 
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a.)  Measured  data,  GT-68,  metal  plate,  6"  depth,  7%  moisture. 


b.)  Cross-correlation  of  reference  signal  and  measured  data 
GT5-68. 

Figure  80.  Cross-correlation  of  Ideal  mine  signature  with  measured 
data  from  a metal  plate. 
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approximately  half  that  of  the  mine  correlation  function.  Thus,  it  is 
obvious  that  in  these  two  examples  a threshold  can  be  used  to  differentiate 
between  the  metal  plate  and  the  type  A mine.  Non-metallic  mines  can  also 
be  differentiated  from  other  buried  objects  such  as  rocks  and  roots  as  well. 

* In  addition,  the  potential  for  differentiating  between  various  types  of  mines 
is  possible  if  the  thickness  of  the  various  mines  differ  significantly. 

The  performance  of  the  cross-correlation  processor  is  indicated  in 
Figures  81  through  86  in  the  form  of  probability  of  detection  versus  prob- 
ability of  false  alarm  curves.  In  Figure  81,  the  performance  of  the  pro- 
cessor is  shown  for  a buried  type  A mine.  The  sub-headings  at  the  top  of 

* the  figure  indicate  that  the  data  base  from  which  these  probabilities  were 
calculated  was  obtained  from  type  A mines  burled  at  both  3 inches  and  6 

, inches  below  the  soil  and  is  based  on  measurements  made  directly  over  the 

mine  and  4 Inches  either  side  of  the  center  of  the  mine.  Also,  the  data 
base  consists  of  soils  with  moisture  contents  of  7,  17,  and  12-20  and  13- 
16%.  False  alarm  is  defined  in  this  figure  as  the  declaration  of  a target, 
namely  a type  A mine,  when  in  fact  there  is  no  mine  buried  in  that  location. 

The  false  alarm  is  likely  generated  by  anomalies  within  the  soil  below  the 
surface  of  the  ground  which,  as  indicated  in  the  figure  sub-headings,  are 
called  background  noise. 

, As  indicated  in  Figure  81,  the  detection  performance  of  the  cross- 

: 

correlation  processor  under  these  particular  circumstances  is  very  poor. 

For  example,  a probability  of  detection  of  90%  would  yield  approximately 

* 

an  80%  chance  of  a false  alarm  when  what  is  desired  is  a false  alarm  prob- 
ability of  something  on  the  order  of  1%  or  less.  It  should  be  noted  that 
the  solid  curve  in  Figure  81  is  a quadratic  mean  square  fit  to  the  data 
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Figure  82.  Cross-correlation  classifier. 
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Figure  83.  Cross-correlation  classifier. 
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Figure  85.  Cross-correlation  classifier 
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points  which  are  represented  by  triangles  within  the  plot. 

The  detection  performance  of  the  correlation  processor  with  a buried  • 

type  B mine  is  shown  in  Figure  82.  This  mine  was  somehow  different  from  the 
type  A and  the  detection  performance  is  likewise  different.  In  general,  the 
detection  performance  is  somewhat  poorer  than  that  for  the  type  A mine.  The 
detection  performance  of  the  cross-correlation  processor  for  another  type  of 
non-metalllc  mine,  the  type  C mine,  as  compared  to  background  noise  from  the  i 

subsurface  soil,  is  shown  in  Figure  83.  It  compares  roughly  to  the  performance  j 

i 

of  the  type  A and  B mines  under  the  same  clrciunstances . j 

i 

Figure  84  shows  the  detection  performance  of  the  cross-correlation  ^ 

processor  where  the  type  A mine  is  compared  to  a metal  plate,  two  sizes  of  ^ | 

i 

rocks,  and  a root  which  are  buried  in  the  soil.  In  this  case,  a false  alarm 
is  defined  as  declaring  a mine  present  when  in  fact  the  metal  plate  or  the 
rocks  or  the  root  are  present.  As  will  be  noted,  the  detection  probability  is 
somewhat  better  than  that  for  the  type  A mine  in  background  noise,  but  not 
much.  The  detection  performance  for  the  type  B mine  as  compared  to  the  metal 
plate,  the  rocks,  and  the  root  is  shown  in  Figure  85;  the  performance  is  very 
poor.  Figure  86  likewise  indicates  the  performance  of  the  processor  when  the 
target  is  a type  C mine  and  the  non-mines  consist  of  the  plate,  the  two  rocks, 
and  the  root.  Performance  here  is  likewise  very  poor. 

The  conclusions  drawn  based  on  the  performance  curves  just  presented 
are  that  the  cross-correlation  processor  has  a very  poor  classification/ 
discrimination  performance  and  does  not  provide  the  performance  required 
for  the  mine  detection  and  classification  system.  After  considerable  investi- 
gation into  the  reasons  for  the  poor  performance,  it  has  been  determined  that 
the  mathematical  model  reference  waveform  does  not  sufficiently  represent  the 
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waveform  that  is  to  be  expected  from  the  mine  return.  This  is  because  the 
waveform  signature  involves  variables  which  are  dependent  upon  soil  moisture 
and  dielectric  constants  of  the  soil.  As  Indicated  previously,  if  the  soil 
moisture  at  the  top  of  the  mine  is  different  from  that  at  the  bottom  of  the 
mine,  the  radar  return  from  the  bottom  of  the  mine  will  be  of  a different 
amplitude  than  that  of  the  top.  However,  the  mathematical  reference  wave- 
form assumes  that  same  amplitude  return  from  both  the  bottom  and  top  of  the 
mine.  Thus,  in  general,  a mismatch  will  occur  between  the  reference  wave- 
form and  the  actual  returning  waveform.  Thus,  the  correlation  between  the  two 
will  not  have  as  large  an  amplitude  as  it  would  if  the  reference  waveform 
were  exactly  matched  to  the  returning  waveform.  The  actual  returns  are  signi- 
ficantly different  from  the  reference  waveform  as  indicated  in  the  discussion 
of  Figure  80. 

Another  type  of  correlation  was  considered  during  this  study  and  is 
illustrated  in  Figures  87  and  88.  This  correlation  scheme  is  called  the 
impulse  correlation  and  is  illustrated  in  Figure  87  with  an  impulse  corre- 
lation circuit  diagram.  The  concept  is  based  on  the  fact  that  ideally  what 
is  desired  is  that  the  reference  temporal  waveform  illustrated  in  Figure  78(a) 
be  correlated  with  measured  data.  A close  look  at  the  reference  waveform 
indicated  in  Figure  78(a)  indicates  that  is  consists  of  two  major  peaks: 
first,  a positive  going  peak  and  then  a negative  going  peak  separated  in 
time  by  approximately  the  thickness  of  the  mine.  A approximation  to  this 
waveform  would  be  a waveform  which  consists  of  two  impulses  as  shown  in 
Figure  88.  One  of  the  impulses  is  negative  and  one  is  positive  and  they  are 
separated  by  a time,  t^,  which  is  approximately  the  time  duration  of  an 

UNCLASSIFIED 

171 


Figure  87.  Impulse-correlation  circuit. 
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Figure  88.  Impulse  response  of  impulse-correlation 
circuit . 
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electromagnetic  wave  traveling  through  a non-metallic  mine.  The  impulse 
correlation  circuit  shown  in  Figure  87  takes  the  measured  data  and  a delayed 
version  of  the  measured  data  and  subtracts  the  two.  The  output  of  this 
circuit  is  equivalent  to  the  correlation  of  the  reference  waveform  approxi- 
mation with  the  actual  data.  The  impulse  response  of  the  circuit  is  what 
is  shown  in  Figure  88.  The  output  response  of  the  impulse  correlation 
circuit  is  given  by  the  convolution  of  the  impulse  response  shown  in  Figure 
88  and  the  input  signal,  s(t).  The  output  then  is  equivalent  to  a signal 
which  is  derived  by  convolving  the  reference  waveform  approximation  with 
the  input  signal,  s(t).  Thus,  the  output  is  a close  approximation  to  the 
correlation  of  the  reference  waveform  and  the  input  data.  The  advantage 
of  this  approach  is  that  less  complex  processing  is  required  as  compared 
to  the  cross-correlation  processor  shown  in  Figure  77.  The  detection  per- 
formance is  found  to  be  comparable  to  that  of  the  cross-correlation  proces- 
sor. Detection  performance  curves  are  given  in  Figures  89  through  94.  A 
comparison  of  this  set  of  curves  with  those  performance  curves  given  in 
Figures  81  through  86  for  the  cross-correlation  processor  will  indicate  a 
very  close  agreement.  Thus,  one  concludes  that  the  impulse-correlation 
processor  yields  roughly  comparable  performance  to  that  of  the  cross- 
correlation processor  and  is,  therefore,  unacceptable  for  the  intended  appli- 
cation. 
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Figure  92.  Impulse-correlation  classifier. 
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SECTION  VI 


DISCRIMINATION/ CLASSIFICATION  PERFORMANCE  AND  COMPARISON 


Discrimination  is  defined  as  the  process  of  deciding  if  a potential 
target  Is  Indeed  a target  (mine)  or  some  other  target-like  object  such  as 
a rock,  root,  piece  of  metal,  etc.  On  the  other  hand,  classification  is 
the  process  of  deciding  to  which  class  of  mines  a particular  target  belongs. 
The  performance  of  both  the  discrimination  and  classification  processors 
considered  in  Section  V are  presented  in  this  section  with  a discussion  of 
how  the  various  processing  algorithms  compare. 

First,  a general  discussion  of  the  format  in  which  the  performance 
curves  are  presented  is  discussed.  Limits  of  performance  are  also  considered 
to  give  the  reader  insight  into  the  measured  results.  Then  the  performance 
curves  for  the  NBS  data  and  the  MERADCOM  short-pulse  radar  data  are  discussed. 
PERFORMANCE  CURVES 

The  curves  in  Figures  96,  98,  and  100  show  the  probability  of  detec- 
tion versus  probability  of  false  alarm  for  a range  of  postulated  discrimina- 
tor processors  whose  performance  ranges  from  ideal  (perfect)  to  very  poor 
(worst  case).  Here  the  probability  of  detection,  Pj^,  is  defined  as  the 
probability  of  having  the  discriminator  decide  "mine"  for  a return  resulting 
from  a type  A,  B,  or  C mine;  the  probability  of  false  alarm,  Pp^^.  is  defined 
as  the  probability  of  having  the  discriminator  decide  "mine"  for  returns 
from  non-mine  targets . 
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The  performance  curves  indicate  not  only  the  relevant  probabilities 
but  also  convey  some  information  about  the  shape  and  overlap  of  the  distri- 
bution of  points  in  feature  space.  This  information  can  be  an  important 
asset  since  the  points  are  contained  in  N-dimensional  space,  and  thus  the 
shape  of  the  distribution  cannot  be  observed  directly.  In  order  to  demon- 
strate the  meaning  of  the  performance  curve  only  two  dimensions  will  be 
used.  Let  the  mines  and  non-mines  be  represented  by  two  continuous  density 
functions  as  shown  in  Figures  95,  97,  and  99.  Then,  the  corresponding 
probability  curves.  Figures  96,  98,  and  100  are  generated  by  sweeping  the 
decision  line  from  right  to  left  and  by  integrating  the  area  under  each 
probability  density  function  on  the  right  hand  side  up  to  S^.  For  widely 
separated  density  functions  as  in  Figure  95 , the  corresponding  probability 
curve  is  the  dashed  line  shown  in  Figure  96 . This  line  is  generated  in  the 
following  way:  when  S is  placed  at  point  S^,  the  Integrated  area  on  the 
right  is  zero  for  both  non-mines  and  mines  so  and  at  are  both 
zero.  When  S is  moved  to  S^,  approximately  1/2  of  the  total  area  under  the 
density  function  for  mines  will  be  summed  whereas  the  area  for  the  non-mines 
will  still  be  zero.  Hence,  at  S„,  the  P„  and  P„. , become  0.5  and  0 respec- 
tively.  At  the  point  the  P^  will  now  be  1.0  but  the  P^,^  will  still  be  0. 
Hence  up  to  this  point,  the  dashed  line  extending  from  (P^,  P^,^  =0,  0)  to 
(Pp,  Pp^  = 1.0,  0)  was  generated.  At  the  integrated  area  for  the  mines 
is  still  the  total  area,  but  now  1/2  of  the  area  for  the  non-mines  must  also 
be  taken  into  account.  Thus,  the  P^  and  P^^  will  have  reached  a value  of  1.0 
and  0.5,  respectively.  Finally  at  both  P^  and  Pp^  will  be  equal  to  1.0. 
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If  the  density  functions  completely  overlap  as  in  Figure  97,  the  resulting 
probability  plot  is  the  dashed  straight  line  as  shown  in  Figure  98. 

Other  distributions  as  illustrated  in  Figure  99  results  in  performance  which 

r 

fall  into  the  shaded  region  as  in  Figure  100. 

The  performance  curves  for  a classification  processor  appear  somewhat 
different  since  one  is  interested  in  the  probability  of  properly  detecting 
one  class  of  targets  (mines)  versus  another.  Again,  decisions  must  be 
made  in  N-dimensional  space.  However,  for  simplicity  of  illustration,  the 
present  discussion  is  limited  to  two-dimensional  space.  The  probability 
of  class  1 versus  probability  of  class  2 performance  curves  are  generated 
in  much  the  same  manner  as  the  discriminator  processor  performance  curves. 

The  decision  line  is  now  swept  from  left  to  right  and  only  the  area 
under  the  class  1 density  function  on  the  right  side  of  the  decision  line 
is  integrated.  The  area  under  the  class  2 function  is  computed  on  the  left 
I side  of  S^.  The  resulting  performances  are  shown  in  Figures  102,  104,  and 

106.  To  indicate  the  correlation  between  the  distributions  of  density 

I 

functions  and  the  probability  curves,  the  probabilities  corresponding  to  the 
decision  line  positions  through  were  marked  in  Figures  101  and  102. 

' Thus,  if  the  decision  line  S is  moved  to  the  P^  for  class  1 will  be  1,0, 

1 but  the  Pp  for  class  2 will  be  0,  because  the  area  under  the  class  2 density 

function  on  the  left  hand  side  of  is  zero.  At  S2,  the  Pp  for  class  1 
will  be  1.0  again,  but  the  Pp  for  class  2 will  be  equal  to  0.5,  etc.  This 
performance  curve  is  for  a perfectly  performing  classifier. 


The  dashed  straight  line  shown  in  Figure  104  corresponds  to  completely 
overlapping  distributions  as  shown  in  Figure  103.  Distributions  which  partially 
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overlap  as  illustrated  in  Figure  105  give  probabilities  which  fall  into 
the  shaded  region  shown  in  Figure  106. 

The  foregoing  discussion  indicates  the  nature  and  form  of  the  per- 
formance curves  which  are  to  be  expected  when  actual  measured  data  is  used 
in  the  probability  calcula  ■.ions . It  must  be  emphasized  that  the  actual 
sample  distributions  were  not  approximated  with  continuous  density  functions. 
As  a result,  the  plots  consist  of  discrete  probability  points.  In  other 
words,  instead  of  integrating  the  area  under  a density  function,  the  indi- 
vidual samples  for  a particular  class  to  the  right  or  left  of  the  decision 
line  were  counted  and  the  sum  divided  by  the  total  number  of  samples  within 
the  respective  class.  Whenever  feasible,  a quadratic  polynomial  was  fitted 
to  the  probability  points. 

Because  of  the  limited  amount  of  data,  the  collection  of  signatures 
was  used  not  only  for  a learning  set  but  also  for  a test  set.  For  this 
reason,  two  methods  were  considered  to  estimate  the  performance  of  the 
learning  set.  The  FLD  and  SP  classifiers  utilized  only  those  samples  that 
were  used  to  design  the  classifier  resulting,  perhaps,  in  optimistic  esti- 
mates for  Pp  and  P^^  (i.e.,  the  estimate  for  P^  may  be  too  high  and  the 
estimate  for  P^,^  may  be  too  low) . The  "leave  one  out"  method  was  applied 
to  the  nearest  neighbor  classifier.  It  uses  n-1  of  n available  data  samples 
to  design  a classifier,  then  classifies  the  left-out  sample.  This  process 
is  repeated  n times  to  give  estimates  for  P^  and  This  method  generally 

gives  estimates  which  are  pessimistic  because  not  all  of  the  data  samples 
are  used  to  design  a given  classifier.  Both  estimates  tend  to  converge 
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toward  the  true  probabilities  as  the  number  of  samples  in  the  learning  set 
increases . 

NBS  DATA 

The  NBS  data  was  used  to  evaluate  algorithm  performance  for  the  condi- 
tions shown  in  Table  6 to  facilitate  comparison  of  performances  between  different 
classifiers,  features,  mines  and  non-mines.  The  numbers  within  the  matrix 
refer  to  the  figure  in  Appendix  C where  the  particular  plot  can  be  found. 
References  to  these  figures  in  this  section  will  Include  the  Appendix 
figure  numbers  in  parenthesis  for  cross-reference.  The  vertical  column 
on  the  left  hand  side  Indicates  the  type  of  classifier  and  the  corresponding 
algorithms  that  were  used.  The  top  row  shows  the  classes  that  were  used 
for  each  curve.  It  is  important  to  note  that  performance  curves  were  gene- 
rated not  only  for  all  mines  versus  non-mines  but  also  for  their  constitu- 
ents to  point  out  the  type  of  mine  which  gives  the  best  classification 
performance  and  the  type  of  non-mine  which  consistently  gives  high  false 
alarm  rates.  Furthermore,  a detailed  analysis  of  the  performance  of  indi- 
vidual non-mines  versus  mines  can  lead  to  identification  of  parameters  which 
influence  signatures  and  hence  to  identification  of  classification  problems 
and  suggestions  for  improvements. 

There  are  many  variables  that  can  influence  the  signatures  during 
the  collection  phase.  These  are  attenuation  of  signal,  shape  of  scatterer, 
and  radar  pulse  shape.  Others,  such  as  the  gating  routine,  can  distort  the 
signatures  during  data  processing.  Finally,  there  are  parameters  which  do 
not  act  on  the  signatures  but  become  Imporatnt  during  the  training  phase  of 
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the  learning  set.  These  are  the  size  of  the  learning  set  and  the  size  of 
the  feature  vectors.  All  the  variables  above  can  Influence,  either  directly 
or  indirectly,  classifier  performance.  Some  were  discussed  in  other  sections 
of  this  report  and  therefore  will  be  mentioned  only  briefly.  All  others  will 
be  examined  in  depth. 

Signal  Attenuation  and  Signature  Gating 

There  are  two  ways  an  attenuated  signal  can  affect  the  signature  of 
a mine  or  non-mine.  First,  if  the  signal-to-nolse  ratio  is  small,  but  a 
part  of  the  signal  is  clearly  visible  above  background,  this  situation 
could  result  in  Improper  gating  of  the  complete  signature.  Second,  if  the 
the  attenuation  is  strong  enough,  the  whole  signature  could  become  indistin- 
guishable from  background. 

The  latter  case  applies  to  NBS  data  taken  under  wet  soil  conditions. 
This  means  that  many  learning  sets  contain  signatures  of  buried  objects  which 
are  indistinguishable  from  background  as  far  as  classification  is  concerned. 
The  worst  cases  occurred  for  soils  with  30%  moisture  content;  these  were  not 
utilized  in  generating  the  learning  sets.  Still,  some  highly  attentuated 
signatures  collected  from  soils  of  lower  moisture  content  were  present  in 
the  learning  sets.  The  fact  that  these  signatures  were  retained  contributed 
to  a significant  increase  in  the  number  of  misclassified  mines  and  non-mines. 
Of  importance  here  is  the  distribution  of  feature  vectors  in  feature  space. 
The  "background"  signatures  (highly  attenuated  signatures)  in  this  case 
will  form  their  own  separate  distribution  from  that  of  the  other  two  classes 
because  the  numbers  that  constitute  the  features  of  these  signatures  are 


UNCLASSIFIED 

189 


i 


UNCLASSIFIED 


generally  much  smaller  compared  to  "non-background  signatures."  Thus,  there 

is  no  reason  to  assume  that  a particular  "background  signature"  would  be  * 

classified  properly.  In  fact,  it  is  likely  that  the  whole  cluster  is 

assigned  to  one  class  or  another  affecting  thereby  not  only  the  but 

also  the  probability  of  classification.  In  the  case  of  the  NN  classifier 

mlsclasslflcatlon  depends,  in  general,  on  the  distribution  of  feature  vectors 

within  the  cluster  of  "background  signatures"  and  whether  the  cluster  is 

closer  to  one  class  or  another.  Clearly,  detailed  Information  on  the  effects 

of  "background  signatures"  on  classification  performance  is  not  available. 

However,  the  point  to  be  emphasized  here  is  that  both  mine  versus  mine  and 
mine  versus  non-mine  classification  performances  are  affected.  ^ 

The  other  way  in  which  an  attenuated  signal  can  affect  classification 
performance  is  through  the  gating  procedure.  The  gating  of  signals  is  a 
convenient  method  by  which  the  "signature"  of  an  object  is  separated  from 
its  ground  return  in  the  time  domain.  The  gate  positions  are  set  automatically 
by  the  computer.  Basically,  the  search  sequence  "looks"  for  peaks  above  aver- 
age background  beyond  the  ground  return.  Therefore,  for  weak  signals,  the 
gates  lock  on  to  the  highest  peak  above  background.  However,  this  does  not 
mean  that  the  gated  time  domain  contains  the  actual  mine  or  non-mine  signa- 
ture. In  fact,  it  might  contain  only  background  or  only  a portion  of  the 
desired  return  which  would  cause  the  signature  in  the  frequency  domain  to  be 
distorted.  If  these  distorted  mine  signatures  are  misclassif led , then  im- 
proper gating  will  also  affect  classifier  performance. 
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Shape  of  Scatterer  and  Radar  Pulse  Width 

% 

Although  the  diameter  of  a scatterer  such  as  a mine  or  a rock  affects 
the  magnitude  of  the  return  and  thus  the  magnitude  of  the  signature,  it  does  I 

not  influence,  to  a significant  extent,  the  character  of  the  signature  (e.g., 

* shape  and  location  of  the  maxima  and  minima) . The  parameter  which  does  play 
a dominant  role  in  the  selection  of  maxima  and  minima  of  the  signatures  is 
the  thickness  of  the  object.  For  pulsewldths  much  smaller  than  the  scatterer' s 
thickness , two  returns  were  observed  provided  the  second  return  was  not 
appreciably  attenuated  within  the  scatterer.  One  return  originated  from 
the  top,  the  first  dielectric  interface;  the  other  from  the  bottom  or  the 
second  dielectric  interface.  Hence,  the  time  delay  between  the  returns  was 
directly  proportional  to  the  scatterer's  thickness.  Therefore,  those  objects 
which  give  approximately  the  same  time  delays  as  mines  will  be  improperly 
classified  as  mines.  Non-mines  are  easily  identified  if  they  have  time  delays 
that  differ  from  those  produced  by  mines. 

The  worst  offenders  in  the  non-mine  category  which  gave  consistently 
high  false-alarm  rates  were  the  rocks  and  roots.  For  example,  the  perform- 
ance shown  in  Figures  107  (C-59  mines  vs.  rocks)  and  108  (C-60  mines  vs. 
roots)  are  inferior  to  that  of  Figure  109  (C-61  mines  vs.  plates).  In  the 
mine  category,  it  is  the  type  A and  type  B which  could  not  be  readily  dlstin- 

* gulshed  from  the  non-mines  as  illustrated  in  Figures  110  (C-48),  111  (C-53) , 

and  112  (C-58) . However,  because  the  rocks  and  roots  were  specially  chosen 

^ to  be  approximately  equal  in  thickness  to  mines,  it  is  suggested  that  the 

classifier  may  exhibit  improved  performance  under  actual  field  tests  where 
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Figure  110.  Nearest  neighbor  classifier. 
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rocks  of  all  sizes  and  shapes  are  available.  Therefore,  the  performances 

given  here  represent  a lower  bound.  In  other  words,  the  P_.  Is  expected  to 

FA 

decrease  and  the  to  Increase  'mder  field  test  conditions. 

A similar  argument  applies  to  the  classification  performances  between 
the  type  C and  type  A mines  which  have  a thickness  ratio  of  1.47.  Regard- 
less of  the  classifier  used,  good  performances  were  consistently  observed 
for  these  mines;  however.  It  was  generally  difficult  to  distinguish  between 
the  different  types  of  type  A mines  because  of  their  similarity  In  thickness 
as  shown  In  Figures  113  (C-41),  114  (C-64) , and  115  (C-33) . Again  this 
clearly  Indicates  that  the  scatterer's  thickness  is  a vital  factor  in  classi- 
fication performance. 

Learning  Set  Size 

Although  the  NBS  data  set  contained  signatures  for  all  parameters  of 
interest,  not  enough  signatures  were  available  to  evaluate  all  of  the  classi- 
fiers that  were  originally  planned  to  be  considered  (i.e.,  the  sample  size 
was  limited).  This  fact  was  not  foreseen  because  the  sample  size  must  not 
only  contain  all  the  information  about  the  environment  but  depends  to  some 
extent  on  the  size  of  the  feature  vectors  used.  Unfortunately  the  selection 
of  features,  and  hence,  determination  of  feature  size,  occurs  after  the  data 
has  been  collected.  The  remedy  for  this  problem  is  to  collect  a preliminary 
and  a final  learning  set.  The  preliminary  set  allows  one  to  find  the  type 
and  number  of  features  required  so  that  the  size  of  the  final  learning  set 
can  then  be  determined. 

As  a result  of  the  limited  size  of  the  learning  set,  we  were  able  to 


use  only  feature  vectors  with  a small  number  of  components  (i.e.,  intuitive 
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Inverse  of  the  within-scatter  matrix  which  measures  the  distribution  of 
feature  vectors  with  respect  to  the  mean  within  a particular  class.  Quite 
often,  this  inverse  does  not  exist  if  the  number  of  samples,  n,  is  of  the 
order  of  or  less  than  the  number  of  dimensions  of  the  feature  vectors  (i.e.. 


the  matrix  is  singular).  Since  each  file,  mine  or  non-mine,  contained  18 
signatures  the  performance  of  the  FLD  with  raw  data  could  be  tested  only 
for  the  case  of  "all  mine  vs.  non  mine"  where  the  number  of  signatures 
exceeded  the  number  of  dimensions  (i.e.,  the  combined  signatures  for  all 
"mines  vs.  all  non-mines"  were  54  and  72  respectively,  whereas  the  number  of 
dimensions  was  equal  to  32) . 

In  the  case  of  the  SP  classifier  the  quadratic  discriminant  function 

A A 

required  an  expanded  feature  vector  d where  d = [ (d+2) (d+1) ] /2  = 561  and 
d = 32.  The  processing  of  a data  array  with  that  size  would  not  only  have 
been  time  consuming  and  prohibitively  expensive,  but  would  also  have  resulted 
in  a greatly  underdetermined  solution  vector  (i.e.,  under  actual  field  condi- 
tions most  unknown  signatures  might  have  been  misclassified) . 

The  nearest  neighbor  classifier  may  also  have  been  affected  by  the 
low  sample  size.  In  fact,  the  type  of  nearest  neighbor  algorithm  that  pro- 
vides the  best  performance  is  determined  by  the  sample  size  and  the  distri- 
bution of  feature  vectors.  For  nearly  separated  distributions,  little  differ- 
ence in  performance  for  the  3,  5,  or  7 algorithm  would  be  expected  if  the 
sample  size  were  much  larger  than  the  number  of  nearest  neighbors  to  be 
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tested  (e.g. , hundreds  of  samples).  On  the  other  hand,  if  the  sample  size 
is  small  (i.e.,  between  18  and  72  as  in  our  case)  then  the  performances  are 
expected  to  depend  to  some  degree  on  the  number  K.  The  reason  is  that  the 
classification  probability  must  decrease  to  0.5  as  K approaches  the  total 
number  of  samples,  n,  provided  the  sample  populations  for  each  class  are 
equal,  as  in  our  case.  This  type  of  behavior  can  be  observed  by  comparing 
Figures  116  (C-48)  through  119  (C-63)  with  Figures  120  (C-71)  through  123 
(C-74),  and  124  (C-79)  through  127  (C-82) . The  plots  show  that  the  3-nearest 
neighbor  classifier  (K  = 3)  is  clearly  superior  to  all  other  NN  algorithms. 
Influence  of  Feature  Size  on  Classification 

In  many  instances,  classifiers  using  raw  data  outperformed  those 
using  intuitive  observables.  For  example.  Figures  116  (C-48)  through  119 
(C-63)  show  better  discrimination  performance  than  Figures  128  (C-67) 
through  131  (C-70)  where  raw  data  was  replaced  by  intuitive  data.  This 
suggests  that  although  only  the  "best"  intuitive  features  were  picked  from 
the  original  package  the  amount  of  information  each  Intuitive  vector  con- 
tained was  evidently  less  than  that  contained  in  the  raw  data. 

Classifier  Performance 

To  save  time  and  money  some  performances  with  the  5-  and  7-nearest 
neighbor  algorithm  performances  indicated  in  Table  6 were  not  evaluated. 

In  order  to  reduce  the  number  of  plots  to  be  analyzed  with  the  nearest 
neighbor  classifier.  Columns  5,  10,  15,  and  20  in  Table  6 (where  mines  are 
compared  to  the  class  of  all  non-mines)  were  tested  first,  on  the  assumption 
that  classification  of  mines  versus  the  constituents  of  non-mines  would  give 
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Figure  116.  Nearest  neighbor  classifier. 
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Figure  117.  Nearest  neighbor  classifier. 
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Figure  119.  Nearest  neighbor  classifier. 

UNCLASSIFIED 


PROBABILITY  OF  DETECTI 


PROBABILITY  OF  DETECTION  (PD) 


UNCUSSIFIED 


NEAREST  NEIGHBOR  CLASSIFIER 

RAW  DRTfl 

NO.  NEAREST  NEIGHBORS  5 
REFERENCE  322 


NONIilNE 

TARGET 

DEPTH 

LOCATION 

MOISTURE 


PLATE  ROCK(l)  R0CK(2)  ROOT 

TYPE  B MINE 
3 6 

+ 4 0 -4 

7 17  12-20  13-16 


o 

o 


CM  ^ 

a' 


o 

o 

'^.00  0.20  o'.  40  o'.  60  o'.SO 

PROBABILITY  OF  FALSE  ALARM  (PFA) 


Figure  121.  Nearest  neighbor  classifier. 
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Figure  123.  Nearest  neighbor  classifier. 
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Figure  124.  Nearest  neighbor  classifier. 
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Figure  127.  Nearest  neighbor  classifier. 
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Figure  129.  Nearest  neighbor  classifier. 
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Figure  131.  Nearest  neighbor  classifier. 
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comparable  or  better  performances.  When  it  was  established  which  algorithm 
would  perform  best,  the  data  in  other  columns  for  that  algorithm  were  also 
tested. 

An  overall  comparison  of  the  performances  for  all  three  classifiers 
indicated  the  3-nearest  neighbor  algorithm  to  be  as  good  as,  or  more  often 
better  than,  other  classifier  techniques.  For  example,  when  the  probability 
of  detection  (P^)  reached  1.0  the  probability  of  false  alarm  all 

mines  versus  non-mines  was  0.6,  0.42,  and  0.58  for  the  FLD,  NN,  and  SP 
classifiers,  respectively,  as  shown  in  Figures  132  (C-36),  133  (C-63) , and 
134  (C-86) . The  superior  performance  of  the  NN  classifier  suggests  that  the 
distributions  for  the  two  classes  may  have  overlapped  to  a considerable  ex- 
tent in  feature  space.  Under  these  circumstances  it  would  have  been  Impossible 
to  completely  separate  the  clusters  with  the  FLD  or  SP  technique.  However, 
better  separation  migjit  have  been  achieved  by  using  a set  of  different 
feature  vectors  which  would  have  had  a different  distribution  in  feature 
space.  The  performance  of  the  NN  (K  =3)  classifier  for  each  of  the  three 
types  of  mines  against  the  non-mines  is  shown  in  Figures  135  (C-48)  through 

137  (C-58) . When  all  three  of  the  mines  are  combined  into  one  class  and 
compared  to  the  non-mines,  the  combined  performance  given  by  Figure  133  (C-63) 
is  obtained. 

When  compared  to  one  another,  the  three  mines  can  be  distinguished  from 
one  another  with  the  NN  (K  = 3)  classification  algorithm  as  shown  in  Figures 

138  (C-64)  through  140  (C-66) . Note  that  the  type  C mine  can  easily  be  dis- 
tinguished from  either  the  type  A or  B mines.  However,  the  type  A and  B 


♦ 


PROBABILITY  OF  DETECTION  (PD) 

.00  0.20  0..40  0,.60  0^.80  1^.00 


UNCLASSIFIED 


NERRLSf  NEIGHBOR  Cl  MGG  I I 1 1 R’ 

RPW  noTR 

NO.  MP.flRESl  NtIGIinORO  O 
REFERENCE  316 

NONMINE  F’LRIE  ROCK(l)  fU)CK(2  1 ROOT 

TRRGET  type  C MINE  TYPE  A MINE  TYPE  B MINE 

DEPTH  3 6 

LOCATION  +4  0 -4 

MOISTURE  7 17  12-20  13-16 


0.60  O.HO 

PROBABILITY  OF  FALSE  ALARM  (F'l 


1 .00 


PROBABILITY  OF  DETECTION  ( PD  I 


UNCLASSIFIED 

SPACE  PARTITION  CLASSIFIER 
(QUADRATIC  DISCRIMINANT) 

INTUITIVE  DATA 
VERSION  2 
REFERENCE  355 

NONMINE  PLATE  ROCK(l)  R0CK(2)  ROOT 

TARGET  type  C MINE  TYPE  A MINE  TYPE  B MINE 

DEPTH  3 6 

LOCATION  +4  0 -4 

MOISTURE  7 17  12-20  13-16 


o 

o 


^.00  0.20  0.40  0.60  0.80  I .00 

PROBABILITY  OF  FALSE  ALARM  ( PFA  ) 

Figure  134.  Space  partition  classifier. 
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Figure  135.  Nearest  neighbor  classifier. 
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Figure  136.  Nearest  neighbor  classifier. 
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Figure  137.  Nearest  neighbor  classifier. 

UNCLASSIFIED 


PROBABILITY  OF  CLASSIFICATION 


UNCLASSIFIED 


NEAREST  NEIGHBOR  CLASSIFIER 

RRW  DRTfi 
VERSION  3 
REFERENCE  280 


TflRGET2 

TflRGETl 

DEPTH 

LOCATION 

MOISTURE 


TYPE  B MINE 
TYPE  A MINE 
3 b 

+ 4 0 

7 17 


-4 

12-20 


13-16 


o 

a 


o 

CO 


o 

CD 

a o' 


CO 

to  o 

CC 

-J  o' 
CJ 


o 

CM 


o 

o 


'^.OO  0.20  CK40  o'.60  0-80 

PROBABILITY  OF  CLASSIFICATION 
CLASS  ONE 

Figure  138.  Nearest  neighbor  classifier. 
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mines  are  not  as  easily  separated.  This  result  is  not  surprising  when  it 

is  remembered  that  the  type  A and  B mines  are  the  same  mine  design  (thick-  ; 

I ’ 

ness)  with  some  minor  (unknown)  differences.  ; 

SHORT-PULSE  RADAR  DATA  « 

The  motivation  for  collecting  and  processing  short-pulse  radar  data 
was  to  verify  the  previous  computerized  analysis  which  was  based  on  CW 
measurements  over  a broad  range  of  frequencies  (NBS  data).  By  virtue  of 
the  wide  spectrum  of  the  transmitted  pulse,  the  short-pulse  radar  system 
represents  the  analog  hardware  implementation  of  the  FFT  technique  used  to 
synthesize  time  domain  data  with  the  NBS  data.  The  short-pulse  radar  tern- 
poral  data  was  processed  by  the  same  discrimination/classification  algorithms 
as  the  NBS  data  and  performance  curves  were  generated  for  comparison  with 
the  NBS  data  performance  plots. 

The  short  pulse  radar  data  was  used  to  evaluate  discriminator/classi- 
fier performance  for  the  conditions  shown  in  Table  7.  Each  entry  in  the  matrix 
refers  to  the  figure  number  in  Appendix  C where  the  performance  plots  can 
be  found.  In  addition,  when  one  of  these  figures  is  referred  to  in  this 
section,  the  Appendix  figure  number  is  included  and  enclosed  in  parenthesis 
for  reference  purposes. 

As  will  be  noted,  this  matrix  does  not  cover  classified  performances 
to  the  same  extent  as  the  previous  NBS  data  matrix.  A comparison  of  the 
NBS  and  short-pulse  radar  gated  spectral  responses  indicated  very  close 

agreement  in  spectral  signatures  for  the  same  type  of  mine.  This  fact  sug-  | 

I 

gested  that  the  performance  associated  with  the  NBS  and  the  short-pulse  radar  j 
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data  might  be  quite  similar.  Consequently,  to  conserve  our  resources,  we 
elected  to  use  only  the  NN  classification  algorithm  which  had  previously 
given  the  best  performance. 

The  short-pulse  radar  data  was  only  taken  in  dry  soil  comparable  to 
the  7%  moisture  soil  used  in  the  NBS  data  collection.  A limited  number  of 
non-mine  objects  were  included  in  the  data  collection — e.g.,  a rock  and  a 
metal  plate  buried  6 inches  in  dry  soil.  Soil-only  data  where  no  mine-like 
objects  were  present  was  used  to  exercise  the  classification  algorithm. 

The  various  combinations  of  mines  and  non-mines  analyzed  is  shown  in  Table  7, 

Short-Pulse  Radar  Data  Performance  Matrix. 

It  should  be  noted  that  the  MERADCOM  short-pulse  radar  data  contained  • 

fewer  measurements  or  signatures  than  the  NBS  data.  For  example,  the  short- 
pulse  data  consisted  of  12  type  C mine  signatures  as  compared  to  18  for  the 
NBS  data.  Consequently,  the  statistical  confidence  level  for  the  short- 
pulse  data  is  not  quite  as  good  as  for  the  NBS  data.  Nevertheless,  the 
short-pulse  data  does  allow  a relative  comparison  of  performance  for  the 
two  measurements  systems. 

In  the  performance  comparison,  first  consider  the  gated  spectral 
signatures  of  the  two  basic  classes  of  non-metallic  mines — the  type  A and  the 
type  C mine.  The  spectral  response  of  the  type  A mine  buried  in  dry  soil 
for  the  NBS  data  and  the  short-pulse  data  is  shown  in  Figures  141  and  142 
respectively.  Comparable  responses  for  the  type  C mine  buried  in  dry  soil 
are  shown  in  Figures  143  and  144.  Note  that  for  each  set  of  data,  the  type 
A mine  is  characterized  by  nulls  separated  by  about  750  MHz  or  approximately 
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Figure  141.  Spectral  and  temporal  responses  of  type  A mines  buried  in  dry 
soil  (7%  moisture  in  loam)  from  NBS  measurements.  Files  54, 

55,  and  56  represent  positions  -4,  0,  and  +4  inches  over  buried 
mine. 
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Figure  143.  Spectral  and  temporal  responses  of  type  C mines  buried  in  dry 

soil  (7%  moisture  in  loam)  from  NBS  measurements.  Files  62,  63, 
and  64  represent  measurement  positions  -4,  0,  and  +4  inches  over 
burled  mine. 
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Figure  144.  Spectral  and  temporal  responses  of  type  C mines  buried  in  dry 
soil  from  short-pulse  radar  measurements.  Files  46,  47,  and 
48  represent  measurement  positions  -4,  0,  and  +4  inches  over 
buried  mine. 


UNCUSSIFIED 


UNCLASSIFIED 


the  reciprocal  of  Its  thickness  time  extent.  Since  the  type  C mine  is  not 

* as  thick  as  the  type  A,  the  nulls  in  its  spectral  response  for  each  set  of 

* data  is  about  900  MHz . The  bandwidth  of  the  short  pulse  data  is  not  as 
large  as  that  of  the  NBS  data;  therefore,  the  spectral  responses  of  the  two 

^ data  sets  do  not  compare  beyond  about  1.8  GHz. 

The  spectral  responses  just  presented  are  but  several  examples  of 
many  such  responses  found  in  Appendix  C.  Reference  to  the  Appendix  will 
confirm  the  conclusion  that  the  spectral  signatures  of  the  NBS  and  short- 
pulse  data  for  the  same  type  of  target  are  in  very  close  agreement  when 
the  soil  conditions  are  similar.  Thus,  one  surmises  that  the  classification 

• algorithms  designed  for  the  NBS  data  should  perform  just  as  well  for  the 
short-pulse  data. 

Detection  probability  versus  false  alarm  curves  for  the  short-pulse 
data  processed  by  the  various  NN  algorithms  are  given  in  Appendix  C.  As 
indicated  previously,  each  of  the  mine,  non-mine  combinations  listed  in 
Table  7 has  a corresponding  detection  performance  curve  in  Appendix  C. 
Several  of  these  performance  curves  are  presented  here  for  consideration. 

In  comparing  the  detection  performance  curves  of  the  NBS  and  short- 
pulse  data,  the  reader  should  remember  that  the  NBS  data  curves  are  for  a 
variety  of  soil  moistures  ranging  from  7%  (relatively  dry  soil)  to  20% 

I (relatively  wet)  whereas  the  short-pulse  data  curves  are  for  dry  soil 

only  (comparable  to  the  7%  moisture  soil).  Thus,  one  expects  the  short- 
pulse  data  detection  curves  to  indicate  better  performance  than  the  NBS 

4 

• data  curves.  This  conclusion  is  verified  by  comparing  Figures  145  (C-59) 
and  146  (C-103)  which  represent  the  performance  of  the  NN  (K  “3)  algorithm 
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Figure  146.  Nearest  neighbor  classifier. 
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with  the  NBS  and  short-pulse  data  of  all-mines  versus  a mine-sized  rock. 
Similar  results  are  obtained  when  all-mines  versus  a metal  plate  are  pre- 
sented to  the  algorithm  as  shown  in  Figures  147  (C-61)  and  148  (C-104) . The 
short-pulse  measurements  yield  a better  discrimination  performance  than 
the  all-moisture  NBS  data. 

When  the  ability  of  the  algorithm  to  distinguish  between  mine  types 
is  considered,  the  short-pulse  measurement  again  yields  a better  perform- 
ance than  the  NBS  data  as  illustrated  by  the  curves  in  Figures  149  (C-65) 
and  150  (C-105).  In  fact,  reference  to  Figure  150  indicates  perfect  classi- 
fication performance.  However,  since  the  number  of  data  points  used  was 
quite  limited  (12  samples  for  each  mine  type),  the  performance  is  not 
likely  to  be  as  good  as  indicated  in  this  figure. 

The  conclusions  to  be  drawn  from  comparing  the  NBS  and  short-pulse 
radar  data  is  that  the  collection  of  buried  mine  data  via  a CW  measurement 
system  (such  as  the  NBS  data)  seems  to  be  equivalent  to  collecting  mine 
data  via  a short-pulse  radar  system.  Therefore,  it  is  proposed  that  all 
of  the  conclusions  drawn  from  the  NBS  data  measurement  program  apply  equally 
well  to  an  equivalent  set  of  data  taken  with  the  MERADCOM  short-pulse  radar. 
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HARDWARE  REQUIREMENTS 

The  mine  detection  system  that  will  achieve  the  objectives  listed  In 
Table  8 Is  presented  In  this  section.  The  emphasis  of  this  study  was  on 
the  signal  processing  required  with  short-pulse  radar  returns.  Thus,  the 
design  of  the  processing  required  to  achieve  these  objectives  Is  the  major 
subject  of  concern  In  this  discussion. 

The  radar  to  be  used  In  the  mine  detection  system  Is  assumed  to  be 
the  MERADCOM  short-pulse  radar  furnished  by  the  government.  A comparison 
of  the  processor  performance  with  the  NBS  data  and  the  short-pulse  radar 
data  Indicates  approximately  equivalent  performance.  Therefore,  It  Is 

proposed  that  the  hardware  requirements  for  the  short-pulse  radar  as  listed  j 

In  the  CALSPAN  report  No.  MA-5366-E-3  entitled  "Vehicular -Mounted  Mine 
Detector  Radar  Operational  Manual,"  June  1976,  are  sufficient  for  the  mine 
detection  system  under  consideration.  i 


Table  8 

MINE  DETECTION  SYSTEM  REQUIREMENTS 


Target  Depth 

3 to  12  Inches 

Target  Type 

Non-metalllc  anti-vehicle 

Number  of  Target  Classes 

3 

Instantaneous  Width  Coverage 

4 meters 

Vechlcle-Mounted  Radar  Speed 

10  km/hr 

False  Alarm  Rate  , . 

Less  than  1 per  1000  meters 

Target  Type  Identification 

Yes 

Recallbratlon  with  Change  In  Location 

‘ No 

Adjustment  In  Detection  Threshold  Setting 
With  Change  In  Soil  Moisture 

Yes 

Standard  Parts  Used  In  Design 

Yes 
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SIGNAL  PROCESSING  ALGORITHM 

In  recommending  a processor  design,  Che  maximum  performance  achieved 
with  the  processing  algorithm  considered  In  this  study  was  deemed  necessary 
to  meet  the  requirements  listed  In  Table  8.  This  Implies  that  the  target 
search  pattern  consists  of  subterranean  interrogations  on  a 1-foot  square 
grid.  That  is,  each  interrogation  is  made  within  one  foot  of  all  other  in- 
terrogations. In  addition,  the  processing  parameters  used  are  those 
assumed  in  the  performance  analysis  presented  elsewhere  in  this  report. 

The  flow  diagram  shown  in  Figure  151  presents  the  processing  algorithms 
suggested  for  the  best  system  performance.  The  first  part  of  the  diagram 
outlines  the  detection  process  presented  earlier  in  this  report.  It  should 
be  remembered  that  internal  to  the  radar,  the  radar  returns  are  processed 
by  special  circuitry  and  thus,  the  bandwidth  of  the  radar  data  is  reduced 
by  10^  to  1.  Thus,  in  this  discussion  when  we  talk  about  a video  signal 
time  expanse  of  1 millisecond,  this  corresponds  to  a real  signal  time  expanse 
of  1 nanosecond.  The  reduced  bandwidth  signals  are  presented  at  the  output 
of  the  radar  as  analog  video  signals  which  then  are  sampled  at  a 13.33  KHz 
rate  which  is  about  six  times  the  nyquist  rate.  Two  hundred  and  fifty-six 
samples,  which  corresponds  to  a time  expanse  of  approximately  19.2  nano- 
seconds, are  stored  in  a buffer.  Each  of  these  256  temporal  samples  is  sub- 
jected to  a threshold  for  detection  purposes.  If  no  detection  occurs,  the 
process  is  terminated  and  the  next  location  is  interrogated  by  the  radar. 

On  the  other  hand,  if  a detection  does  occur,  the  256  samples  are  subjected 
to  the  next  process — discrimination. 
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Figure  151.  Flow  diagram  of  signal  processing  algorithm. 
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In  this  process,  the  location  of  the  target  temporal  signature  must 
first  be  precisely  located.  This  Is  accomplished  by  first  placing  limits 
on  the  signature  gate  position.  In  other  words,  the  gate  must  not  Include 
the  ground  return;  thus,  the  ground  return  Is  first  located  and  the  lower 

boundary  of  the  gate  Is  not  allowed  to  Include  any  of  this  return.  The  j 

gate  Is  located  as  follows : j 

(1)  Search  for  largest  negative  peak  In  the  vicinity  of  the  ground  j 

J 

return  (search  among  30  points) ; ! 

(2)  Search  for  second  zero  crossing  after  the  largest  negative  peak; 

(3)  The  second  zero  crossing  represents  lower  bound  on  the  left  j 

* j 

most  gate  position;  j 

(4)  Upper  bound  on  right  most  gate  position  Is  200  points  beyond 
the  negative  peak  located  In  step  1; 

(5)  Smooth  the  magnitude  of  the  data  between  the  upper  and  lower 
bounds  defined  In  steps  3 and  4 with  a uniform  amplitude  window 
of  5 points; 

(6)  Search  the  smooth  data  for  the  largest  peak; 

(7)  Locate  the  closest  null  In  the  smoothed  array  on  either  side  of 
the  largest  peak  found  In  the  previous  step.  These  nulls  must 
satisfy  the  following  criteria:  (a)  amplitude  of  the  null  must 
be  less  than  or  equal  to  the  maximum  divided  by  4;  (b)  the  first 
null  location  must  be  less  than  or  equal  to  the  maximum  location 
minus  seven  points;  (c)  the  second  null  location  must  be  greater 
than  the  maximum  location  plus  seven  points.  The  gate  extremities 
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are  set  at  the  two  null  locations  which  satisfy  the  criteria 
above . 

(8)  If  the  null  values  do  not  satisfy  the  criteria  outlined  in 
step  7,  set  the  gate  limits  to  the  upper  and  lower  bounds 
determined  in  steps  3 and  4 above. 

After  locating  the  temporal  target  signature  and  gating  out  all  other 
parts  of  the  radar  return,  the  remaining  samples  representing  the  temporal 
signature  are  transformed  to  the  frequency  domain  by  the  use  of  a 64  point 
Fast  Fourier  Transform.  The  number  of  time  samples  is  assumed  to  be  less 
than  32  samples,  or  equivalently,  they  are  assumed  to  occupy  less  than  2.4 
nanoseconds.  Only  the  magnitude  of  the  first  32  points  of  the  frequency 
domain  is  used  in  the  remainder  of  the  processor.  As  indicated  in  the  flow 
diagram  of  Figure  151,  these  frequency  samples  are  denoted  by  (X^;  i = 1,  2,... 
32) . In  a later  discussion  these  32  samples  denoted  by  the  vector  X will 
represent  one  vector  in  a 32-dimensional  space. 

The  nearest  neighbor  (K  = 3)  algorithm  is  used  in  the  discrimination 
processor  since  this  algorithm  yielded  the  best  performance  in  the  data 
analysis.  It  is  assumed  that  a learning  set  of  data  has  been  collected  and 
placed  in  permanent  storage.  Furthermore,  it  is  assumed  that  the  learning 
set  consists  of  126  different  spectral  measurements,  each  consisting  of  32 
frequency  samples  consistent  in  the  frequency  span  with  the  measured  spectral 
data.  That  is,  the  32  samples  represent  frequencies  from  0 to  3.33  GHz. 

This  array  size  was  the  one  used  in  the  data  analysis  previously  discussed 
and  was  assumed  here  for  convenience.  A larger  array  can  be  used  with  possibly 
better  performance  and  also  at  a higher  hardware  complexity  and  cost.  The 
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learning  set  consists  of  a 32  x 126  point  array  whose  members  are  denoted  by 
; i = 1,...32;  j = 1,...126).  Another  way  of  viewing  this  learning  set 
is  that  it  represents  the  location  of  126  vectors  or  points  in  32  dimensional 
space.  Some  of  these  126  points  belong  to  the  mine  category  or  class  and 
some  of  them  belong  to  the  non-mine  category  or  class.  As  indicated  by  the 
flow  diagram,  the  next  step  in  the  discrimination  process  is  to  calculate  the 
distance  between  the  measured  vector  or  point  X and  each  of  the  126  points 
located  in  the  learning  set.  An  example  of  how  the  distance  is  calculated 
for  the  point  in  the  learning  set  is  as  follows : 

32 

k^,  distance  = "S,  (X  - Y ) (50) 

th  1 Ik 

Each  distance  is  associated  with  one  of  132  vectors.  Thus,  each  distance  is 
associated  with  either  the  mine  class  or  the  non-mine  class.  The  discrimina- 
tion decision  is  made  by  selecting  the  three  shortest  distances  which  have 
been  calculated  and  determining  which  of  the  two  classes  are  in  the  majority. 
This  is  equivalent  to  determining  which  three  of  the  learning  set  vectors 
are  closest  to  the  measured  vector.  The  measured  vector  is  assumed  to  be  of 
the  same  class  as  the  majority;  thus,  the  measured  vector  is  designated  as 
a mine  or  a non-mine.  If  it  is  called  a non-mine,  the  processing  is  termi- 
nated and  another  location  is  interrogated  by  the  radar.  On  the  other  hand, 
if  the  measurement  is  decided  to  represent  a mine,  the  classification  process 
I is  Initiated. 

The  classification  process  determines  which  type  of  mine  is  most  likely 
represented  by  the  radar  return  under  investigation.  Since  the  result  of  an 
error  in  mine  type  classification  was  not  deemed  as  disastrous  as  an  error 
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In  discrimination,  a simpler  and  less  complex  classification  algorithm  was 
used  in  this  processing,  namely  the  Fisher's  Linear  Discriminant.  Its  per- 
formance with  the  NBS  data  was  not  quite  as  good  as  that  of  the  nearest 
neighbor  algorithm  but  was  deemed  adequate.  A classification  algorithm 
described  by  the  remainder  of  the  flow  diagram  of  Figure  151  provides  the 
classification  between  the  three  types  of  non-metallic  anti-vehicular  mines. 

Thus,  the  Fisher's  Linear  Discriminant  algorithm  is  applied  twice. 

The  procedure  utilizes  the  32  dimension  vector  or  point  representing 
the  gated  spectrum,  (X^,  i = 1,2... 32)  as  follows: 

(1)  Find  maximum  X^.  ] 

(2)  Locate  nulls  adjacent  to  the  maximum.  \ 

i 

(3)  Define  the  maximum  and  the  two  nulls  as  X^^,  X^,  X^. 

(4)  Form  weighted  sum 

S = Wi  Xi  + W2  X2  + X^  (51) 

The  weighted  sum,  S,  is  subjected  to  a fixed  threshold,  S^.  If  the  threshold 
is  exceeded  by  the  weighted  sum,  a type  1 mine  is  declared  and  the  processing 
is  terminated.  If  the  threshold  is  not  exceeded,  the  classification  weights, 

W^,  W^,  Wj,  used  with  the  Fisher's  Linear  Discriminant  are  changed  to  another 
pre-selected  set  and  the  resulting  sum  is  subjected  to  another  threshold,  S^. 

If  this  threshold  is  exceeded,  a type  3 mine  is  declared.  If  this  threshold 
is  not  exceeded,  a type  2 mine  is  assumed. 

SIGNAL  PROCESSING  EQUIPMENT 

The  signal  processing  equipment  used  in  the  mine  detection  system  must 
be  capable  of  handling  13  mine  detection  radars  and  must  solve  the  processing 
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algorithms  in  111.1  milliseconds.  These  requirements  are  based  on  an  overall 
operational  requirement  that  the  system  cover  a single  look  area  1 foot  long 
by  13.1  feet  wide  with  the  vehicle  moving  at  a speed  of  10  km/hr.  These 
requirements  place  a large  burden  on  the  signal  processing  hardware  both  in 
terms  of  data  handling  requirements  and  processing  time.  The  signal  pro- 
cessing approach  that  will  be  discussed  emphasizes  low  technical  risk  over 
system  cost.  Hardware  trade-offs  will  be  discussed  which  lead  to  a lower 
system  cost  if  the  processing  time  requirements  can  be  met. 

Systems  Configuration 

Figure  152  illustrates  the  suggested  signal  processing  hardware  con- 
figuration. A minicomputer  serves  as  the  central  processor  and  handles  the 
majority  of  the  processing  load.  The  data  from  the  13  radars  could  be  multi- 
plexed into  the  minicomputer  but  this  would  place  a considerable  burder  on 
the  computer.  Instead,  it  is  proposed  that  a preprocessor  be  used  for  each 
radar.  This  reduces  the  data  flow  problem  and  the  processing  time  require- 
ments. In  addition  to  the  minicomputer,  the  system  contains  a cassette  loader 
for  software  loading,  the  computer  control  box,  special  operator  displays, 
and  a high  speed  Fast  Fourier  Transform  module.  All  the  equipment  is  capable 
of  meeting  military  specifications. 

The  mine  detection/classification  algorithms  involve  a three  step 
process:  detection,  discrimination,  and  classification.  The  microprocessors 
would  handle  the  detection  processing.  In  addition,  the  microprocessors 
handle  the  A/D  conversion  processing  and  putting  the  post-detection  data 
into  a format  compatible  with  the  minicomputer  input  requirements.  Data 
would  be  inputted  into  the  minicomputer  in  a Direct  Memory  Access  (DMA)  mode. 
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As  will  be  discussed  later,  the  processing  time  of  the  minicomputer 
Is  too  slow  to  satisfy  the  processing  requirements  associated  with  the  dis- 
crimination and  classification  algorithms.  Therefore,  separate  high  speed 
processing  equipment  would  be  used  to  handle  the  Fast  Fourier  Transform  (FFT) 
algorithms.  This  saves  substantial  processing  time  and  provides  a high 
speed  multiplier  for  use  In  solving  other  processing  algorithms.  There 
are  several  such  high  speed  processors  on  the  market. 

The  minicomputer  Is  controlled  via  a remote  control  box  which  Is  used 
to  start  up  the  system.  The  displays  on  the  control  box  are  not  very  useful 
to  alert  the  operator  to  a detection.  Therefore,  a special  purpose  display 
would  be  used  to  signal  a detection.  The  display  would  consist  of  a set  of 
13  Indicator  lights  and  a buzzer.  The  Indicator  lights  give  a visual  Indi- 
cation of  a detection  and  the  buzzer  gives  an  audible  warning. 

Prior  to  use  of  the  system,  the  software  must  be  loaded  into  the  com- 
puter memory.  A simple  way  to  do  this  Is  to  use  a cassette  loader.  The 
software  is  stored  on  standard  cassette  cartridges  and  loaded  into  the  com- 
puter. An  alternate  approach  is  to  hardwire  the  software  into  the  computer 
memory.  Even  if  the  software  is  hardwired  in,  the  cassette  loader  would 
still  be  desirable  to  load  system  diagnostic  software  for  testing  of  the 
computer.  An  advantage  to  having  the  system  software  on  tape  is  that  changes 
to  the  software  can  be  made  by  simply  replacing  the  cassette. 

The  signal  processing  approach  outlined  in  Figure  152  Is  only  one  possible 
technical  approach.  There  are  many  other  technical  approaches  which  can  be 
used  and  a detailed  trade-off  study  would  have  to  be  conducted  to  select  the 
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best  technical  approach.  Acceptable  system  cost  would  have  to  be  an  input 
to  such  a trade-off  study. 

Minicomputer 

There  are  many  minicomputers  which  could  be  used  in  the  signal  pro- 
cessing hardware.  One  that  Georgia  Tech  has  considerable  experience  with 
is  the  Rolm  1602  processor  which  is  militarized.  The  Rolm  computer  has  a 
1 ps  cycle  time  and  permits  up  to  64K  or  memory.  In  addition,  a DMA  option 
is  available  at  a maximum  input  rate  of  666K  words/second.  Some  of  the  im- 
portant execution  times  are  listed  in  Table  9;  these  will  be  important  in 
the  following  discussions. 


Table  9 

1602  EXECUTION  TIMES 

Operation 

Execution  Time 

Core  Cycle  Time 

1 ps 

Add 

1 ps 

Subtract 

1 ps 

Multiply 

5.4  ps 

Divide 

12.6  ps 

The  1602  processor  is  physically  7.62"  x 10.12"  x 19.56"  and  weighs 
47  pounds  (with  64K  memory) . This  does  not  include  the  control  box  which 
is  7.75"  X 10.25"  X 5.35"  and  weighs  9 pounds.  The  computer  can  be  mounted 
to  the  bulkhead  of  a military  vehicle.  The  power  requirements  of  the  unit 
are  47-400  Hz,  115  vac  power.  The  total  power  requirement  is  up  to  250  watts 
(64K  memory). 
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The  software  for  the  Rolm  computer  operates  under  a Real  Time  Operating 
System  (RTOS)  executive.  The  executive  is  task-oriented  with  a priority 
assigned  to  each  task.  This  permits  multiple  tasks  to  be  run  at  the  same 
time.  The  Rolm  computer  has  extensive  software  support  which  eases  the  soft- 
ware development  problem. 

Preprocessors 

There  are  many  microprocessors  which  can  be  used  for  this  application. 
The  new  Bit  Slice  microprocessors  are  quite  fast  and  provide  a great  deal  of 
flexibility  in  their  architecture.  Georgia  Tech  has  considerable  experience 
with  these  devices  and  has  developed  a standard  microprocessor  for  radar 
applications  around  the  devices.  There  are  other  devices  such  as  the  INTEL 
8080  or  8085  which  can  be  used  and  a trade-off  study  would  have  to  be  con- 
ducted to  select  the  best  microprocessor  for  this  application.  The  Bit  Slice 
microprocessors  are  a good  first  start. 

The  execution  times  for  the  Bit  Slice  microprocessor  are  comparable 
with  the  Rolm  computer  and  can  be  faster.  An  add  takes  approximately  .8  ys 
and  a multiply  takes  3.2  ys  for  a 5 MHz  clock.  In  the  following  analysis, 
the  execution  times  of  the  Rolm  computer  were  used  for  all  the  execution 
time  estimates. 

One  of  the  advantages  of  the  Bit  Slice  microprocessors  is  that  they 
use  a 16-bit  word  which  is  compatible  with  the  Rolm  computer.  Many  of  the 
microprocessors  use  an  8-blt  word  which  requires  reformatting  of  the  data 
words  and  may  force  the  use  of  double  precision  calculations.  There  are 
16-bit  microprocessors  which  can  be  used  in  this  application  as  alternatives 
to  the  Bit  Slice  microprocessors. 
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High  Speed  Processing  Equipment 

F The  signal  processing  algorithms  Involve  solving  an  FFT  which  places 

^ a substantial  burden  on  the  processing  equipment.  Since  FFT's  are  commonly 

encountered  In  signal  processing  applications,  a number  of  special  purpose, 
high  speed  signal  processors  have  become  available  on  the  market.  One  unit 

I 

Is  capable  of  solving  a 64  point  real  FFT  In  .14  milliseconds.  This  Is  a 
considerable  reduction  over  the  time  that  It  would  take  with  the  Rolm  com- 
puter. Since  the  Rolm  computer  operates  on  a task  basis  with  a priority 
assignment,  the  computer  can  send  data  to  the  FFT  module  and  process  other 
data  while  It  Is  waiting  for  a solution. 

Processing  Requirements 

Two  of  the  primary  data  processing  constraints  are  memory  require- 
ments and  processing  time.  The  signal  processor  memory  requirements  are  not 
that  substantial.  The  computer  does  have  to  store  a learning  set  which  con- 
sists of  an  array  of  dimensions  32  x 6 x 126  which  corresponds  to  4032  words 
of  memory.  The  same  learning  set  Is  used  for  each  antenna  return.  A total 
of  32K  memory  should  be  sufficient  to  handle  the  operating  system,  task  soft- 
ware, and  data  storage  requirements. 

The  more  difficult  problem  Is  processing  time.  This  problem  will  be 
addressed  by  looking  at  each  step  of  the  signal  processing  procedure.  First, 
It  Is  helpful  to  estimate  all  the  basic  processing  times.  This  has  already 
been  done  for  the  basic  arithmetic  operations.  The  two  missing  operations 
«re  aquare  root  and  FFT.  The  FFT  operation  can  be  handled  using  special 
■...r (>r >?<*■•  Ini  equipment  **  ha*  already  been  noted.  If  the  FFT  Is  solved 


>ii<  reke  I '1- 1 7 milltercond*  which  becomes  a significant 
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part  of  the  processing.  This,  coupled  with  the  fact  that  an  FFT  only  has 
to  be  solved  when  a detection  is  obtained,  provides  a strong  case  for  the 
processing  configuration  shown  in  Figure  152, 

The  square  root  operation  can  be  handled  in  a number  of  ways . One 
approach  is  to  use  the  approximation: 

2 2 1/? 

(x  + y ) = max  (x,  y)  + 1/2  min  (x,  y)  (52) 

This  involves  a compare,  an  add  and  a shift  operation  which  can  be  com- 
pleted in  3 MS  or  less.  A second  approach  is  to  use  the  Newton-Raphson 
approach  which  is  an  iterative  technique.  Each  iteration  would  require  1 
multiply,  1 divide,  2 adds  and  one  shift  for  a total  of  21  pseconds.  If 
five  iterations  are  performed,  .105  milliseconds  are  required  to  perform 
a square  root  operation.  It  is,  therefore,  advantageous  to  use  the  previous 
approximation  technique  which  requires  only  3 ps. 

Detection  Processing 

The  detection  processing  consists  of  a CFAR  calculation  and  a gating 
procedure.  The  CFAR  calculation  involves  forming  a sum  of  the  form 
256 

S = 2 s (53) 

i=l  ^ 

This  sum  is  multiplied  by  a constant  to  form  a threshold  which  is  compared 
with  all  256  range  bins.  Table  10  summarizes  the  computational  requirements. 

The  total  execution  time  of  the  CFAR  process  takes  approximately  2 milliseconds . 
In  addition  to  the  CFAR  computations,  a time  gate  must  be  put  around 

the  target  return.  The  lower  position  of  the  gate  is  the  second  zero  cross- 

\ 

ing  after  the  largest  negative  peak  in  the  vicinity  of  the  ground  return. 
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Table  10 
CFAR  PROCESSING 


Operation 

Number  of  Operations 

Execution  Time 

Multiply 

257 

1387.8  ys 

Add 

255 

255.0  ys 

Compare 

256 

256.0  ys 

Memory  Cycle 

10 

10.0  ys 

The  upper  end  of  the  gate  is  200  points  from  the  lower  end.  The  data  within 
the  200  point  window  is  smoothed  using  a five  point  uniform  amplitude  average. 
The  largest  peak  in  the  smoothed  data  is  found  along  with  the  closest  nulls 
on  both  sides  of  the  peak.  The  nulls  must  satisfy  the  requirements: 

(1)  Their  amplitude  must  be  less  than  one  quarter  the  peak  value; 

(2)  The  location  of  the  first  null  must  be  less  than  or  equal  to 
the  location  of  the  peak  minus  7 , and 

(3)  The  location  of  the  second  null  must  be  less  than  or  equal  to  the 
location  of  the  peak  plus  seven. 

These  two  nulls  set  the  gate  position.  If  no  such  pair  of  nulls  can  be  found, 
the  previously  computed  gate  positions  are  used. 

The  computational  requirements  of  the  gate  setting  search  procedure 
are  a function  of  the  data.  A worst  case  estimate  is  given  in  Table  11. 

This  gives  a worst  case  computational  time  of  3.75  milliseconds;  the  total 
computational  time  for  the  detection  processing  is  5.66  milliseconds. 

) 

I 


I 
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Table  11 

GATE  SETTING  PROCESSING 

Operation 

Number  of  Operations 

Execution  Time 

Divide 

200 

2520  ps 

Add 

398 

398  PS 

Compare 

831 

831  PS 

Discrimination  Processing 

The  FFT  portion  of  the  discrimination  processing  has  already  been 
discussed.  Once  the  32  point  complex  FFT  has  been  performed  on  the  gated 
data,  the  following  computations  are  performed  between  the  FFT  data  and  a 
learning  set  consisting  of  126  sets  of  32  points: 

(1)  Calculate  the  distance  between  the  FFT  data  and  each  member  of 
the  learning  set; 

(2)  Select  3 largest  distances;  and 

(3)  Determine  which  class  Is  In  the  majority. 

The  class  type  is  associated  with  the  member  of  the  learning  set. 

Table  12  lists  the  computational  requirements  of  the  discrimination 
algorithms.  The  total  discrimination  process  would  take  26.54  milliseconds; 
this  assumes  a high  speed  FFT  module. 

Classification  Processing 

The  classification  process  classifies  the  mines  as  to  1 of  3 types. 
First,  the  maximum  of  the  32  FFT  samples  is  found.  Then  the  minimum  on  either 
side  of  the  maximum  which  is  less  than  1/2  the  maximum  is  found.  The  maximum 
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Table  12 

DISCRIMINATION  PROCESSING 

Operation 

Number  of  Operations 

Execution  Time 

FFT 

1 

14  ps 

Multiply 

4096 

22118.4  PS 

Add 

3938 

3938  ps 

Square  Root 

32 

96  ps 

Compare 

374 

374  ps 

TOTAL 

us 

and  two  minima  form 

a set  of  three  numbers,  x^.  The  weighted  sum 

3 

s - 2 

”i  ’‘i 

(54) 

i-1 

is  formed  where  the 

w^'s  are  weight  factors.  The  sum  is 

compared  to  a 

threshold  which  determines  if  it  is  type  1.  If  the  mine 

Is 

not  type  1,  the 

weighted  sum  Is  recomputed  using  a second  set  of  weights 

and  compared  to  a 

second  threshold. 

This  specifies  whether  it  Is  type  2 or 

3. 

The  computational  requirements  associated  with  the 

classification  pro' 

cess  are  summarized 

In  Table  13.  The  total  computational 

time  required  for 

the  classification 

processing  Is  1.33  milliseconds;  this 

Is 

quite  small  com- 

pared  to  the  discrimination  processing. 

Table  13 

CLASSIFICATION  PROCESSING 

Operation 

Number  of  Operations 

Execution  Time 

Multiply 

7 

37.8  ps 

Add 

4 

4 ps 

Compare 

91 

91  ps 

TOTAL 

132.8  us 

J 
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Total  Processing  Requirements 

As  can  be  seen  from  the  preceedlng  analysis , the  discrimination  pro- 
cessing takes  the  largest  amount  of  processing  time.  This  would  be  even  more 
the  case  if  a high  speed  FFT  processor  is  not  used.  The  total  processing 
time  depends  on  number  of  detections,  number  of  mines,  and  number  of  non- 
mines  encountered  during  any  given  pass  of  the  radar.  If  one  mine  is  detected, 
discriminated  and  classified,  the  total  processing  time  is  32.33  milliseconds. 
This  is  about  a third  of  the  total  time  available;  in  addition,  is  the  execu- 
tion time  of  the  software  executive  and  other  miscellaneous  operations.  For 
the  target  false  alarm  requirements  set  of  the  system,  this  should  provide 
sufficient  time  for  the  worst  case  processing  load.  It  may  be  possible  to 
reduce  the  32.33  milliseconds  by  using  the  high  speed  multiply  in  the  FFT 
module  if  available. 

Alternate  Hardware  Configurations 

It  is  felt  that  the  processing  system  outlined  in  Figure  151  is  the 
Lest  processing  approach  for  the  mine  detection  system  outlined  in  this  report. 
The  primary  trade-off  from  a cost  savings  standpoint  would  be  to  use  a micro- 
processor in  place  of  the  minicomputer.  This  would  save  about  $30,000  in 
hardware  cost  but  would  add  to  software  development  cost  since  the  minicom- 
puter based  system  has  higher  level  software  which  reduces  the  software  devel- 
opment problems.  This  trade-off  would  have  to  be  looked  at  in  more  detail 
during  the  system  development.  The  other  hardware  trade-offs  center  around 
the  selection  of  the  specific  hardware  to  be  used  and  the  exact  specifications 
that  the  hardware  must  meet. 
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signal  Processing 

* At  this  time.  It  Is  difficult  to  go  Into  much  more  detail  on  the  hard- 

• ware  Interfaces  shown  In  Figure  152.  Each  microprocessor  shown  In  Figure  152 
must  have  an  A/D  converter  on  the  Input  and  the  A/D  conversion  process  must 
be  synchronized  with  the  radar  clock  (main  bang) . The  data  transfer  Into 
the  minicomputer  Is  via  the  DMA  mode  and  a Data  Channel  Controlled  would  be 
Included  In  the  minicomputer.  The  minicomputer  would  Interrogate  each  micro- 
processor sequentially  to  check  If  the  microprocessor  has  data  to  transfer 
(l.e.,  a detection).  The  minicomputer  Interfaces  with  the  operator  via  a 
small  display  box  which  consists  of  a set  of  13  LED's  and  a buzzer.  All 

the  other  Interfaces  are  standard  computer  peripheral  Interfaces. 

I It  Is  difficult  at  this  early  stage  of  design  to  estimate  the  system 

I 

j size,  weight,  and  cost.  Table  14  Is  a rough  estimate  of  the  system  break- 

j do%m.  The  total  cost  represents  only  the  hardware  cost.  The  cost  to  proto- 

! type  a unit  Is  estimated  to  be  $200,000  plus  the  hardware  cost.  The  cost  of 

I 

i an  engineering  unit  Is  very  difficult  to  estimate  and  depends  on  the  produc- 

^ tlon  size.  No  attempt,  therefore,  was  made  to  estimate  the  engineering  cost. 


Table  14 

SYSTEM  CHARACTERISTICS 


Unit 

Weight 

Volume 

Cost 

i 

Minicomputer 

47  lbs. 

1.12  ft^ 

$ 45,000 

•t  • 

1 

13  microprocessor 

150  lbs. 

1.6  ft^ 

104,000 

Cassette  loader 

7 lbs . 

.16  ft^ 

5,000 

FFT  processor 

20  lbs. 

.5  ft^ 

20,000 

Display 

5 lbs. 

.12  ft^ 

2,000 

1 , 

r 

Power  Supply 

30  lbs. 

1 ft^ 

5.000 

i' 

1 

TOTALS 

229  lbs. 

$181.000 
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(2)  the  attenuation  and  dielectric  characteristics  of  the  object's  Interior. 

A mathematical  model  including  the  effects  of  these  characteristics  as  well 
as  the  soil  environment  was  derived  and  used  to  verify  this  conjecture.  The 
math  model  yielded  temporal  and  spectral  responses  almost  identical  to  those 
measured  with  the  actual  objects.  Each  of  the  objects  Included  in  the  data 
collection  program  had  different  internal  attenuation  and  dielectric  char- 
acteristics and  some  had  different  thicknesses.  Thus,  different  temporal 
and  spectral  responses  were  obtained.  Based  on  these  different  responses, 

1 

the  objects  were  distinguished  from  one  another.  Even  for  objects  of  the 
same  thickness,  but  of  a different  Interior,  the  objects  were  distinguished 
from  one  another. 

Based  on  the  spectral  characteristics  which  were  different  for  each  j 

buried  object  of  concern,  algorithms  were  designed  to  automatically  operate  , i 

on  the  received  radar  spectral  signatures  to  indicate  if  a mine  was  present  I 

and  to  which  mine  category  it  belonged.  The  process  of  initially  detecting 
a buried  object,  determining  whether  the  detected  object  was  a mine,  then 
Identifying  the  type  of  mine  was  divided  into  three  steps.  The  first  step 
was  to  separate  radar  returns  due  to  potential  targets  from  radar  returns 
due  to  soil  anomalies.  For  this  study,  a potential  target  was  considered  to 
be  any  buried  object  which  caused  a radar  return  that  was  significantly  dif- 
ferent from  radar  returns  caused  by  the  soil  alone.  Some  examples  of  buried 
non-mine  objects  that  might  be  encountered  in  a typical  field  are  rocks, 
roots,  metal  scraps,  etc. 

The  discrimination/classification  steps  which  follow  the  detection 
process  involve  a relatively  complicated  and  time-consuming  process;  thus, 
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not  every  location  in  the  field  can  be  subjected  to  all  three  steps.  The 
detection  processor  was  designed  to  drastically  reduce  the  number  of  loca- 
tions in  a field  that  must  be  scrutinized  by  the  remaining  sections  of  the 
processor.  The  discrimination  processor,  the  second  processing  section, 
was  designed  to  determine  xvhether  the  potential  targets  are  Indeed  targets 
(i.e.,  mines).  The  classification  processor,  the  third  and  final  processing 
section,  was  designed  to  identify  the  class  of  mine  to  which  the  target 
belongs,  choosing  from  several  classes  of  mines. 

Several  algorithms  were  designed  fdr  each  section  of  the  three  pro- 
cessor concept.  The  algorithms  were  tested  using  the  collected  data,  and 
the  performance  of  each  algorithm  was  calculated  and  presented  as  a curve. 

The  measure  of  performance  used  was  the  probability  of  detection  versus  the 
probability  of  false  alarm.  The  performance  of  the  first  processor  in  the 
three  step  processor,  the  detector,  was  ideal  for  dry  soil  (i.e.,  the 
detection  probability  was  1.0,  and  the  false  alarm  probability  was  0).  Since 
the  number  of  data  points  used  was  somewhat  limited,  detection  performance 
that  will  be  obtained  in  a practical  application  may  be  somewhat  different. 
Very  good  performance,  however,  is  highly  likely.  As  the  amount  of  moisture 
was  increased  to  15  or  20%  in  loam  or  clay  soil,  the  detection  performance 
deteriorated  noticeably. 

The  discrimination  processor  separated  non-mines  (i.e.,  rocks,  roots, 
etc.)  from  mines  reasonably  well  in  dry  soil.  For  example,  with  the  best 
algorithm  considered,  the  mine  detection  probability  was  about  0.9  and  the 
false  alarm  probability  (declaring  the  rock,  root,  etc.,  as  a mine)  was  about 

0.1. 
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If  the  area  being  searched  Is  relatively  free  of  rocks,  roots,  or 
metal  scraps,  the  false  alarm  probability  for  the  detector  and  discriminator 
In  cascade  Is  the  product  of  the  Individual  false  alarm  probabilities,  and  a 
discriminator  false  alarm  probability  of  0.1  Is  acceptable.  But  If  there  Is 
a large  number  of  non-mines,  false  alarm  probability  will  be  that  of  the 
discriminator  (e.g.,  0.1).  If  a fair  proportion  of  the  Interrogated  loca- 
tions are  expected  to  Include  non-mines,  one  concludes  from  the  foregoing 
that  the  false  alarm  probability  assigned  to  the  discrimination  processor 
should  be  low.  To  meet  a false  alarm  probability  of  10  ^ for  the  overall 
processor,  the  false  alarm  probability  for  the  discriminator  should  be 
approximately  10  Thus,  for  a large  number  of  non-mine  targets,  the  detec- 
tion processor  algorithms  considered  In  this  study  do  little  In  the  way  of 
reducing  the  data  the  discriminator  has  to  process.  If  most  of  the  locations 
only  contain  soil,  however,  the  detection  processor  significantly  reduces 
the  data  to  be  handled  by  the  discriminator/ classifier,  and  the  overall  pro- 
cessor performance  In  dry  soil  Is  reasonable. 

The  last  stage  of  the  three  step  processor  also  performed  reasonably 
well  In  dry  soil.  The  probability  of  correctly  classifying  a type  1 or  type  2 
mine  was  on  the  order  of  0.9  for  the  best  algorithm  considered. 

The  conclusions  described  above  are  primarily  based  on  the  NBS  col- 
lected data  which  was  much  more  abundant  than  the  MERADCOM  short-pulse  radar 
data.  The  few  temporal  and  spectral  plots  and  performance  curves  made  with 
the  limited  short-pulse  radar  data  were  compared  to  the  NBS  data  plots.  The 
conclusion  drawn  in  comparing  the  NBS  data  with  short-pulse  radar  data  was 
that  the  collection  of  mine  data  by  a GW  measurement  system  such  as  the  NBS 
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system  is  equivalent  to  collecting  mine  data  by  a short-pulse  radar  system. 
Therefore,  all  of  the  conclusions  drawn  from  the  NBS  data  measurement 
program  should  apply  equally  well  to  an  equivalent  set  of  data  taken  with 
the  MERADCOM  short-pulse  radar. 

The  hardware  requirements  for  a mine  detection  system  as  considered 
in  this  report  were  determined.  The  emphasis,  however,  was  on  the  signal 
processing  required  for  short-pulse  returns.  The  radar  to  be  used  in  the 
postulated  mine  detection  system  was  assumed  to  be  the  MERADCOM  short-pulse 
radar  furnished  by  the  government.  A comparison  of  the  processor  perform- 
ance using  the  NBS  data  and  the  short-pulse  radar  data  indicated  approximately 
equivalent  performance.  Therefore,  the  hardware  requirements  for  the  short- 
pulse  radar  as  listed  in  the  CALSPAN  Report  No.  MA-5366-E-3  entitled  "Vehicular 
Mounted  Mine  Detection  Radar  Operational  Manual,"  June  1976,  should  be  suffi- 
cient for  the  mine  detection  system  under  consideration. 

In  recommending  a processor  design,  the  maximum  performance  achieved 
with  the  processing  algorithms  considered  in  this  study  was  deemed  necessary 
to  meet  system  requirements.  This  implies  a target  search  pattern  with  sub- 
terranean interrogations  on  a one-foot  square  grid  (i.e.,  each  interrogation 
is  made  within  one  foot  of  all  other  interrogations) . In  addition,  the 
processing  parameters  used  are  those  assumed  in  the  performance  analysis  pre- 
sented elsewhere  in  this  report. 

The  hardware  requirements  analysis  indicated  that  discrimination  pro- 
cessing takes  the  largest  amount  of  processing  time.  This  would  be  even  more 
the  case  if  the  proposed  high  speed  FFT  processor  were  not  used.  The  total 
processing  time  depends  on  the  number  of  detections,  number  of  mines,  and 
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number  of  non-mines  encountered  during  any  given  pass  of  the  radar.  If  one 
mine  is  detected,  discriminated,  and  classified,  the  tocal  processing  time 
is  approximately  32  milliseconds  for  the  processor  proposed.  This  is  about 
a third  of  the  total  time  available.  Some  of  the  total  time  available,  how- 
ever, must  be  reserved  for  execution  time  of  the  software  executive  and  other 
miscellaneous  operations  which  might  require  another  10  milliseconds.  For 
the  target  false  alarm  requirement  set  for  the  system,  this  should  provide 
sufficient  time  for  the  worst  case  processing  load.  It  may  be  possible  to 
reduce  the  40  milliseconds  by  using  a high  speed  multiply  in  the  FFT  module. 
RECOMMENDATIONS 

A study  of  the  spectral  characteristics  of  the  two  different  sizes  of 
non-metallic  mines  used  in  this  program  suggests  that  only  part  of  the  0.1  GHz 
to  3.0  GHz  frequency  band  is  required  for  successful  discrimination/classifica- 
tion. For  example,  if  only  the  band  from  0.5  to  1.5  GHz  is  used,  much  of 
the  target  signature  is  still  visible  and  the  spectrum  signatures  are  suffi- 
ciently different  to  distinguish  between  the  two  mines.  The  two  inert  mines 
used  are  representative  of  the  mine  sizes  expected.  Thus,  all  other  mines 
of  approxlnuitely  the  same  size  will  have  signatures  in  the  same  spectral 
band,  and  the  signature  characteristics  should  be  sufficiently  different  to 
permit  classification. 

The  benefits  to  be  gained  by  designing  a radar  system  for  this  narrow 
frequency  band  become  clear  when  such  a system  is  compared  to  the  present 
MERADCOM  short-pulse  radar  which  was  designed  for  a much  larger  bandwidth. 

In  addition  to  reducing  the  requirements  on  the  radar  transmitter  and  receiver, 
the  antenna  can  be  significantly  reduced  in  length.  The  elimination  of  the 
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0.12  to  0.5  GHz  band  allows  for  a shorter  antenna  and  a simpler  design.  Thus, 
the  weight  and  bulkiness  of  the  antenna  can  be  reduced  from  that  of  the 
MERADCOM  short-pulse  system.  Of  course,  the  final  size  and  weight  of  the 
antenna  will  depend  on  the  particular  design  selected. 

It  Is  recommended  that  such  an  antenna  be  designed  and  fabricated  for 
Installation  on  the  MERADCOM  short-pulse  radar.  Limited  data  should  then  be 
collected  with  this  new  antenna  design  and  subjected  to  the  algorithms  selected 
and  used  In  this  present  program.  The  performance  of  the  short-pulse  radar 
with  the  newer  and  smaller  antenna  should  then  be  compared  to  the  data  already 
taken  and  analyzed.  This  suggested  follow-on  program  will  allow  for  perform- 
ance verification  of  the  new  and  smaller  radar  mine  detection  system. 
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SPECTRAL  AND  TEMPORAL  PLOTS 
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After  receipt  of  the  data  tapes  at  Georgia  Tech,  the  contents  were 
transferred  to  permanent  disk  storage.  The  data  were  later  converted  from 
ASCII  format  to  more  compact  binary  codes.  Occasionally,  the  pre-burial 
background  data  were  recorded  on  the  same  tape  used  for  target  data,  and 
in  those  cases  the  tape  contents  were  appropriately  split  into  two  separate 
files  to  reduce  access  time.  The  background  runs  are  not  included  in  this 
appendix;  hence,  the  data  file  names  are  not  contiguous. 

Files  GTS,  GT6,  and  GT7  were  used  to  store  the  data  collected  at 
the  NBS  site,  and  GT9  and  GTIO  for  those  collected  at  Gun  Barrel  Hill. 

The  moisture  content  at  the  second  site  varied  somewhat  with  depth;  hence, 
a range  of  values  is  listed.  The  range  given  is  felt  to  be  descriptive 
of  the  soil  directly  over  the  target.  Table  A-1  lists  the  pertinent  para- 
meters associated  with  each  file.  Note  that  the  higher  moisture  content 
(17%  or  18%)  was  not  included  in  the  measurement  matrix  for  the  3-inch 
burial  depth.  Furthermore,  a lapse  of  four  days  occurred  between  the  mea- 
surements of  GT6  and  GT7  during  which  the  targets  were  not  disturbed.  Al- 
though the  test  lane  had  been  covered  with  tarpaulins  against  desiccation, 
heavy  rains  fell  in  the  interim  and  increased  the  moisture  content  by  about 
a percentage  point. 

Each  target  and  its  environs  were  therefore  probed  eight  times  for 
the  various  combinations  of  soil  type,  soil  moisture,  and  burial  depth.  The 
data  are  presented  in  seven  groups  (one  group  per  target) , with  each  group 
consisting  of  the  eight  combinations.  The  locations  and  arrangement  of  the 
data  sets  in  this  report  are  summarized  in  Table  A-2  for  reference. 
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Group 

Data  File 

Depth 

Moisture 

Site 

Soil  Type 

1 

GTS 

3" 

6.8% 

NBS 

Loam 

2 

GTS 

6" 

6.8% 

NBS 

Loam 

3 

GT6 

6" 

17.0% 

NBS 

Loam 

4 

GT7 

6" 

17.8% 

NBS 

Loam 

5 

GT9 

3" 

12%  to  20% 

Gun 

Barrel 

Hill 

Loam 

6 

GT9 

6" 

13%  to  16% 

Gun 

Barrel 

Hill 

Bentonite 

7 

GTIO 

6" 

26%  to  30% 

Gun 

Barrel 

Hill 

Bentonite 

8 

GTIO 

3" 

30% 

Gun 

Barrel 

Hill 

Bentonite 

Table 

A-1.  Data 

file  measurement  conditions. 

File 

GTS 

GTS 

GT6 

GT7 

GT9 

GT9 

GTIO 

GTIO 

Moisture 

7% 

7% 

17% 

18% 

12-20% 

13-16% 

26-30% 

30% 

Depth 

3" 

6" 

6" 

6" 

3" 

6" 

6" 

3" 

Type  A 

Mine 

A-8 

A-9 

A- 10 

A-11 

A-12 

A-1 3 

A-14 

A-IS 

Type  B 

Mine 

A-16 

A-1 7 

A-1 8 

A-19 

A-20 

A-21 

A-22 

A-2  3 

Type  C 

Mine 

A-24 

A-2S 

A-26 

A-27 

A-28 

A-29 

A-30 

A-31 

Plate 

A-32 

A-33 

A- 34 

A-3S 

A-36 

A-37 

A-38 

A-39 

Rock  1 

A-40 

A-41 

A-42 

A-43 

A-44 

A-4S 

A-46 

A-47 

Rock  2 

A-48 

A-49 

A-SO 

A-Sl 

A-S2 

A-S3 

A-S4 

A-SS 

Root 

A-S6 

A-S7 

A-S8 

A-S9 

A-60 

A-61 

A-62 

A-63 

Table  A-2.  Arrangement  of  the  data.  Number  listed  Is  the  page 
on  which  the  data  may  be  found. 
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Table  A-3  lists  additional  information  about  the  data  sets,  including  the 


data  run  numbers  as  numbered  in  Georgia  Tech  numerical  data  files.  j 

The  data  plots  are  presented  on  56  pages,  all  having  the  format  shown 
in  Figure  A-1.  The  plots  on  each  page  are  arranged  in  three  columns  and 
five  rows,  each  row  being  the  data  for  one  of  the  five  positions  over  or 
near  the  target.  The  first  plot  in  each  row  is  the  measured  spectrum;  this 
spectrum  was  weighted  and  transformed  to  the  time  domain,  as  explained  ear- 
lier, to  generate  the  time  domain  response  shown  in  the  second  plot  of  each 
row.  The  data  within  the  range  gate  were  then  transformed  back  to  the  fre- 
quency spectrum  and  unweighted.  The  Georgia  Tech  data  file  is  shown  in  the 
upper  right  corner  and  the  data  set  number  is  printed  along  the  right  hand 
side  of  each  sequence  of  plots. 

Figure  A-2  explains  the  scales  on  the  data  plots.  Both  the  measured 
and  deconvolved  spectra  run  from  DC  to  3048  MHz , and  note  that  the  measured 
spectra  have  been  zeroed  for  the  first  five  frequencies.  The  maximum  ordi- 
nates for  the  measured  spectra  vary  from  plot  to  plot  and  were  typically  in 
the  range  from  0.02  to  0.03.  The  ordinates  for  the  deconvolved  spectra 
also  varied  from  plot  to  plot,  but  were  much  smaller,  typically  of  the  order 
of  0.002.  The  time  domain  plots  span  an  interval  of  20-5/6  nsec  and  were 
normalized  to  a maximum  amplitude  of  unity.  The  range  gate  position  is 
indicated  by  a pair  of  vertical  lines. 
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f! 


142 

-8 

13 

54 

13 

13 

13 

54 

13 

54 

146 

-4 

14 

55 

14 

14 

14 

55 

14 

55 

TYPE  A 150 

0 

15 

56 

15 

15 

15 

56 

15 

56 

MINE  154 

4 

16 

57 

16 

16 

16 

57 

16 

57 

1S8 

8 

17 

58 

17 

17 

17 

58 

17 

58 

322 

-8 

31 

71 

31 

31 

31 

72 

31 

72 

326 

-4 

32 

72 

32 

32 

32 

73 

32 

73 

TYPE  B 330 

0 

33 

73 

33 

33 

33 

74 

33 

74 

MINE  334 

4 

34 

74 

34 

34 

34 

75 

34 

75 

338 

8 

35 

75 

35 

35 

35 

76 

35 

76 

202 

-8 

19 

60 

19 

19 

19 

60 

19 

60 

TYPE  c 
MINE 

206 

-4 

20 

61 

20 

20 

20 

61 

20 

61 

210 

0 

21 

62 

21 

21 

21 

62 

21 

62 

214 

4 

22 

63 

22 

22 

22 

63 

22 

63 

218 

8 

23 

64 

23 

23 

23 

64 

23* 

64 

262 

-8 

25 

65 

25 

25 

25 

66 

25 

66 

266 

-4 

26 

66 

26 

26 

26 

67 

26 

67** 

PLATE 

270 

0 

27 

67 

27 

27 

27 

68 

27 

68 

274 

4 

28 

68 

28 

28 

28 

69 

28 

69 

278 

8 

29 

69 

29 

29 

29 

70 

29 

70 

82 

-8 

7 

48 

7 

7 

7 

48 

7 

48 

86 

-4 

8 

49 

8 

8 

8 

49 

8 

49 

ROCK  1 90 

0 

9 

50 

9 

9 

9 

50 

9 

50 

94 

4 

10 

51 

10 

10 

10 

51 

10 

51 

98 

8 

11 

52 

11 

11 

11 

52 

11 

52 

382 

-8 

37 

77 

37 

37 

37 

78 

37 

78 

386 

-4 

38 

78 

38 

38 

38 

79 

38 

79 

ROCK  2 

390 

0 

39 

79 

39 

39 

39 

80 

39 

80 

394 

4 

40 

80 

40 

40 

40 

81 

40 

81 

398 

8 

41 

81 

41 

41 

41 

82 

41 

82 

22 

-8 

1 

42 

1 

1 

1 

42 

1 

42 

26 

-4 

2 

43 

2 

2 

2 

43 

2 

43 

ROOT 

30 

0 

3 

44 

3 

3 

3 

44 

3 

44 

34 

4 

4 

45 

4 

4 

4 

45 

4 

45 

38 

8 

5 

46 

5 

5 

5 

46 

5 

46 

* Noisy  dais 

**  Contains  a faulty  datum 


Table  A- 3.  Tape  and  run  identication  numbers. 
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All  of  the  sites  for  the  short-pulse  radar  measurements  were  selected 
at  random  in  the  grounds  surrounding  the  Georgia  Tech  Baker  Building.  The 
soil  was  dry  red  clay  and  was  sandy  in  nature.  All  the  mines,  metal  plates, 
and  rocks  were  buried  with  their  top  surfaces  approximately  six  inches  be- 
neath the  soil  surface.  All  large  rocks  and  other  debris  were  removed  during 
the  burying  process.  In  some  cases,  the  surface  grass  was  replaced  for  the 
measurement.  Only  one  of  the  three  radar  modules  was  used  in  the  data  col- 
lection; the  same  radar  was  used  for  all  measurements.  The  radar  antenna 
was  placed  approximately  four  inches  above  the  soil  surface  in  each  case. 

A matrix  of  measurements  made  with  the  short-pulse  radar  is  given  in 
Table  B-1.  Since  the  major  purpose  of  the  short-pulse  measurements  was  to 
verify  that  this  measurement  system  is  effective  for  target  discrimination/ 
classification,  all  of  the  measurements  were  made  in  dry  soil.  All  objects 
were  buried  6 inches  below  the  soil  surface  in  a red  clay  soil.  The  objects 
buried  were  two  types  of  non-met^lic  mines,  a metallic  mine*,  a metal  plate, 
and  a rock  about  the  same  size  as  the  mines.  In  addition,  some  data  were 
taken  in  ground  without  buried  objects  and  is  labeled  in  Table  B-1  as  back- 
ground data.  Each  object  was  burled  in  4 different  locations,  and  for  each 
location,  three  radar  measurements  were  taken.  One  measurement  was  taken 
with  the  radar  antenna  centered  over  the  object,  and  measurements  4 inches  on 
either  side  of  the  center.  Thus,  a total  of  12  measurements  per  object  were 

*The  metallic  mine  measurements  were  taken  to  investigate  the  radar  signatures 
obtained  from  metallic  targets.  This  effort  was  not  in  direct  support  of  the 
non-metallic  mine  detection  program;  however,  these  measurements  are  included 
for  the  sake  of  completeness. 
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Target 

Position 

Location  1 

Location  2 

Location  3 

Location  4 

-4 

1 

19 

37 

55 

Rock 

0 

2 

20 

38 

56 

+4 

3 

21 

39 

57 

-4 

4 

22 

40 

58 

Metal  Plate 

0 

5 

23 

41 

59 

+4 

6 

24 

42 

60 

-4 

7 

25 

43 

61 

Type  A Mine 

0 

8 

26 

44 

62 

+4 

9 

27 

45 

63 

-4 

10 

28 

46 

64 

Type  C Mine 

0 

11 

29 

47 

65 

+4 

12 

30 

48 

66 

-4 

13 

31 

49 

67 

Type  D Mine 

0 

14 

32 

50 

68 

+4 

15 

33 

51 

69 

-4 

16 

34 

52 

70 

Background 

0 

17 

35 

53 

71 

+4 

18 

36 

54 

72 

Table  B-1.  Short-pulse  radar  matrix. 
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taken.  Table  B-1  also  Indicates  the  data  run  numbers  assigned  in  the  Georgia 
Tech  data  files. 

The  format  of  the  data  plots  is  shown  in  Figures  B-1  and  B-2.  The 
plots  on  each  page  are  arranged  in  two  columns  and  three  rows,  each  row 
being  the  data  for  one  of  the  three  positions  over  or  near  the  target.  The 
first  plot  in  each  row  is  the  measured  time  response.  This  data  is  the 
sampled  video  signal  out  of  the  short-pulse  radar.  Most  of  the  ground  return 
has  been  eliminated  by  circuitry  in  the  radar  system.  The  left  most  extremity 
of  the  plot  represents  the  radar  return  just  below  the  ground  surface.  The 
time  scale  on  the  plot  is  0 to  7 nanoseconds  (7  milliseconds  of  video) . The 
two  vertical  lines  mark  the  extent  of  the  target  signature  and  were  positioned 
manually. 

The  second  plot  on  each  row  is  the  spectrum  of  the  gated  time  signa- 
ture obtained  by  use  of  the  Fast  Fourier  Transform.  The  frequency  scale  used 
is  from  0 Hz  to  3.3  GHz.  This  is  the  target  signature  used  for  classification 
and  discrimination.  The  number  in  the  upper  right-hand  corner  of  each  plot 
is  the  value  of  the  vertical  scale  and  is  included  to  accommodate  amplitude 
comparisons.  The  integer  numbers  in  the  right  column  are  the  data  file  num- 
bers listed  in  Table  B-1. 
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NBS  AND  SHORT-PULSE  RADAR 
DISCRIMINATION/CLASSIFICATION 
PERFORMANCE  CURVES 
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Discrimination  is  defined  as  the  process  of  deciding  If  a potential 
target  is  indeed  a target  (mine)  or  some  other  target-like  object  such  as  a 
rock,  root,  piece  of  metal,  etc.  On  the  other  hand,  classification  is  the 
process  of  deciding  to  which  class  of  mines  a particular  target  belongs.  The 
performance  curves  of  both  the  discrimination  and  classification  processors 
considered  in  Section  V are  presented  in  this  Appendix.  First,  a general 
discussion  of  the  format  in  which  the  performance  curves  are  presented  is 
discussed.  Limits  of  performance  are  also  considered  to  give  the  reader 
insight  into  the  measured  results.  Then,  the  performance  curves  for  the  NBS 
data  and  the  MERADCOM  short-pulse  radar  data  are  presented. 

DESCRIPTION  OF  PERFORMANCE  CURVES 


The  curves  in  Figures  C-1,  C-4,  and  C-6  show  the  probability  of  detec- 
tion versus  probability  of  false  alarm  for  a range  of  postulated  discrimina- 
tor processors  whose  performance  range  from  ideal  (perfect)  to  very  poor 


(worst  case).  Here  the  probability  of  detection,  P^,  is  defined  as  the 


probability  of  having  the  discriminator  decide  "mine"  for  a return  resulting 
from  a type  A,  type  B,  or  type  C mine;  the  probability  of  false  alarm,  Pp^, 
is  defined  as  the  probability  of  having  the  discriminator  decide  "mine"  for 
returns  from  non-mine  targets . 

In  order  to  demonstrate  the  meaning  of  the  performance  curve  only  two 
dimensions  will  be  used.  Let  the  mines  and  non-mines  by  represented  by  two 
continuous  density  functions  as  shown  in  Figures  C-1,  C-3,  and  C-5.  Then, 
the  corresponding  probability  curves.  Figures  C-2,  C-4,  and  C-6,  are  gene- 
rated by  sweeping  the  decision  line  from  right  to  left,  and  integrating  the 
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area  under  each  probablllCy  density  function  on  the  right  hand  side  up  to 
S^.  For  widely  separated  density  functions  as  in  Figure  C-1,  the  corres- 
ponding probability  curve  is  the  dashed  line  shown  in  Figure  C-2.  This 


line  is  generated  in  the  following  way:  when  S is  placed  at  point  the 
integrated  area  on  the  right  is  zero  for  both  non-mines  and  mines  so  P 


and  P at  S are  both  zero.  V/hen  S is  moved  to  S , approximately  1/2  of 
r A X 

the  total  area  under  the  density  function  for  mines  will  be  summed  whereas 

the  area  for  the  non-mines  will  still  be  zero.  Hence,  at  S^,  the  P^^  and 

Pp^  become  0.5  and  0 respectively.  At  the  point  the  P^  will  now  be  1.0 

but  the  Ppy^  will  still  be  0.  Hence,  up  to  this  point,  the  dashed  line 

extending  from  (P  , P = 0,0)  to  (P  , P »•  1.0,0)  was  generated.  At  S_ 

D rA  U rA  j 

the  integrated  area  for  the  mines  is  still  the  total  area,  but  now  1/2  of 

the  area  for  the  non-mines  must  also  be  taken  into  account.  Thus,  the  P^ 

and  Pp^  will  have  reached  a value  of  1.0  and  0.5  respectively.  Finally  at 
both  Pp  and  Pp^  will  be  equal  to  1.0.  If  the  density  functions  completely 
overlap  as  in  Figure  C-3,  the  resulting  probability  plot  is  the  dashed 
straight  line  P^  = Pp^  as  shown  in  Figure  C-4.  Other  distributions  as 
illustrated  in  Figure  C-5  result  in  performances  which  fall  into  the  shaded 
region  as  in  Figure  C-6. 

The  performance  curves  for  a classification  processor  appear  somewhat 
different  from  that  of  a discrimination  processor  since  one  is  interested  in 
the  probability  of  properly  detecting  one  class  of  targets  (mines)  versus 
another.  Again,  decisions  must  be  made  in  n-dimensional  space;  however, 
for  simplicity  of  illustration,  the  present  discussion  is  limited  to  one- 
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dimensional  space.  The  probability  of  class  1 versus  probability  of  class  2 

• performance  curves  are  generated  In  much  the  same  manner  as  the  dlscrlmina- 

♦ tion  processor  performance  curves. 

The  decision  line  Is  now  swept  from  left  to  right  and  only  the  area 
^ under  the  class  1 density  function  on  the  right  side  of  the  decision  line 

is  Integrated.  The  area  under  the  class  2 function  is  computed  on  the  left 
side  of  S^.  The  resulting  performances  are  shown  in  Figures  C-8,  C-10,  and 
C-12.  To  indicate  the  correlation  between  the  distributions  of  density 
functions  and  the  probability  curves , the  probabilities  corresponding  to 
the  decision  line  positions  through  were  marked  in  Figures  C-7  and 
C-8.  Thus,  if  the  decision  line  S is  moved  to  the  for  class  1 will 
be  1.0,  but  the  for  class  2 will  be  0,  because  the  area  under  the  class  2 
^ density  function  on  the  left  hand  side  of  is  zero.  At  S2,  the  P^  for 

class  1 will  be  1.0  again,  but  the  P^  for  class  2 will  be  equal  to  0.5,  etc. 
This  performance  curve  is  for  a perfectly  performing  classifier. 

The  dashed  straight  line  shown  in  Figure  C-10  corresponds  to  completely 
overlapping  distributions  as  shown  in  Figure  C-9.  Distributions  which  par- 
tially overlap  as  illustrated  in  Figure  C-11  give  probabilities  which  fall 
into  the  shaded  region  shown  in  Figure  C-12. 

The  foregoing  discussion  indicates  the  nature  and  form  of  the  per- 
formance curves  which  are  to  be  expected  when  actual  measured  data  is  used 
in  the  probability  calculations.  It  must  be  emphasized  that  the  actual  sam- 
ple distributions  were  not  approximated  with  continuous  density  functions. 

As  a result,  the  plots  consist  of  discrete  probability  points.  In  other 
words,  instead  of  integrating  the  area  under  a density  function,  the  individual 
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samples  for  a particular  class  to  the  right  or  left  of  the  decision  line  were 


counted  and  the  sum  divided  by  the  total  number  of  samples  within  the  respec- 


tive class.  Whenever  feasible  a quadratic  polynomial  was  fitted  to  the 


probability  points. 


Because  of  the  limited  amount  of  data,  the  collection  of  signatures 


was  used  not  only  for  a learning  set  but  also  for  a test  set.  For  this 


reason,  two  methods  were  considered  to  estimate  the  performance  of  the 


learning  set.  The  FLD  and  SP  classifiers  utilized  only  those  samples  that 


were  used  to  design  the  classifier,  resulting  perhaps  In  optimistic  estimates 


for  Pp  and  P^,^  (l.e.,  the  estimate  for  P^^  may  be  too  high  and  the  estimate 


for  P„.  may  be  too  low).  The  "leave  one  out"  method  was  applied  to  the 
FA 


nearest  neighbor  classifier.  It  uses  n-1  of  n available  data  samples  to 


design  a classifier,  then  classifies  the  left-out  sample.  This  process  Is 


repeated  n times  to  give  estimates  for  P^  and  This  method  generally 


gives  estimates  which  are  pessimistic  because  not  all  of  the  data  samples 


are  used  to  design  a given  classifier.  Both  estimates  tend  to  converge 


toward  the  true  probabilities  as  the  number  of  samples  In  the  learning  sets 


Increase. 


NBS  DATA 


the  NBS  data  performance  plots  are  Indexed  In  the  matrix  outline 


shown  In  Table  C-1  to  facilitate  comparison  of  performance  between  different 


classifiers,  features,  mines,  and  non-mines.  The  numbers  within  the  matrix 


refer  to  the  figure  where  the  particular  plot  can  be  found.  The  performance 


plots  are  numbered  Figures  C-13  through  C-98.  The  vertical  column  on  the 


left  hand  side  Indicates  the  type  of  classifier  and  the  corresponding 
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Non-Mine  includes  rocks,  root  and  plate 


UNCLASSIFIED 


algorithms  that  were  used.  The  top  row  shows  the  classes  that  were  used 
to  generate  the  performance  curves. 

A typical  descriptive  title  block  for  each  performance  curve  is  shown 
in  Table  C-2.  At  the  top  of  each  figure  is  a classifier  descriptor  which 
corresponds  to  one  of  the  descriptors  listed  in  Table  C-1.  Next  is  a 
list  of  the  ojbects  included  in  each  of  the  two  classes  considered.  The 
last  three  lines  in  the  title  block  give  information  on  the  object  and  the 
object  environment.  For  example,  the  title  block  shown  in  Table  C-2  indi- 
cates that  measurements  from  targets  and  non-mines  burled  at  both  3 inches 
and  6 inches  below  the  ground  surface  were  used  as  the  two  class  data.  Mea- 
surements taken  4 inches  on  either  side  of  the  object  center  as  well  as 
directly  over  the  object  were  also  used.  In  addition,  measurements  of 
objects  for  a variety  of  surrounding  soil  moistures  were  used:  specifically, 
for  soils  with  moisture  content  of  7%,  17%,  12-20%,  and  13-16%. 

FLD/BAYES  CLASSIFIER 
Intuitive  Data 


Contents  of  the 

/ Non-mine 

Plate 

Rock  (1) 

Rock  (2) 

Root 

Two  Classes 

* ( Target 

Type  C 

Information  on 

(Depth  (") 

3 

6 

Object  and  - 

- » < Location 

(")  +4 

0 

-4 

Environment 

f Moisture 

(%)  7 

17 

12-20 

13-16 

Table  C-2 

A Typical  Descriptive  Title  BLock  for  the  Performance  Curves 


Classifier  Description 
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The  performance  curves  for  the  NBS  data  are  followed  by  performance 
plots  for  the  short-pulse  radar  data  that  begins  on  Page  C-99. 

SHORT-PULSE  RADAR  DATA 

The  motivation  for  collection  and  processing  short-pulse  radar  data 
was  to  verify  the  previous  computerized  analysis  which  was  based  on  CM  mea- 
surements over  a broad  range  of  frequencies  (NBS  data) . By  virtue  of  the 
wide  spectrum  of  the  transmitted  pulse,  the  short-pulse  radar  system  repre- 
sents the  analog  hardware  implementation  of  the  FFT  technique  used  to  syn- 
thesize time  domain  data  with  the  NBS  data.  The  short-pulse  radar  temporal 
data  was  processed  by  the  same  discrimination/classification  algorithms  as 
the  NBS  data  and  performance  curves  were  generated  for  comparison  with  the 
NBS  data  performance  plots. 

The  results  of  the  discriminator /classifier  performance  are  organized 
and  displayed  in  matrix  form  as  shown  in  Table  C-3.  Each  entry  in  the 
matrix  refers  to  the  figure  number  in  this  Appendix  where  the  performance  plot 
can  be  found. 

As  will  be  noted,  this  matrix  does  not  cover  discriminator/classifier 
performances  to  the  same  extent  as  the  previous  NBS  data  matrix.  A compari- 
son of  the  NBS  and  short-pulse  radar  gated  spectral  responses  indicated  very 
close  agreement  in  spectral  signatures  for  the  same  type  of  mine.  This  fact 
suggested  that  the  performance  associated  with  the  NBS  and  the  short-pulse 
radar  data  might  be  quite  similar.  Consequently,  to  conserve  our  resources, 
we  elected  to  use  only  the  NN  classification  algorithm  which  had  previously 
given  the  best  performance. 

(continued  on  Page  C-97) 
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Figure  C-13.  FLD/Bayes  classifier. 
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Figure  C-15.  FLD/Bayer  classifier. 
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Figure  C-18.  FLD/Bayes  classifier. 

UNCLASSIFIED 


I 


UNCLASSIFIED 

FLD/BRYES  CLRSSIFIER 

INTUITIVE 
VERSION  2 

DATA 

REFERENCE 

341 

NONMINE 

ROOT 

TARGET 

TYPE  B MINE 

DEPTH 

3 6 

LOCATION 

+ 4 0 

-4 

MOISTURE 

7 17 

12-20  13-11 

PROBRBILITY  OF  FOLSE  OLORM  (PFO) 

Figure  C-19.  FLD/Bayer  classifier. 
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Figure  C-20.  FLD/Bayes  classifier. 
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Figure  C-21.  FLD/Bayer  classifier. 
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Figure  C-22.  FLD/Bayes  classifier. 

UNCLASSIFIED 

C-20 


J 


UNCLASSIFIED 


FLD/BRYES  CLASSIFIER 

INTUITIVE  DRTR 
VERSION  2 
reference  344 

I NONMINE  ROCK( 1 ) ROCK( 2 ) 

\ * TRRGET  TYPE  B MINE 

DEPTH  3 6 

LOCRTION  +4  0 -4 

MOISTURE  7 17  12-20  13-16 


PROBABILITY  OF  DETECT 


UNCLASSIFIED 


FLD/BnYES  CLRSSIFIER 

INTUITIVE  DRTn 
VERSION  2 
REFERENCE  346 


NONMINE 

PLATE 

TARGET 

type  c mine 

DEPTH 

3 6 

LOCATION 

+ 4 0 

-4 

MOISTURE 

7 17 

12-20 

o 

o 


Qo 
Q_  CO 


2: 

CD 


cjS 

LlJ  •- 

U_l 

CD 

La_^ 

>-o' 


CDS 
cc  . 

ODo" 

CD 

QC 

Q_ 

o 

o 


A 

s. 


13-16 


*^•00  0.20  0.40  0'.60  o'.80 

PROBABILITY  OF  FALSE  ALARM  (PFA) 


.00 


Figure  C-25.  FLD/Bayes  classifier. 
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Figure  C-27.  FLD/Bayes  classifier. 
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Figure  C-28.  FLD/Bayes  classifier. 
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Figure  C-30.  FLD/Bayes  classifier. 
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Figure  C-31.  FLD/Bayes  classifier. 
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Figure  C-32.  FLD/Bayes  classifier. 
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Figure  C-37.  Nearest  neighbor  classifier. 
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Figure  C-38.  Nearest  neighbor  classifier. 
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Figure  C-40.  Nearest  neighbor  classifier. 
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Figure  C-41.  Nearest  neighbor  classifier. 
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Figure  C-42.  Nearest  neighbor  classifier. 
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Figure  C-43.  Nearest  neighbor  classifier. 
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Figure  C-A4.  Nearest  neighbor  classifier. 
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Figure  C-45.  Nearest  neighbor  classifier. 
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Figure  C-46.  Nearest  neighbor  classifier. 
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Figure  C-47.  Nearest  neighbor  classifier. 
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Figure  C-50.  Nearest  neighbor  classifier. 
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Figure  C-51.  Nearest  neighbor  classifier. 
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Figure  C-53.  Nearest  neighbor  classifier. 
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Figure  C-54 . Nearest  neighbor  classifier. 
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Figure  C-55.  Nearest  neighbor  classifier. 
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Figure  C-56.  Nearest  neighbor  classifier. 
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Figure  C-57.  Nearest  neighbor  classifier. 
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Figure  C-59.  Nearest  neighbor  classifier. 
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Figure  C-60.  Nearest  neighbor  classifier. 
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Figure  C-61.  Nearest  neighbor  classifier. 
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Figure  C-62.  Nearest  neighbor  classifier. 
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Figure  C-63.  Nearest  neighbor  classifier. 
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Figure  C-65.  Nearest  neighbor  classifier. 
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Figure  C-66.  Nearest  neighbor  classifier. 
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Figure  C-67.  Nearest  neighbor  classifier. 
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Figure  C-68.  Nearest  neighbor  classifier. 
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Figure  C-69.  Nearest  neighbor  classifier. 
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Figure  C-73.  Nearest  neighbor  classifier. 
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Figure  C-74.  Nearest  neighbor  classifier. 
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Figure  C-75.  Nearest  neighbor  classifier. 
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Figure  C-76.  Nearest  neighbor  classifier. 
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Figure  C-77.  Nearest  neighbor  classifier. 
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Figure  C-78.  Nearest  neighbor  classifier. 
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Figure  C-79.  Nearest  neighbor  classifier. 
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Figure  C-80.  Nearest  neighbor  classifier. 
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Figure  C-81.  Nearest  neighbor  classifier. 
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Figure  C-82.  Nearest  neighbor  classifier. 
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Figure  C-84.  Nearest  neighbor  classifier. 
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Figure  C-87.  Cross-correlation  classifier. 
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Figure  C-94.  Cross-correlation  classifier. 
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Figure  C-95.  Cross-correlation  classifier. 
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Figure  C-96.  Cross-correlation  classifier. 
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Also,  short-pulse  radar  data  was  only  taken  in  dry  soil  comparable  to 
the  7%  moisture  soil  used  in  the  NBS  data  collection.  In  addition,  a limited 
number  of  non-mine  objects  were  included  in  the  data  collection — e.g.,  a rock 
and  a metal  plate  burled  six  Inches  in  dry  soil.  Soil-only  data  where  no 
mine-like  objects  were  present  were  also  used  to  exercise  the  classification 
algorithm.  The  various  combinations  of  mines  and  non-mines  analyzed  is  shown 
in  Table  C-3,  Short-Pulse  Radar  Performance  Matrix. 

It  should  also  be  noted  that  the  MERADCOM  short-pulse  radar  data  con- 
tained fewer  measurements  or  signatures  that  the  NBS  data.  For  example,  the 
short-pulse  data  consisted  of  12  type  C mine  signatures  as  compared  to  18 
for  the  NBS  data.  Consequently,  the  statistical  confidence  level  for  the 
short-pulse  data  is  not  quite  as  good  as  for  the  NBS  data.  Nevertheless, 
the  short-pulse  data  does  allow  a relative  comparison  of  performance  for  the 
two  measurements  systems . 

The  performance  curves  are  numbered  Figures  99  through  128. 
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Figure  C-101.  Nearest  neighbor  classifier. 
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Figure  0-104.  Nearest  neighbor  classifier. 
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Figure  C-106.  Nearest  neighbor  classifier. 
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Figure  C-109.  Nearest  neighbor  classifier. 
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Figure  C-110.  Nearest  neighbor  classifier. 
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Figure  C-111.  Nearest  neighbor  classifier. 
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Figure  C-112.  Nearest  neighbor  classifier. 
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Figure  C-116.  Nearest  neighbor  classifier. 
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Figure  C-117.  Nearest  neighbor  classifier. 
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Figure  C-119.  Nearest  neighbor  classifier. 
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Figure  C-120.  Nearest  neighbor  classifier. 
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Figure  C-124.  Nearest  neighbor  classifier. 
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Figure  C-125.  Nearest  neighbor  classifier. 
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Figure  C-126.  Nearest  neighbor  classifier. 
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Figure  C-127.  Nearest  neighbor  classifier. 
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Figure  C-128.  Nearest  neighbor  classifier. 
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A process  called  homomorphic  deconvolution  was  Investigated  as  part 
of  the  study  for  applicability  to  the  mine  detection  problem.  This  approach 
theoretically  allows  one  to  trade  radar  bandwidth  for  processing  complexity. 

In  other  words,  additional  processing  might  synthetically  provide  greater 
range  resolution  than  one  could  conventionally  achieve  with  a given  radar 
system.  The  investigation  of  this  technique  was  a low  priority  task  and 
was  not  completed;  however,  some  preliminary  conclusions  were  reached  based 
on  the  work  that  was  accomplished.  These  conclusions  and  a description  of 
the  process  are  given  in  this  Appendix. 

Functionally,  homomorphic  deconvolution  accomplishes  the  same  task  as 
inverse  filtering — it  accentuates  the  extremeties  of  the  spectrum  to  yield  a 
larger  3-dB  bandwidth.  Inverse  filtering  is  illustrated  in  Figure  D-1.  The 
spectra  of  the  signal  at  the  input  and  output  of  the  inverse  filter  is  illus- 
trated in  Figures  D-l(a)  and  (d) . The  inverse  filter  is  designed  to  have  a 
frequency  response  which  is  the  inverse  of  the  fundamental  component  of  the 
input  spectrum  as  illustrated  in  Figures  D-l(b)  and  (c) . In  this  filtering 
process,  it  is  assumed  that  the  information  of  interest  is  the  smaller,  faster 
varying  component  of  the  input  spectrum  which  is  preserved  in  the  inverse 
filtering  process.  In  a radar  application  where  a transmitted  signal  is 
reflected  from  multiple  surfaces  or  points  extended  in  range  (or  depth  for 
a mine  detection  radar) , the  slowly  varying  spectra  component  represents  the 
transmitted  signal  spectra  while  the  faster,  more  rapidly  varying  component 
yields  information  about  the  separation  between  scatterers.  It  is  this  latter 
information  which  is  of  interest  and  which  is  preserved  in  the  inverse  filter- 
ing process. 
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Figure  D-1.  Inverse  filtering:  (a)  spectrum  before  filtering;  (b)  an  estimate  of  the 
fundamental  component  of  the  spectrum;  (c)  the  inverse  filter  response; 
and  (d)  the  spectrum  at  the  output  of  the  inverse  filter. 
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One  of  the  problems  in  using  inverse  filtering  is  that  the  portion  of 
the  input  spectrum  which  represents  the  transmitted  signal  must  be  known 
precisely.  Including  magnitude,  in  order  to  perform  inverse  filtering.  The 
output  of  this  type  of  filter  is  very  sensitive  to  filter  mismatch.  Thus, 
in  practice,  the  theoretical  performance  of  Inverse  filtering  is  seldom 
achieved . 

One  of  the  proported  advantages  of  homomorphic  deconvolution  is  that 
the  transmitted  signal  does  not  need  to  be  known.  The  process  separates  the 
radar  return  spectrum  into  two  parts — that  part  representing  the  transmitted 
slngal  and  that  portion  which  yields  information  about  scatterer  position. 

The  equivalent  operation  in  the  time  domain  is  called  deconvolution.  Radar 
returns  can  be  mathematically  represented  by  the  convolution  of  the  trans- 
mitted temporal  signal  with  a series  of  point  or  plane  scatterers  repre- 
sented by  impulse  functions.  If  we  can  successfully  deconvolve  the  trans- 
mitted signal  and  the  impulse  functions,  then  the  range  resolution  has  been 
enhanced  and  is  theoretically  zero. 

The  non-linear  process  called  homomorphic  deconvolution  is  illustrated 
in  Figure  D-2  where  the  transmitted  signal  is  denoted  by  x(t)  in  the  time 
domain  and  XCe^'^)  in  the  frequency  domain.  The  response  h(t)  is  the  Impulse 
function  representing  the  point  scatterer  distribution  in  time  or  range.  The 
radar  return  is  denoted  as  y(t)  and  Y(e^'^),  respectively.  As  indicated  in 
Figure  D-2,  the  process  involves  the  use  of  the  forward  and  inverse  Fourier 
transform  F and  F respectively.  The  complex  natural  logarithm  and  inverse 
natural  logarithm  (exponential)  is  also  utilized.  The  time  domain  function 
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Figure  D-2.  The  homomorphic  deconvolution. 
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formed  by  the  logarithmic  process  Is  called  the  cepstrum.  Several  types  of 

\ 

linear  filtering  can  be  used  in  the  cepsturm  domain  to  separate  x(t)  and 
h(t),  which  are  non-linear  versions  of  x(t)  and  x(t).  The  filtering  techni- 
que often  used  is  to  simply  low-pass  and  high-pass  filter;  that  is,  dividing 

y(t)  at  some  convenient  value  of  t as  shown  in  Figure  D-2.  Each  of  these 
\ 

components  of  y(t)  are  then  transformed  back  into  the  frequency  and  time 

domains  by  an  inverse  non-linear  process.  The  result  is  a deconvolution 
A A, 

if  h(t)  and  x(t)  are  separable  in  the  cepstrum  domain.  This  depends  some- 
what on  the  separation  and  relative  amplitude  of  the  scatterers . The  decon- 
volved components  x(t)  and  y(t)  are  shown  in  Figure  D-2. 

In  practice,  the  deconvolution  process  just  described  is  not  as 
simple  as  it  sounds.  One  must  be  careful  to  insure  that  the  logarithmic 
transformation  provides  a unique  transformation.  A process  not  indicated  in 
Figure  D-2,  called  "phase  unwrapping"  must  be  performed.  The  complex  log 
function  is  module  2 it  in  its  imaginary  component  which  component  also  repre- 
sents the  phase  variation  of  the  frequency  spectrum.  Thus,  the  imaginary  com- 
ponent must  be  "unwrapped"  so  that  is  is  not  module  2it.  In  addition  to  this 
special  consideration,  one  must  also  remove  the  linear  component  of  the  spec- 
trum phase  before  transforming  into  the  cepstrum  domain.  This  linear  phase 
must  be  replaced  after  the  divided  cepstruins  are  transformed  back  into  the 
frequency  domain.  A flow  diagram  which  itemizes  the  functions  to  be  performed 
in  the  homomorphic  deconvolution  process  is  shown  in  Figure  D-3.  In  this 
diagram,  only  the  high-passed  portion  of  the  cepstrum  is  processed  since  only 
the  scatterer  distribution  is  of  interest.  Detailed  discussions  concerning 
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the  practical  considerations  of  homomorphic  deconvolution  are  given  in  Refer- 
ences 7 through  23. 

A computer  program  to  perform  the  homomorphic  deconvolution  described 
by  Figure  D-3  was  written  and  exercised.  It  was  Initially  Intended  to  test 
the  program  with  synthetically  derived  radar  returns  and  then  to  use  the 
program  on  measured  NBS  data.  The  program  was  debugged  and  tested  on  syn- 
thetic data;  however,  the  mine  detection  program  terminated  before  measured 
data  could  be  utilized. 

Some  results  of  the  deconvolution  program  on  synthetic  data  are  shown 
in  Figures  D-4  and  D-5.  In  Figure  D-4(a),  it  is  assumed  that  two  unequal 
size  scatterers  separated  in  range  reflect  a transmitted  pulse  similar  to 
the  NBS  short  pulse.  The  second  return  is  one-half  the  amplitude  of  the  first 
return.  The  spectrum  of  this  composite  radar  return  is  shown  in  Figure  D-4(b). 
The  period  of  the  rapidly  varying  component  of  this  spectrum  is  inversely 
proportional  to  the  separation  of  the  two  scatterers.  The  slowly  varying 
component  of  the  spectrum  is  the  spectrum  of  the  transmitted  waveform.  The 
upper  portion  of  the  cepstrum  is  shown  in  Figure  D-4(c);  the  lower  portion 
has  been  eliminated.  The  deconvolved  return  is  shown  in  Figure  D-4(d). 

Note  that  the  time  spread  (resolution)  of  the  two  returns  is  smaller  than 
in  Figure  D-4 (a).  Thus,  the  resolution  has  been  Increased  via  this  process. 

Figure  D-5  depicts  a similar  situation  except  the  second  radar  return 
is  now  almost  as  large  as  the  first.  The  point  to  be  made  in  this  example 
is  illustrated  in  Figure  D-5(c).  Note  that  a decaying  train  of  impulses 
occurs  in  the  cepstrum  domain.  The  rate  of  decay  depends  on  the  relative 
size  of  the  two  radar  returns.  If  the  two  returns  are  equal  in  amplitude. 
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Figure  D-4.  Results  of  deconvolution  on  synthetic  data  where  the  second  scatterer  is 
one-half  the  size  of  the  first. 
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Figure  D-5.  Results  of  deconvolution  on  synthetic  data  where  the  second  scatterer  is 
95%  of  the  size  of  the  first. 
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the  Impulse  train  In  the  cepstrum  domain  does  not  decay  but  consists  of  an 
Infinite  number  of  Impulses  of  equal  amplitude.  If  the  second  r^ular  return 
Is  larger  than  the  first,  the  Impulse  train  In  the  cepstrum  domain  diverges 
without  limit.  Since,  In  practice,  the  second  return  can  easily  be  larger 
than  the  first,  one  must  consider  ways  of  solving  this  processing  problem. 

One  approach  Is  to  weight  the  radar  return  response  to  force  the  second  re- 
turn to  be  smaller  than  the  first  and  then  later  remove  the  weighting  factor. 

The  foregoing  example  considered  was  for  a simple  two  scatterer  prob- 
lem. If  more  scatterers  are  considered,  the  process  could  become  more  com- 
plicated In  a practical  application. 


